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Single-bubble sonoluminescence occurs when an acoustically trapped and periodically driven gas 
bubble collapses so strongly that the energy focusing at collapse leads to light emission. Detailed 
experiments have demonstrated the unique properties of this system: the spectrum of the emitted light 
tends to peak in the ultraviolet and depends strongly on the type of gas dissolved in the liquid; small 
amounts of trace noble gases or other impurities can dramatically change the amount of light 
emission, which is also affected by small changes in other operating parameters (mainly forcing 
pressure, dissolved gas concentration, and liquid temperature). This article reviews experimental and 
theoretical efforts to understand this phenomenon. The currently available information favors a 
description of sonoluminescence caused by adiabatic heating of the bubble at collapse, leading to 
partial ionization of the gas inside the bubble and to thermal emission such as bremsstrahlung. After 
a brief historical review, the authors survey the major areas of research: Section II describes the 
classical theory of bubble dynamics, as developed by Rayleigh, Plesset, Prosperetti, and others, while 
Sec, III describes research on the gas dynamics inside the bubble. Shock waves inside the bubble do 
not seem to play a prominent role in the process. Section IV discusses the hydrodynamic and chemical 
stability of the bubble. Stable single-bubble sonoluminescence requires that the bubble be shape 
stable and diffusively stable, and, together with an energy focusing condition, this fixes the parameter 
space where light emission occurs. Section V describes experiments and models addressing the origin 
of the light emission. The final section presents an overview of what is known, and outlines some 
directions for future research, 
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1. INTRODUCTION 


A. The discovery of single-bubble sonoluminescence 


Single-bubble sonoluminescence was discovered in 
1989 by Felipe Gaitan, then a graduate student at the 
University of Mississippi working with Larry Crum 
(Gaitan, 1990; Gaitan and Crum, 1990; Gaitan et al., 
1992). Crum had seen hints of light emission from a 
single bubble in 1985 (Crum and Reynolds, 1985), and 
Gaitan's objective for his thesis was to search systemati- 
cally for it. Gaitan was carrying out a set of experiments 
on the oscillation and collapse of bubbles, using a flask 
of liquid lined with transducers tuned to set up an acous- 
tic standing wave at the resonant frequency at the jar. 
When the pressure amplitude P, of the sound waves is 
larger than the ambient pressure Pọ=1 bar, the pressure 
in the flask becomes negative, putting the liquid under 
tension. At large enough tension, the liquid breaks apart. 
(cavitation), creating unstable bubble clouds in which 
the bubbles often self-organize into dendritic structures 
(streamers; see Neppiras, 1980). These cavitation clouds 
collapse with enormous force, powerful enough to do 
serious damage to the surfaces of solid bodies in their 
vicinity. 

In his search for single-bubble sonoluminescence, 
Gaitan at some point found a regime with a moderate 
forcing pressure P,/Po~1.2—1.4 and with the water de- 
gassed to around 20% of its saturated concentration of 
air. He then observed that "as the pressure was in- 
creased, the degassing action of the sound field was re- 
ducing the number of bubbles, causing the cavitation 
streamers to become very thin until only a single bubble 
remained. The remaining bubble was approximately 20 
um in radius and [...] was remarkably stable in position 
and shape, remained constant in size and seemed to be 
pulsating in a purely radial mode. With the room lights 
dimmed, a greenish luminous spot the size of a pinpoint 
could be seen with the unaided eye, near the bubbles 
position in the liquid" (Gaitan et al., 1992). The experi- 
ment is shown in Fig. 1, a sketch of a typical experimen- 
tal setup for single-bubble sonoluminescence in Fig. 2. 

At the time of Gaitan's experiment, all previous work 
with light-emitting bubbles involved many unstable 
bubbles being simultaneously created and destroyed. 
Using Mie scattering to track the volumetric contrac- 
tions and expansions of the bubbles (Gaitan, 1990; 
Gaitan and Crum, 1990; Gaitan et al., 1992) Gaitan and 
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Single-bubble sonoluminescence 





FIG. 1. A sonoluminescing bubble. The dot in the center of 
the jar is the bubble emitting light. From Crum, 1994. 


co-workers demonstrated that their setup indeed gener- 
ated a single bubble, undergoing its oscillations at a 
fixed, stable position at a pressure antinode of the ultra- 
sound field in the flask. The oscillation frequency f is 
that of the sinusoidal driving sound (typically 20-40 
kHz), but the dynamics of the bubble radius is strongly 
nonlinear. Once during each oscillation period, the 
bubble, whose undriven (ambient) radius Ro is typically 
around 5 um, collapses very rapidly from its maximum 
radius R,,—50 um to a minimum radius of Rmin 
~0.5 um, changing its volume by a factor of 1x10" 
(Barber et al., 1992). Figure 3 shows the radius, forcing 
pressure, and light intensity (top to bottom) during this 
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FIG. 2. Sketch of a typical setup for generating sonoluminesc- 
ing bubbles. 
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FIG. 3. Radius R(r), driving pressure P(t), and light intensity 
I(f) from Crum (1994), as measured by Gaitan et al. (1992). A 
negative driving pressure causes the bubble to expand; when 
the driving pressure changes sign, the bubble collapses, result- 
ing in a short pulse of light, marked SL. 


process (Crum, 1994). The bubble expansion caused by 
the negative pressure is followed by a violent collapse, 
during which light is emitted. The process repeats itself 
with extraordinary precision, as demonstrated by mea- 
surements of the phase of the light emission relative to 
the driving. 

Light emission from collapsing ultrasound-driven 
bubbles had long been dubbed sonoluminescence (SL). 
Researchers were familiar with the energy-focusing 
power of cavitation clouds, and it was therefore not sur- 
prising when Frenzel and Schultes (1934) demonstrated 
that these cavitation clouds emitted a low level of light 
[slightly earlier, Marinesco and Trillat (1933) had found 
indirect evidence when photographic plates fogged in an 
ultrasonic bath]. After all, if the cloud collapses violently 
enough to break the molecular bonds in a solid, causing 
cavitation damage (Leighton, 1994), there is no reason 
why photons should not also be emitted. The energy- 
focusing power of the cavitation cloud was understood 
to arise from a singularity occurring when a bubble col- 
lapses in an ambient liquid (Rayleigh, 1917): inertial 
forces combined with mass conservation lead to bubble- 
wall velocities that become supersonic during the col- 
lapse, causing rapid heating of the bubble interior. To 
the engineering community of the time, the fluid me- 
chanics of this process were much more interesting than 
the character of the radiation produced. This was for a 
very practical reason: people wanted to understand how 
to prevent cavitation damage, or how to harness its 
energy-focusing power. Although historically the light 
emission has played a useful role in measuring proper- 
ties of cavitation [Flint and Suslick (1991b) used the 
spectrum to measure the temperature in a cavitating 
bubble cloud], it was not considered of intrinsic impor- 
tance until Gaitan's discovery of what is now known as 
single-bubble sonoluminescence (SBSL). 

The brightness of Gaitan's single, isolated bubble 
caused great excitement in the scientific community; it is 
visible to the naked eye! Though the light emission from 
conventional cavitation clouds [now called multibubble 
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sonoluminescence (MBSL); see Kuttruff, 1962; Walton 
and Reynolds, 1984; Brennen, 1995] is also visible as 
diffuse glowing, in that case no individual, stable 
bubbles can be identified. The excitement about single- 
bubble sonoluminescence was driven in large part by a 
set of experiments by Seth Putterman’s group at UCLA 
from 1991 to 1997, which exposed further peculiarities, 
making single-bubble sonoluminescence seem very dif- 
ferent from MBSL (the experiments of the UCLA 
group are reviewed by Barber et al., 1997 and Putterman 
and Weninger, 2000). Was new physics (beyond that im- 
plied by the collapse mechanism of Lord Rayleigh in 
1917) responsible for this difference? Many people were 
also excited by the fact that single-bubble sonolumines- 
cence appeared to be much more controllable than its 
multibubble counterpart, bringing expectations of both 
good careful scientific studies and the possibility of new 
technologies, including the harnessing of the energy- 
focusing power of SBSL. 

It is natural that the excitement at first caused specu- 
lation about very exotic conditions inside the bubble, 
such as extremely high temperatures and pressures 
Even Hollywood caught on to the excitement, producing. 
a movie in which the central character created a fusion 
reactor using a single sonoluminescing bubble. As the 
field matured over time and the models were refined, 
the results became more down to earth; for instance, the 
commonly believed maximum temperature at the 
bubble collapse has been revised downward during a de- 
cade of research from early estimates of 10* K to the 
more modest present-day estimates which cluster 
around 10* K. 

In the years since SBSL was discovered, much has 
been learned about how and why it occurs. The goal of 
this review is to clarify the basic ideas that have proven 
necessary for a quantitative understanding of single- 
bubble sonoluminescence and to present an overview of 
the current state of the field, of what is known and what 
is yet to be fully understood. 


B. Structure of the review 


The structure of this review is as follows: The remain- 
der of this Introduction presents an overview of the sa- 
lient historical and experimental facts and qualitatively 
describes the ideas and issues that have been shown to 
be important for understanding the phenomenon. This 
overview will illustrate the enormous variety of physical 
processes taking place inside this simple experiment, 
ranging from fluid dynamics, to acoustics, to heat and 
mass transfer, to chemical reactions, and finally to the 
light emission itself. We shall then spend the next four 
sections following in detail the sequence of events that 
happen to a sonoluminescing bubble, beginning with the 
motion of the flask and liquid and proceeding to the 
dynamics of the bubble wall and interior. Figure 4 shows 
the radius R(t) of the bubble as a function of time dur- 
ing a single cycle of the driving; the inset blows up the 
innermost ~60 ns around the cavitation event, where 
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FIG. 4. Classical bubble dynamics calculation for a driving 
pressure amplitude P, 12 atm, frequency f=26.5 kHz, and 
ambient bubble radius Ry=4.5 um. One oscillation cycle of 
R(t) is shown. The bubble expands to nearly ten times its 
ambient radius, then collapses extremely quickly, leading to 
adiabatic heating of the gas inside the bubble. The collapse is 
followed by afterbounces with roughly the eigenfrequency of 
the bubble. The vertical dashed lines and small-print numbers. 
indicate the intervals 1-10 (summarized in Sec. VI) at which 
different physical processes are important, which are discussed 
throughout the review. The inset shows the innermost 60ns 
around the time ¢* of maximum compression and highlights 
the bubble radius during Rayleigh cavitation collapse, where 
the light is emitted. 


the bubble temperature rises rapidly due to adiabatic 
compression and light is emitted. 

Section II reviews classical studies of the hydrody- 
namics of bubble motion, showing, for example, how to 
derive the equation for the bubble radius leading to Fig. 
4, and also discussing the hydrodynamics of Lord Ray- 
leigh's cavitation collapse (Fig. 4, inset). Section III de- 
scribes the fluid dynamics of the bubble’s interior, focus- 
ing mainly on what happens to the gas during the 
cavitation event, but also discussing water evaporation, 
heat transfer, and chemistry. Section IV discusses the 
physical processes that fix the ambient size Ry of the 
bubble, including the diffusive and chemical processes of 
mass exchange between bubble and liquid as well as me- 
chanical stability constraints. Finally, Sec. V discusses 
the light emission itself, which occurs when the bubble is 
in its maximally compressed state. The discussion will 
emphasize the mechanisms that are consistent with the 
current experimental data. In the final section, we give a 
brief summary and present our opinions on the current 
state of the field as well as the areas of activity with the 
brightest outlook for future work. 


C. Historical overview 


After Gaitan's discovery, the initial goal of research 
was to quantify how much more efficiently a single 
bubble focuses energy than a bubble cloud. To address 
this question, Barber and co-workers (Barber and 
Putterman, 1991; Barber etal, 1992) measured the 
width of the light pulse, by studying the response of a 


Rev. Mod. Phys., Vol. 74, No. 2, April 2002 





300 500 700 
‘wavelength (nm) 


FIG. 5. Spectrum of single-bubble sonoluminescence, for wa- 
ter at 22°C. The data points are redrawn from Fig. 1 of Hiller 
et al. (1992). Fits to a blackbody spectrum can be attempted 
for different temperatures, with best results for about 40 000 K 
(solid line), higher than the 25000 K suggested by Hiller et al. 
(1992). 


single photomultiplier tube to the sonoluminescent flash. 
It was concluded that the width of the light pulse was 
less than 50 ps. The importance of the measurement was 
that this upper bound for the pulse width was much 
smaller than the time during which the bubble remained 
in its most compressed state. Roughly, the time scale of 
bubble compression is given by the time it takes a sound 
wave to travel across the minimum radius of the bubble. 
With a sound velocity of c~1000 m/s, one obtains a 
ballpark estimate of Rymin/c~10~° s, far in excess of the 
measured pulse-width limit, Lord Rayleigh’s cavitation 
mechanism implies that the energy focusing is coupled 
to the bubble collapse: this discrepancy suggests that in 
SBSL the light emission is decoupled from the bubble 
dynamics. 

The gauntlet was thus thrown, and a search for the 
correct mechanism began. An influential early idea [in- 
troduced independently by Greenspan and Nadim 
(1993), Wu and Roberts (1993), and Moss et al. (1994)] 
was that the energy focusing in the bubble was caused 
by a converging spherical shock. It had been known 
since the seminal work of Guderley (1942) (see also 
Landau and Lifshitz, 1987) that such shocks focus en- 
ergy, and in the absence of dissipation the temperature 
of the gas diverges to infinity. In fact, Jarman (1960) had 
already suggested converging shocks as the source of 
multibubble sonoluminescence. This mechanism neatly 
solved the upper-bound problem for the width of the 
light pulse (since in this picture the light originates from 
a much smaller region in the center of the bubble) and 
proposed an elegant mechanism for energy focusing 
compounding Lord Rayleigh bubble-collapse mecha- 
nism. Simulations by Wu and Roberts (1993) had the 
maximum temperature approaching 105 K, very hot in- 
deed. 

For several years, experimental information accumu- 
lated about the properties of sonoluminescing bubbles. 
Hiller et al. (1992, 1994, 1998) measured the spectrum of 
a sonoluminescing air bubble in water and demonstrated 
that it increases toward the ultraviolet (Fig. 5). The ap- 
parent peak in some spectra is due to the strong absorp- 
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FIG. 6. MBSL (thin line) and SBSL (thick line) spectra in a 
0.1M sodium chloride solution. Each spectrum was normalized 
to its highest intensity. Note the prominence (MBSL) and ab- 
sence (SBSL, see the inset for an enlargement) of the sodium 
line near 589 nm. Figure reproduced from Matula et al. (1995). 


tion of wavelengths below ~200 nm by the water in the 
flask. In sharp contrast to the spectrum of MBSL, single- 
bubble sonoluminescence shows a smooth continuum, 
without spectral lines (see Fig. 6). The presence of spec- 
tral lines points to lower temperatures, since the atomic 
transitions leading to lines tend to be overwhelmed by 
continuous emission processes at high temperatures. By 
fitting the observed spectra to that of a blackbody emit- 
ter (Fig. 5), Hiller et al. (1992) concluded that the tem- 
perature of the gas was at least 25000 K. 

Barber et al. (1994) demonstrated that both the light 
intensity and amplitude of the oscillations of the bubble 
depend sensitively not only on the forcing pressure am- 
plitude, but also on the concentration of the gas dis- 
solved in the liquid, the temperature of the liquid, or 
small amounts of surface active impurities (Weninger 
et al., 1995; Ashokkumar et al., 2000; Toegel, Hilgen- 
feldt, and Lohse, 2000). As an example, Fig. 2 of Barber 
et al. (1994) shows the dependence of R(t) and the total 
light intensity on the increasing drive level for an air 
bubble in water. As the forcing is increased, the bubble 
size abruptly decreases, and then the light turns on (see 
Fig 7). For some years, the precise reasons for this sen- 
sitivity (observed repeatedly in experiments) were diffi- 
cult to understand, mostly because varying one of the 
experimental parameters, such as the water tempera- 
ture, would tend to change others as well. 

Perhaps most surprisingly, Hiller et al. (1994) found a 
sensitive dependence on the type of gas within the 
bubble: when the air dissolved in the liquid was replaced 
with pure nitrogen, the characteristically stable SBSL 
disappeared, With a gas composed of 80% nitrogen and 
20% oxygen, there was still no sonoluminescence. Only 
when the inert gas argon was added did SBSL light 
emission return. Figure 8 shows a plot of the intensity of 
sonoluminescence as a function of the percentage of in- 
ert gas doped in nitrogen. For both argon and xenon, the 
intensity peaks around 1%, the concentration of argon 
in air. 
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FIG. 7. The ambient bubble radius as a function of forcing 
pressure P, for a gas mixture of 5% argon and 95% nitrogen 
at a pressure overhead of 150 mm. For sonoluminescing 
bubbles the symbols are filled; for nonglowing bubbles they are 
open. Note the abrupt decrease in bubble size right before the 
sonoluminescence threshold. The figure is a sketch from Fig. 
38 of Barber et al. (1997). In that paper the ambient radius is 
obtained from a fit of the Rayleigh-Plesset equation to the 
R(t) curve. In that fit heat losses are not considered explicitly, 
but material constants are considered as free parameters. 
‘Therefore the values for Ry are only approximate; see the dis- 
cussion in Sec. ILE. 


SBSL can be achieved with a pure noble gas as well, 
but in a vastly different range of gas concentrations: In 
the original experiment with air, Gaitan (1990) observed 
stable light emission when degassing using a partial pres- 
sure of p2”/Py~0.2-0.4; i.e., the water contained 20— 
40 *6 of the air it would contain if in saturation equilib- 
rium with a pressure of Pọ=1 bar. Barber et al. (1995) 
demonstrated that, when using pure argon gas, the de- 
gassing has to be 100 times stronger, requiring partial 
pressures as low as p"/Pp~0.002—0.004 to obtain stable 
SBSL. The pressures p. are the partial gas pressures 
used in experiment when preparing the degassed liquid. 
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FIG. 8. Dependence of the sonoluminescence intensity (nor- 
malized to that of air) in water as a function of the percentage 
(mole fraction) of noble gas mixed with nitrogen. Two noble 
gases are shown: xenon (@) and argon (B). Both give maxi- 
mum light intensity around 1% dissolution, as does helium 
(not shown). The figure is a sketch from Fig. 22 of Barber et al. 
(1997). 
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During this time the predominant belief in the field 
was that shocks (see, e.g., Barber et al., 1994, 1997) were 
somehow important for the energy focusing and light 
emission of sonoluminescence. However, there was little 
agreement as to the details of how this worked, and 
many other physical mechanisms were suggested, includ- 
ing dielectric breakdown of the gas (Garcia and Le- 
vanyuk, 1996; Lepoint et al., 1997; Garcia and Hasmy, 
1998), fracture-induced light emission (Prosperetti, 
1997), bremsstrahlung (Moss, 1997; Frommhold, 1998), 
collision-induced emission (Frommhold and Atchley, 
1994; Frommhold, 1997; Frommhold and Meyer, 1997), 
and even the quantum-electrodynamical Casimir effect 
(Eberlein, 19962, 1996b), an idea pioneered in this con- 
text by Schwinger (1992). 

The difficulty in evaluating these ideas was that they 
required probing the bubble collapse in greater detail 
than was experimentally possible. This led Robert Apfel 
to pose a "challenge to theorists" in a session on sonolu- 
minescence at the annual meeting of the Acoustical So- 
ciety of America in Honolulu in 1996, The challenge was 
to make concrete, experimentally testable predictions. 
Many creative ideas were collected at this meeting, only 
a fraction of which still survive today. [One of the early 
casualties includes the acoustic-resonator theory devel- 
oped by the present authors speculating on energy stor- 
age in the bubble (Brenner, Hilgenfeldt, et al., 1996).] 

Meanwhile, it was equally clear that at least some of 
the experimental facts of sonoluminescence were direct 
consequences of the classical theory of bubble dynamics, 
having nothing to do with light emission per se. The time 
scale of the light emission is so much shorter than a 
complete cycle of the acoustic driving that bubble dy- 
namics goes a long way towards explaining issues of 
bubble stability and constraints for driving parameters. 
Since Lord Rayleigh's characterization of cavitation col- 
lapse (Rayleigh, 1917), bubble dynamics had become 
well understood,” but, although the theory was formally 
quite mature, it had never been put to work in the pre- 
cise regime of single-bubble sonoluminescence. 

The application of classical bubble dynamics to SBSL 
substantially clarified the experimental situation. The 
first contribution in this regard was made in the original 
paper of Gaitan ef al. (1992), which demonstrated that 
the radius of the bubble as a function of time observed 
experimentally exhibits the same behavior as solutions 
to the Rayleigh-Plesset equation (to be derived in Sec. 
1I); subsequently, studies by Löfstedt er al. (1993, 1995) 
confirmed and elaborated on this conclusion. The 
Rayleigh-Plesset theory is remarkably simple, and it 
captures many important features of single-bubble 
sonoluminescence. To practitioners of classical bubble 





"This was primarily due to the contributions of Plesset, 1949, 
1954; Epstein and Plesset, 1950; Plesset and Zwick, 1952; 
Plesset, 1954; Plesset and Mitchell, 1956; Eller and Flynn, 1964; 
Eller, 1969; Eller and Crum, 1970; Prosperetti, 1974, 1975, 
1977a, 197d; Plesset and Prosperetti, 1977; Prosperetti and 
Lezzi, 1986; Prosperetti et al., 1988, 
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FIG. 9, Phase diagram in the p/P) vs P,/P, parameter 
space, according to the hydrodynamic/chemical theory of 
Hilgenfeldt et al. (1996) and Lohse et al. (1997). The driving 
frequency is f=33.4 kHz. The three phases represent stable 
SL, unstable SL, and no SL. The symbols represent measure- 
‘ments by Ketterling and Apfel (1998), either stable sonolumi- 
nescing bubbles (#) or stable, nonsonoluminescing bubbles 
(O), showing good agreement with the earlier theoretical pre- 
dictions. 


dynamics, the excellent agreement was particularly sur- 
prising because this theory has long been known to show 
large quantitative discrepancies even for bubbles that 
are more weakly forced than in the case of SBSL (Pros- 
peretti er al., 1988). In the SBSL parameter regime, the 
periodic forcing of the pressure waves in the container 
leads to a periodic bubble response, with a cavitation 
collapse happening exactly once per cycle [chaotic mo- 
tion as in Lauterborn (1976) and Lauterborn and Suchla 
(1984) is notably absent]. The qualitative and even most 
quantitative features of bubble oscillations agree with 
the experimental observations. The solution also has the 
courtesy to predict its own demise: at cavitation collapse 
the speed of the bubble wall approaches or surpasses the 
speed of sound in the liquid, contradicting one of the 
essential assumptions of the theory. The total time dur- 
ing which the bubble wall is supersonic is a tiny fraction 
of a cycle; the errors that accumulate in this regime do 
not substantially affect the rest of the cycle. 

If the solutions to the Rayleigh-Plesset equation ex- 
plain the experimental measurements of the bubble ra- 
dius, then their stability must constrain the parameter 
space where SBSL can occur (Brenner et al., 1995; Bren- 
ner, Hilgenfeldt, et al., 1996, Hilgenfeldt et al., 1996). 
There are three major instabilities of the bubble that 
(i) the bubble must not change 
ies; Brenner et al., 1995; Hilgen- 
feldt et al., 1996); (ii) the average number of gas mol- 
ecules in the bubble must not increase or decrease over 
time (diffusive instability; Brenner, Lohse, et al., 1996; 
Hilgenfeldt er al., 1996); (iii) the bubble must not be 
ejected from the acoustic trap where it is contained 
(Bjerknes instability; Cordry, 1995; Akhatov er al., 1997; 
Matula et al., 1997). All of these constraints must be sat- 
isfied in a parameter regime where the bubble oscilla- 
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FIG. 10. Phase diagram for air at p,./Py=0.20 in the Ry-P, 
space. The arrows denote whether the ambient radius grows or 
shrinks at this parameter value. Curve A denotes the equilib- 
rium for an air bubble; on curve C the bubble contains only 
argon. The intermediate curve B necessarily exists because of 
the topology of the diagram and represents an additional 
stable equilibrium. The thin line indicates where the (approxi- 
mate) threshold temperature of nitrogen dissociation 
(79000 K) is reached. From Lohse et al. (1997). 





tions become nonlinear enough for sonoluminescence to 
occur. The allowable parameter space of SBSL is thus 
severely limited to a narrow range of relative gas con- 
centrations c../co=p./Pq and forcing pressure ampli- 
tudes P, (see Fig. 9). 

While the regime of stable sonoluminescence in argon 
gas is in good agreement with that predicted by the hy- 
hydrodynamic stability calculations of Hilgenfeldt er al. 
(1996), Barber et al. (1995) found that the ranges of dis- 
solved gas concentrations for stable SBSL were lower by 
a factor of 100 in pure argon gas than in air. Löfstedt 
et al. (1995) pointed out that a sonoluminescing bubble 
cannot possibly be in diffusive equilibrium for these pa- 
rameters and postulated another “anomalous mass 
flow,” whose mechanism would be “the key to SL in a 
single bubble.” 

To account for these discrepancies to classical bubble 
dynamics, Lohse et al. (1997) proposed that the extra 
mass-ejection mechanism of Léfstedt et al. (1995) is of a 
chemical nature. The gas in the bubble is hot enough 
upon collapse to allow for significant dissociation of Nz 
and O;. The dissociated nitrogen and oxygen, as well as 
some radicals from dissociated water vapor, will undergo 
chemical reactions, whose products are very soluble in 
water and are expelled from the bubble. Only inert, non- 
reactive gases (such as argon) remain inside according to 
this “argon rectification hypothesis." This idea immedi- 
ately resolves the apparent discrepancy between the 
measured and predicted parameter regimes for stable 
SBSL in air: if the bubble ends up filled with argon gas 
only, then only the argon dissolved in the liquid has to 
be in diffusive equilibrium with the bubble. As air con- 
tains 1% of argon, the effective dissolved gas concentra- 
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FIG. 11. Experimental phase diagram in the Ro-P, parameter 
space for air at px/Py=020. The driving frequency is 
20.6 kHz. Arrows indicate whether the bubbles grow or 
shrink. Three equilibrium curves A, B, and C can be recog- 
nized. In between curves B and C there is a "dissolution is- 
land,” The shaded area shows the shape-stable parameter do- 
main (sce Sec. IV.D). Figure adopted from Holt and Gaitan 
(1996). 


tion for diffusive stability of argon is 100 times smaller, 
and explains the hundredfold difference between ob- 
served concentrations for air and argon bubbles. The 
phase diagram in the Ro-P, space resulting from that 

















6 7 
dm m 
able 
4 í 
mo 
€ stable no SL 
á stable SL 
0. 
To 72 p Da 


P, (bar) 


FIG. 12. Experimental phase diagram for air saturated in wa- 
ter to 20%. Each data point represents the P, and Ro found 
from a single R(t) curve and is indicated to be luminescing 
and/or stable. The curves in the plot are lines of diffusive equi- 
librium for a given gas concentration c.,/¢y=0.2 (solid line) 
and c, /co70.002 (dashed line). The range of P, where danc- 
ing bubbles were observed is indicated, as are regions of 
bubble growth (g) and dissolution (d) relative to each equilib- 
rium curve. The stable no-SL points (8) correspond to a stable. 
chemical equilibrium which would lie above the c, /¢y=0.2 
curve if plotted. From Ketterling and Aptel (1998). 
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theory is shown in Fig. 10. In particular, the theory pre- 
dicts a new stable branch (called "B" in Fig. 10) on 
which mass losses from chemical reactions and growth 
from rectified diffusion just balance. 

Experiments by Holt and Gaitan (1996) on bubble 
stability published contemporaneously with the theoret- 
ical work indeed showed this extra regime of bubble 
stability predicted from the argon rectification hypoth- 
esis (see Fig. 11). Ketterling and Apfel (1998, 2000, 
2000b) later showed the stability predictions to be quan- 
titatively correct. Figure 12 shows experimental mea- 
surements of a phase diagram in comparison with theo- 
retical predictions. One consequence of the interplay of 
diffusive and shape instabilities is indicated in this figure: 
bubbles can “dance” due to the recoil when they un- 
dergo fragmentation (see Sec. IV.E). 

Phase diagrams such as Figs. 9-12 help us to under- 
stand the limitations of the parameter space for sonolu- 
minescence, and in particular the crucial role of noble 
gases for SBSL stability. The same theoretical concepts 
could be applied to explain the pronounced increase in 
the intensity of emitted light with decreasing water tem- 
perature (Hilgenfeldt, Lohse, et al, 1998), and the 
quenching of light due to small concentrations of surfac- 
tants, both of which were shown to be in agreement with 
experiments (Ashokkumar et al., 2000; Matula, 2000; 
Toegel, Hilgenfeldt, et al., 2000). 

‘There was, however, still the nagging problem of the 
light emission itself. In contrast to the bubble dynamics, 
the available experimental information was insufficient 
to constrain the theories. The breakthrough contribution 
was made by Gompf et al. (1997), who measured the 
width of the light pulse using time-correlated single- 
photon counting (TC-SPC). This technique has a much 
higher resolution for measuring flash widths than a 
single photomultiplier tube, because it measures time 
delays in arrivals of single photons. The measurement of 
the delay time between the two photons reaching the 
two different photomultiplier tubes is repeated many 
times so that the width of the flash can be reconstructed. 
Gompf et al. (1997) discovered that the width of the 
light pulse is actually of the order of a few hundred pi- 
coseconds (see Fig. 13), much longer than the previous 
50-ps upper bound measured by Barber and Putterman 
(1991). Moreover, since Gompf et al, (1997) could now 
resolve the shape of the light pulse, it was possible to 
study the dependence of the width on external param- 
eters (the forcing pressure and dissolved gas concentra- 
tion; see Fig. 14). 

After this paper was published, at a meeting on 
sonoluminescence at the University of Chicago, two 
other groups announced that they had confirmed its 
findings: Moran and Sweider (1998) and Hiller er al. 
(1998) also used TC-SPC. At the same time, Gompf's 
group succeeded in obtaining an independent confirma- 
tion of the much longer duration of the light pulse using 
a streak camera for direct measurement of the pulse 
width (Pecha et al., 1998). A previous attempt by Moran 
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FIG. 13. First measurement of SBSL pulse widths. The param- 
eters were P,=1.2 bars, f=20 kHz, and the gas concentration 
was 1.8-mg/l Oz. Both the width in the red and the ultraviolet. 
spectral range were measured. The indistinguishable widths 
rule out blackbody radiation, but not a thermal emission pro- 
cess in general. From Gompf et al. (1997). 





et al. (1995) employing a streak camera had yielded only 
a tentative upper bound for pulse width, which again 
proved too small. 

The increased experimental resolution of TC-SPC and 
the subsequent discovery of a long flash width put all of 
the theories of light emission and energy focusing, which 
required ultrashort flash widths, out of business. More- 
over, as was emphasized by Gompf et al. (1997) in their 
seminal paper, the measurement restored hope that a 
variant of the simplest possible theory for the light emis- 
sion might be correct the cavitation collapse of the 
bubble is so rapid that heat cannot escape from the 
bubble. Therefore, the bubble heats up, leading to light 
emission. Figure 15 shows the heating as calculated by 
Gompf eral (1997), by solving a variant of the 
Rayleigh-Plesset equation for the bubble radius and as- 
suming adiabatic heating (ratio of specific heats T 
/3) near the collapse. Although the calculation con- 
tains some severe approximations, the agreement is 
quite reasonable. 

This idea was buttressed by an earlier numerical simu- 
lation of Vuong and Szeri (1996), which, when reinter- 
preted with the new experiments in mind, questioned 
the notion that strong shocks are important for single- 
bubble sonoluminescence. Vuong and Szeri included dis- 
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FIG. 14. Dependence of the full width at half maximum of the. 
SBSL pulse on the driving pressure and the gas concentration 
at room temperature. f=20 kHz. From Gompf et al. (1997). 
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FIG. 15. Calculated shape of temperature pulse using a simple 
model based on the Rayleigh-Plesset equation, assuming the 
gas temperature and density are uniform throughout the col- 
lapse. Figure reproduced from Gompt et al. (1997). 


sipative effects and showed that the strong shocks pre- 
dicted by Wu and Roberts (1993) and Moss et al. (1994, 
1996, 1997) were absent in noble gas bubbles, and were 
replaced by gentler inhomogeneities. The predicted 
maximum temperatures in the bubble were therefore 
much lower, several 10* K, compared with the 10° K 
previously announced by Wu and Roberts (1994). More- 
over, the hot spot was not highly localized in the bubble 
center These arguments were elaborated upon by 
‘Vuong et al. (1999); these models are much closer to the 
simple picture of adiabatic heating and thermal light 
emission than the shock-wave scenario. The tempera- 
ture profiles and motions of Lagrangian points as com- 
puted by Vuong and Szeri (1996) are shown in Fig. 16: 
‘The characteristic scale over which temperature varies is 
of the order of the bubble radius. 

Since the experimental resolution of the flash, re- 
searchers have focused on trying to determine which 
variant of the thermal light-emission model is correct. Is 
the interior of the bubble uniform? Is the radiation 
blackbody, bremsstrahlung, or some other process? Is 
the bubble optically thin or thick? What physical mecha- 
nism is suppressing spectral lines? Since experiments are 
now able to measure both the shape of the light pulse 
and the spectrum independently and accurately, it is pos- 
sible to determine how these quantities depend on ex- 
perimental parameters like forcing pressure, gas concen- 
tration, etc. The power of these measurements is that 
they provide severe constraints for theories of SBSL 
light emission that did not exist when the pulse width 
was believed to be very short. Moreover, since the 
bubble dynamics itself is well understood, closer exami- 
nation of these parameter dependencies makes it pos- 
sible to focus attention on subtle details of the light- 
emitting process. Single-bubble sonoluminescence has 
thus become a rather sophisticated testing ground for 
the ability of mathematical models and numerical simu- 
lations to explain detailed experimental data from a 
complicated physical process. 

Although there are still open questions about the de- 
tails of the light emission, considerable progress has 
been made. When Gompf et al, (1997) resolved the light 
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FIG. 16. Motion and temperature in a bubble shortly before 
collapse: (a) motion history of 20 Lagrangian points inside a 
Ro=45 um bubble driven at P,-13 atm and f=26.5 kHz. 
Strong wavy motion occurs inside the bubble, but no shock 
waves develop. (b) Temperature profiles in the bubble for vari- 
ous times around the bubble collapse. The profiles span a time 
interval of ~170 ps near the collapse. The temperature at the 
center increases monotonically, until the maximum tempera- 
ture is reached at the last snapshot. Note that the temperature. 
profile is smooth, without any discontinuity that would be 
present with a shock. From Vuong and Szeri (1996). 





pulses, they also made measurements of the dependence 
of the width on optical wavelength. Strikingly, such a 
dependence was found to be absent, contradicting a 
simple blackbody emission model, which demands that 
the width increase with the wavelength. 
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FIG. 17. Emission spectra from rare gases at room tempera- 
ture. The dotted lines are calculations based on the theoretical 
model of Hammer and Frommhold (2000). The only adjust- 
able parameters in the comparison are the ambient radii and 
forcing pressures of the bubbles. From Hammer and Fromm- 
hold (2001). 


A resolution for this conundrum was hinted at in nu- 
merical simulations by Moss et al. (1994, 1997, 1999), 
who realized that the temperature-dependent photon 
absorption coefficients of the gas must be taken into ac- 
count. The size of the bubble and thus the size of the 
light-emitting region are so small that the bubble is 
nearly transparent for its own photons: the bubble is a 
volume emitter, not a surface emitter like an ideal black- 
body. Among other things, Moss et al. (1999) used this 
idea to rationalize the qualitative shape of the emission 
spectrum in noble gases. 

Hilgenfeldt et al. (19992, 1999b) used varying absorp- 
tion coefficients to explain the wavelength-independent 
pulse widths: Both the absorptivity and emissivity of the 
bubble drop precipitously directly after collapse for all 
wavelengths, since they depend exponentially on tem- 
perature, but only weakly on wavelength. Combining 
this model of thermal radiation with the parameter de- 
pendencies predicted by the stability constraints on the 
bubble, they also found agreement with the observed 
parameter dependencies of the pulse width, number of 
photons per burst, and spectral shape. Hammer and 
Frommhold (2000a, 2000b) demonstrated that this 
model could be refined with ab initio quantum- 
mechanical calculations of electron-neutral bremsstrah- 
lung, further improving the agreement with experi- 
ments. Examples of their spectra are shown in Fig. 17. 

‘An important aspect of bubble thermodynamics, 
which has been pointed out by Kamath et al. (1993), Ya- 
sui (19975), Colussi and Hoffmann (1999), Moss et al. 
(1999), Storey and Szeri (2000, 2001); Toegel, Gompt, 
et al. (2000), Hilgenfeldt er al. (2001), and Putterman 
et al. (2001), is the presence of water vapor inside the 
bubble. Upon bubble expansion, vapor invades the 
bubble. At collapse, it cannot completely escape (con- 
dense at the bubble wall) because the diffusion time 
scale is much slower than the time scale of the collapse. 
Therefore water vapor is trapped inside the bubble (Sto- 
rey and Szeri, 2000). It limits the maximum temperature 
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FIG. 18. Dependence of the spectra of argon SBSL (for a 
partial pressure of ~150 torr at 25 °C) on the forcing pressure. 
Spectra are shown for five levels of overall brightness. The OH 
line is vanishing in the thermal bremsstrahlung spectrum with 
increasing forcing pressure P, . From Young et al. (2001). 


in the bubble due to its lower polytropic exponent (com- 
pared to inert gases) and above all because of the endo- 
thermic chemical reaction H,O-+OH+H, which eats up 
the focused energy. Within the model of Storey and 
Szeri (2000), taking water vapor and its chemical reac- 
tions into account leads to calculated maximum tem- 
peratures in the bubble of only around 6000 K. This 
seems to contradict experiments, in that thermal light 
emission would be strongly suppressed below the mea- 
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FIG. 19. Light-emission spectra from moving SBSL bubbles in 
adiponitrile. The driving pressure amplitude increases from 
bottom to top, between 1.7 bars and 1.9 bars. The spectral line 
at ~400 nm corresponds to an excitation of CN. From 
Didenko et al. (2000b). 
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sured values, and is an indication that the modeling 
overestimates the amount of water vapor in the bubble. 

In two very recent experiments, the signatures (char- 
acteristic lines) of the liquid or liquid vapor were de- 
tected in the spectrum, finally closing the gap between 
MBSL and SBSL. In both cases the lines belong to con- 
stituents of vapor molecules. 

Young et al. (2001) discovered spectral lines for SBSL 
in water by decreasing the driving pressure very close to 
the threshold for SBSL. In this regime, the light pulse is 
so weak that Young ef al. (2001) had to collect photons 
over several days. Figure 18 shows how, as the forcing 
pressure is increased, the OH line vanishes behind the 
enhanced continuum contribution to the spectrum. 

Didenko et al. (2000b) found spectral lines of SBSL in 
organic fluids (see Fig. 19). These tend to require larger 
driving to show SBSL, because the vapor molecules 
have more rotational and vibrational degrees of free- 
dom, leading to a weaker temperature increase at 
bubble collapse. 

We believe that the observation of spectral lines her- 
alds a new era of research on single-bubble sonolumi- 
nescence, one in which it will be possible to use SBSL to 
study chemical reactions. Such studies have long been 
conducted for multibubble cavitation, and indeed Sus- 
lick and collaborators (Suslick ef al., 1986; Flint and Sus- 
lick, 1991b; Didenko er aL, 1999) have used the widths 
and intensities of spectral lines in multibubble sonolumi- 
nescence to deduce the temperature of cavitation. The 
great advantage of using single-bubble sonolumines- 
cence in these studies is that, in contrast to MBSL, the 
mechanics of SBSL is well understood and character- 
ized. It thus seems possible that one will be able to use 
SBSL to carefully study chemical reactions under exotic 
conditions of high temperatures and extreme densities. 


Il. FLUID DYNAMICS OF THE FLASK 


The very existence of a sonoluminescing bubble de- 
pends critically on a subtle balance of hydrodynamic and 
acoustic forces inside the flask. During sonolumines- 
cence, a diverse array of physical effects influences this 
balance: the pressure becomes low enough that the 
liquid-air interface vaporizes, and temperatures rise so 
high that the gas inside the bubble emits light. Gas is 
continually exchanged between the bubble and the sur- 
rounding liquid, causing the number of molecules in the 
bubble to vary. In a small part of the cycle, the bubble- 
wall velocity may become supersonic, During all of these 
processes there is no a priori reason for the shape of the 
bubble to remain spherical, so this must be accounted 
for as well. 

Although the equations of motion governing these ef- 
fects were written in the nineteenth century, it is a tri- 
umph of twentieth-century applied mathematics that all 
of them can be accounted for simultaneously in a precise 
and controlled way. This is the theory of classical bubble 
dynamics, started by Lord Rayleigh (1917) during his 
work for the Royal Navy investigating cavitation dam- 
age of ship propellers. The formalism was substantially 
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refined and developed by Plesset, Prosperetti, and oth- 
ers over a span of several decades. A review of early 
work is presented by Plesset and Prosperetti (1997); a 
later overview is given by Prosperetti (1998). The 
present section summarizes this theory with a view to- 
wards its application to experiments on single-bubble 
sonoluminescence. Our discussion will highlight the va- 
lidity of the approximations made when the theory is 
applied to SBSL, and will also underscore how and why 
the theory works when it does. The presentation of this 
section was greatly influenced by the excellent recent 
review by Prosperetti (1998). 


A. Derivation of the Rayleigh-Plesset equation 


The ultrasonic forces in the liquid are caused by the 
oscillating transducers on the container walls, which are 
tuned to excite an acoustic resonance mode of the con- 
tainer, often the lowest. The Q factor of a typical flask is 
710", so the resonance is quite sharp. Its frequency is 
about 20 kHz for a container a few centimeters across, 
mercifully above the range of human hearing. The driv- 
ing pressure amplitude at the center of the flask is 
around P,~1.2~1.4 bars when SBSL occurs. 

The equations governing the sound waves in the liq- 
uid are the compressible Navier-Stokes equations 


Vp nVut (VV-u, (en) 
. Q 


where u is the fluid velocity, p the density, p the pressure 
(as specified by an equation of state), 7 the shear viscos- 
ity, and 7 the bulk viscosity of the liquid. In writing these 
equations, we have assumed that the liquid is isothermal 
and so have neglected the equation for the fluid tem- 
perature. As an approximation, the bubble's extension 
compared to that of the flask and that of the sound wave 
is neglected, as it is orders of magnitude smaller. 

The forces on the bubble depend on where it is 10- 
cated in the flask. In general there will be both an iso- 
tropic oscillatory pressure (causing volumetric oscilla- 
tions) and, in addition, pressure gradients, quadrupole 
components, etc. In practice, for small bubbles, all that 
matters are the isotropic volumetric oscillations and the 
pressure gradients, which can create a net translational 
force on the bubble. The translation can vanish only at 
pressure maxima or minima, We shall see below that 
these forces cause sonoluminescing bubbles to be 
trapped at a pressure antinode of the sound field. 

To compute the magnitude of the forces it is necessary 
first to characterize the volumetric oscillations, for which 
the sound field around the bubble is purely radial. The 
velocity can then be represented by a potential, with u 
=V¢. Equations (1) and (2) then become 
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Efforts to scale up sonoluminescence have ventured into the 
lower-frequency regime of audible sound. Bad luck for the ex- 
perimentalist. 
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Note the assumption that the flow field is purely radial 
and therefore viscous stresses are not important, 

To proceed we need to combine Eqs. (3) and (4) into 
a single equation for 4. Defining the enthalpy dH 
—dplp, and using dp=(dpldp)dp=cdp (with c the 
speed of sound in the liquid) implies 
u 1 
Vidc [021 —2,H) + rh, 6 
where u— 2,4 is the radial velocity field. As long as the 
fuid velocity is much smaller than c, the square- 
bracketed terms are negligible. The linear c 22d term 
is only negligible close to the bubble: at distances on the 
order of the sound wavelength away from the bubble, 
this term will become important. 

We would like to solve Eq. (5) for the velocity of the 
bubble wall dR/dt, caused by the resonant oscillation of 
the container. We proceed in two steps: near the bubble 
the velocity potential obeys the Laplace equation, V?76 
—0. The solution satisfying the boundary condition at 


the bubble wall 4,6(r=R)=R is 








RR 
g=-——+A(0, (6) 
where A(t) is a free constant, This free constant is de- 
termined by matching the solution (6) onto the pressure 
field far from the bubble. Neglecting the sound radiated 
by the bubble itself, the velocity potential far from the 
bubble is a standing wave—the acoustic mode that is 
excited by the transducer. For our present purposes, we 
do not require the entire spatial structure of this mode, 
but only the field close to the bubble. Since the bubble is 
much smaller than the sound wavelength, this sound 
field will be independent of r, so that $= @.(1). Match- 
ing the near field and the far field implies A=... The 
pressure in the neighborhood of the bubble is then p= 
~pajde=Po+P(t), ie, the sum of the background 
static pressure P=1 bar and the sinusoidal driving pres- 
sure P(t) - — P, sin wt. 

The velocity field in the liquid around the bubble then. 
follows as 








RR? 
w=. 


0o 


We now use this to solve for the dynamics of the bubble 
wall. To this end, we use the force balance on the bubble 
surface, which gives 


n4) Xubr- RG)- pf) PER] 2 92,07 R) 





(8) 
where X,, is the radial component of the stress tensor in 
the liquid, ø is the surface tension of the gas-liquid in- 
terface, and p, is the pressure in the gas, assumed to be 


Roc 
—p(0-pIRG]-42 8 72. 
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spatially uniform. Evaluating this formula using Eq. (3) 
for the pressure in the liquid gives 


(9) 





Equation (9) is the celebrated Rayleigh-Plesset equa- 
tion. The left-hand side of the equation was known to 
Lord Rayleigh (though never written). A historical re- 
view of the development of this equation is given by 
Plesset and Prosperetti (1977). 

Closing the equation requires knowing the pressure in 
the gas. Roughly speaking, when the bubble wall moves 
slowly with respect to the sound velocity in the gas, the 
pressure in the gas is uniform throughout the bubble. In 
this regime, how strongly the pressure depends on the 
bubble volume depends on the heat transfer across the 
bubble wall (Prosperetti etal., 1988). The pressure- 
volume relation is given by 
(R$- mr 


2 
puo Pot RRO ao) 





Here Ry is the ambient radius of the bubble (ie., the 
radius at which an unforced bubble would be in equilib- 
rium), and A is the van der Waals hard-core radius de- 
termined by the excluded volume of the gas molecules. 

If the heat transfer is fast (relative to the time scale of 
the bubble motion), then the gas in the bubble is main- 
tained at the temperature of the liquid, and the pressure 
is determined by an isothermal equation of state with 
y=1. On the other hand, if the bubble wall moves very 
quickly relative to the time scale of heat transfer, then 
heat will not be able to escape from the bubble, and the 
bubble will heat (cool) adiabatically on collapse (expan- 
sion). For a monatomic (noble) gas, this implies that y 
5/3. The dimensionless parameter that distin- 
guishes between these two regimes is the Péclet number, 








, a1) 


where x, is the thermal diffusivity of the gas. 

This idea about heat transfer is based on a more care- 
ful version of this argument by Kamath et al. (1993) and 
Prosperetti et al. (1998). They showed that the tempera- 
ture 7, at the bubble surface is basically the water tem- 
perature: Conservation of energy at the bubble interface 
requires continuity of the heat flux, 


Kyu T-Kj)Ty, (12) 
with the thermal conductivities K, and K; of gas and 
liquid. The gradients are estimated via the thermal 


boundary layer thicknesses 5, and & in and around the 
bubble, 








T= T, 
iT a EE 
where T, is the temperature at the bubble center. The 
diffusion lengths can be estimated with the relevant time 
scale At of the bubble oscillation and the respective 
thermal diffusivity y, namely, 6~ /xAf. With the con- 


(13) 
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nection between thermal conductivity and diffusivity, K 
=xpCp , where C, is the specific heat per unit mass, one 
obtains the final result, 


XzPzC, 
xe, 





(14) 


T, , 
Since the density and the specific heat of water are so 
much larger than the respective values for gas, the right- 
hand side of Eq. (14) is typically of the order of 


1073-1072, Therefore the temperature drop basically 
occurs inside the bubble, and the temperature at the sur- 
face basically equals the water temperature, 

If the rate of heat transfer is intermediate between 
adiabatic and isothermal, the situation is more compli- 
cated. Here, a correct calculation requires solving the 
heat conduction problem throughout the bubble cycle 
and using the computed temperature in the bubble to 
evaluate the pressure in the pas (through its equation of 
state). This is quite a difficult task. Over the years, sev- 
eral methods have been proposed that amount to vary- 
ing y continuously between the isothermal value and the 
adiabatic value (Plesset and Prosperetti, 1977; Prosper- 
etti et al., 1988; Kamath et al., 1993) depending on the 
Péclet number. This approach can yield quantitatively 
incorrect results, as shown by Prosperetti and Hao 
(1999), in large part because energy dissipation from 
thermal processes is neglected. 





B. Extensions of the Rayleigh-Plesset equation 


So far we have not considered damping of the bubble 
dynamics by the sound radiated by the bubble itself. The 
most complete and elegant derivation of this effect is 
due to Lezzi and Prosperetti (1987; Prosperetti and 
Lezzi, 1986). 

In arriving at Eq. (9), we asserted that the velocity 
potential of the sound field in the liquid far from the 
bubble is the same as in the absence of the bubble, $ 
=¢.(t). The radial sound wave emitted from the 
bubble introduces a modification, 


1 1 F(t) 
$= $.(0)— > Flt—rle)~balt)— =F) +—, (15) 


where we have estimated the velocity potential at small 
r. As above, this now must be matched to the near-field 
velocity potential Eq. (6). The matching yields F(t) 








=R?R and A(f)=¢.+Fie. Substituting this into the 
pressure jump condition one obtains 
ig um 
p( RR 4R )-D. 7 Po- P(0]- 40g 720 
pd R 
+E Fal (16) 





The sound radiation term is of order R/c times the other 
terms in the equation. When the bubble-wall motion is 
slow it is therefore negligible. 
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When |R|/c—1, sound radiation is important. For- 
mally, sound radiation raises the order of the Rayleigh- 
Plesset equation from second order to third order. At 
first glance, this seems strange, because physically initial 
conditions are given for both R and R, but not Ë. The 
discrepancy arises because Eq. (16) has a spurious un- 
stable solution which grows exponentially in time. This 
is unphysical; the initial condition on & must be chosen 
to suppress this solution. 

As emphasized by Prosperetti et al. (1988; Prosperetti 
and Hao, 1999), this procedure is inherently impractical, 
as numerical errors will always excite the spurious solu- 
tion. A better way to take care of this is to calculate the 
d'|di? (R?R) term using the Rayleigh-Plesset equation 
itself. A standard way of doing this was invented by 
Keller and co-workers (Keller and Kolodner, 1956; 
Keller and Miksis, 1980) and leads to the Keller equa- 
tion (Prosperetti and Lezzi, 1986; Brennen, 1995) 





-| 1+ Alp Po ron 


R 2e 


Ang (7) 





As discussed by Prosperetti et al. (1988; Prosperetti and. 
Lezzi, 1986), the precise form of this equation is not 
unique: There is a one-parameter family of equations 
that can be consistently derived from Eq. (16), namely, 





R 3 R 
[1-002 ores Rp 1-04 53 


R 
-| 10-07 |ip,-2a- r0 


R R 20 


tob ARR (18) 


where the parameter value A=0 recovers the Keller 
equation, and A=1 results in the formula used by Her- 
ring (1941) and Trilling (1952). Introducing higher-order 
terms leads to variations like the form derived by Flynn 
(1975a, 1975b), but Prosperetti and Lezzi (1986) have 
shown that the higher order does not, in general, guar- 
antee higher accuracy of the formula. Other well-known 
forms of Rayleigh-Plesset derivatives are compared by 
Lastman and Wentzell (1981, 1982). Prosperetti and 
Lezzi (1986) demonstrate that, for a number of relevant 
examples, the Keller equation yields results in closest 
agreement with full partial differential equation numeri- 
cal simulations. 

An “odd cousin” of Eq. (18) is the Gilmore equation 
(Gilmore, 1952; Brennen, 1995), 
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(19) 


whose derivation relies on the Kirkwood-Bethe approxi- 
mation (Kirkwood and Bethe, 1942). In Gilmore's equa- 
tion, the key quantity is the enthalpy H, and not the 
pressure. In this approach, the speed of sound C is not a 
constant, but depends on H. According to Gompf and 
Pecha (2000; Pecha and Gompf, 2000), this allows one to 
model the increase of the speed of sound with increasing 
pressure around the bubble, which leads to significantly 
reduced Mach numbers at bubble collapse. 

The breakdown of the Rayleigh-Plesset variants when 
|R|/c approaches unity is reflected in unphysical singu- 
larities when |R|/c—1 in the major terms of the equa- 
tions. Since equations with different à lead to similar 
results, one solution to this problem is to delete all the 
prefactors in parentheses containing R/c. We thus arrive 
at a popular form in the context of sonoluminescence 
(see, for example, Löfstedt er al., 1995; Barber et aL, 
1997), 


R 1 
o(RR+ 3R?)=[p,— Po PQ]-49 g -20 


Rd 2 
toat) (20) 

For very strong forcing, these different equations de- 
viate in the small time interval of bubble collapse, 
though they are in near-perfect accord for the rest of the 
driving cycle. Therefore they can be expected to pro- 
duce quantitative discrepancies for the properties of the 
collapsed bubble (e.g., the minimum radius, maximum 
gas pressure, etc.). These discrepancies are a principal 
source of modeling error for theories of SBSL. Another 
is the treatment of heat exchange via an effective poly- 
tropic exponent in Eq. (10). Simple refinements for heat 
exchange have been employed by Yasui (1995), though 
the only infallible solution is a direct calculation of the 
heat transfer. This was first carried out in numerical 
simulations by Vuong and Szeri (1996) and more re- 
cently by Moss et al. (1999). 

Given these difficulties, it is surprising that solutions 
to Rayleigh-Plesset-type equations still provide a quan- 
titatively accurate representation of the mechanics of a 
sonoluminescing bubble and of many of its accompany- 
ing effects. Recently, Lin et al. (2001) achieved a better 
understanding of why finite Mach number corrections to 
Rayleigh-Plesset-type equations are relatively unimpor- 
tant. They showed that the Rayleigh-Plesset equation is 

ite accurate even with significant spatial inhomogene- 
ities in the pressure field inside the bubble. This extends 
the utility of the Rayleigh-Plesset equation into the re- 
gime where the Mach number for the gas M,=Ric, 
(where c, is the speed of sound in the gas) is no longer 
small. Lin et al. (2001) show that the relevant condition 
is not [M,|<1, but |ej| c1, where 
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FIG. 20. Solutions to the modified Rayleigh-Plesset Eq. (20) at 
forcing pressures P,=1.0, 1.1, 12, and 1.3 atm. The ambient 
bubble radius is Ro~2 um, the frequency f=1/Ty=26.5 kHz. 
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.e., what is relevant is the bubble-wall acceleration. So 
even in the sonoluminescence regime, Lin et al. (2001) 
find excellent agreement when comparing their full gas- 
dynamical partial differential equation simulations with 
the solutions to the Rayleigh-Plesset ordinary differen- 
tial equation with the assumption of a uniform pressure 
inside. They also developed an approximation for the 
internal pressure field, taking into consideration first- 
order corrections from pressure inhomogeneity. 

In the remainder of this section, we present calcula- 
tions and experiments on bubble dynamics during a 
cycle of the driving, discussing the various physical ef- 
fects that are important away from the bubble collapse. 
Later sections will describe our present knowledge of 
the collapse itself. 


C. The bubble's response to weak and strong driving 


First, to give some feeling for solutions to the 
Rayleigh-Plesset equation, we study small oscillations of 
the bubble about its ambient radius Rọ. A straightfor- 
ward calculation (Brennen, 1995) shows that such a 
bubble oscillates at the resonant frequency 


2mfy— 


A typical sonoluminescing bubble has Ry~5 um, corre- 
sponding to a resonant frequency of f~0.5 MHz, much 
higher than the frequency of the driving f~20 kHz. 
Figure 20 shows solutions to the modified Rayleigh- 
Plesset Eq. (20) for a bubble at different forcing pres- 
sures, At low forcing, the bubble undergoes almost sinu- 
soidal oscillations of relatively small amplitude, with a 
period equal to that of the external forcing f. Here, the 
oscillations are essentially “quasistatic,” because’ the 
resonant frequency is so much larger than f: the oscilla- 


Q2) 
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tory pressure forcing is balanced by the gas pressure 
(Löfstedt et al., 1993; Hilgenfeldt, Brenner, et al, 1998), 
with inertia, surface tension, and viscosity playing a neg- 
ligible role. At a critical pressure around P,~Po, such 
quasistatic oscillations are no longer possible, resulting 
in a nonlinear response of the bubble. The critical P, 
depends slightly on Ro, and is referred to as the (dy- 
namical) Blake threshold (Blake, 1949; see also Hilgen- 
feldt, Brenner, etal, 1998). Beyond this threshold, 
sonoluminescence can occur. 

In the SBSL regime, the solution to Eq. (20) in this 
regime can be divided into several different stages. 


* Expansion: During the negative half-cycle of the driv- 
ing, the applied tension makes the bubble expand. 
Since f<f, the expansion continues until the applied 
pressure becomes positive. The time scale of this re- 
gime is thus set by the period of the driving pressure 
wave and is typically ~20 us for sonoluminescence 
experiments. This is sufficient to increase the bubble 
radius by as much as a factor of 10. 

* Collapse: When the driving changes sign, the ex- 
panded bubble is "released" and collapses inertially 
over a very short time (~1 ns for SBSL bubbles). The 
solution during collapse is well described by the clas- 
sical solution of Lord Rayleigh. SBSL light emission 
occurs at the end of the collapse. 

* Afterbounces: After the collapse, the bubble spends 
the remaining half of the cycle oscillating about its 
ambient radius at roughly its resonant frequency fo 
>f, giving rise to characteristic “afterbounces.” 


It is worthwhile at this point to comment on the roles 
of surface tension and viscosity. The surface tension 
term is dynamically important when it is as large as the 
external forcing pressure, implying that o/R~P,. This 
occurs when the bubble radius is smaller than Ry 
=o/P,. For water, this corresponds to a radius of 
~0.7 umi (P, /bar). We shall see below that this length 
scale plays an important role in determining the stability 
of the solutions to the Rayleigh-Plesset equation with 
respect to both dissolution and breakup. 

Viscous effects are important when the viscous damp- 
ing time scale is of the order of the time scale of bubble 
motion, roughly v/Rĝ~fo, with the kinematic viscosity 
v= n/p. For water, this does not occur; for more viscous 
fluids it can be important (Hilgenfeldt, Brenner, et al., 
1998). 


D. The Rayleigh collapse 


Now we turn to the behavior of the bubble radius 
near the collapse. As emphasized above, this is the re- 
gime in which the Rayleigh-Plesset description is in dan- 
ger of breaking down. The approach to the collapsed 
state, however, can be captured very well by the equa- 
tion, and is given by a classical solution of Lord Ray- 
leigh. 

Lord Rayleigh (1917) imagined a bubble dynamics for 
which only liquid inertia mattered, with gas pressure, 
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FIG. 21. Measured R(t) (with Mie scattering, dots) and a fit to 
these data based on the Keller equation (solid curve). The thin 
curve shows the driving pressure P(t). From Matula (1999). 


surface tension, and viscosity all negligible—in other 
words, the collapse of a void. The equation for the wall 
motion of the bubble/void is then RR+3/2R2=0 and 
can be directly integrated. The solution is of the form 
R()=Rol(t,.—/t, P^, with the remarkable feature of 
a divergent bubble-wall velocity as £ approaches the 
time £, of total collapse. Lord Rayleigh pointed out that 
this singularity is responsible for cavitation damage, and 
it is also the central hydrodynamic feature responsible 
for the rapid and strong energy focusing that leads to 
sonoluminescence. 

Clearly, something must stop the velocity from diverg- 
ing. For the Rayleigh-Plesset Eq. (9) to capture sonolu- 
minescence, it must contain the physical effect that does 
this. Viscous stresses 4 7R/Re-(t,—1)~} and surface ten- 
sion forces o/ R(t, — 1)-?5 diverge at slower rates than 
the inertial terms [2(t, —1)-55] and are therefore too 
weak. What about the gas pressure? The collapse rate is 
eventually so fast that the heat does not have time to 
escape the bubble. The pressure in the gas then obeys 
the adiabatic equation of state, which diverges as p, 
R^ «(t,—1) ? (for a monatomic ideal gas with T 
=5/3), which is stronger than the inertial acceleration. 
This effect is therefore capable of stopping the collapse, 
Modifications from the ideal gas law, e.g., van der Waals 
forces [see Eq. (10)], do not affect this conclusion. 

Although the gas pressure can halt Rayleigh collapse, 
it turns out that the most strongly divergent term in Eq. 
(20) is the last one, associated with sound radiation into 
the liquid during the last stages of collapse; it diverges as 
(t,—1)-P5 (Hilgenfeldt, Brenner, etal., 1998), and 
overwhelms the other terms. Up to 50% of the kinetic 
energy in the collapse may end up as a radiated pressure 
wave (Gompf and Pecha, 2000). 


E. Comparison to experiments 


Of course, it is crucial to compare solutions of 
Rayleigh-Plesset equations to experimental data on the 
bubble radius as a function of time. However, neither 
the ambient bubble radius Rp nor the driving pressure 
P, is known a priori. Ro changes through gas diffusion 
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as well as evaporation/condensation of water vapor (see 
Sec. III), and the (local) driving pressure P, is very sen- 
sitive to perturbations of the flask geometry, such as 
might be caused by a small hydrophone attempting to 
measure P,,. In addition, the precision of such a hydro- 
phone is limited to roughly 0.05 bar. 

The standard procedure has been to measure R(t) 
with Mie scattering? and then to fit the data to Rayleigh- 
Plesset-type dynamics by adjusting Ry and P,. A typical 
trace for a sonoluminescing bubbles radius during a 
cycle of the drive is shown in Fig. 21. The filled circles 
represent experimental measurements, and the solid line 
is a solution to the Keller equation under the assump- 
tion of isothermal heating (y=1). Superimposed as a 
thin line is the applied forcing pressure, 

The problem with these fits is that Ry and P, sensi- 
tively depend on model details. In particular, if one ad- 
justs Ro and P, such that the bubble's maximum is well 
fitted, the afterbounces are always overestimated (see 
Fig. 21). Better fits can be achieved by allowing more 
parameters, e.g., by allowing the material constants such 
as the viscosity or the surface tension to vary, Barber 
et al. (1992), for example, used seven times the usual 
value of the viscosity of water to achieve a fit to the 
afterbounces. As clarified by Prosperetti and Hao 
(1999), the larger viscosity effectively parametrizes other 
damping mechanisms not captured in simple Rayleigh- 
Plesset-type models. In particular, Prosperetti and Hao 
(1999) included thermal losses, following Prosperetti 
(1991), reducing the size of the afterbounces. Yasui 
(1995) had some success by introducing thermal bound- 
ary layers as well. 

Another effect that must be considered when fitting 
experimental R(t) curves to Rayleigh-Plesset models is 
the invasion of water vapor at bubble maximum. This 
leads to a varying ambient radius Ro over the bubble 
cycle, being largest at maximum radius, Since many 
early fits of R(t) curves (summarized by Barber et al., 
1997) did not consider these effects, the resulting values 
for Ro and P, are only approximate. 

Mie scattering data near the collapse are also notori- 
ously difficult to interpret because of the unknown index 
of refraction inside the compressed bubble and because 
the bubble radius R becomes of the order of the light 
wavelength. The simple proportionality of Mie intensity 
and R?, valid for larger R, gets lost and the relation 
even becomes nonmonotonic (Gompf and Pecha, 2000). 
Moreover, at collapse, the light is reflected not only from 
the bubble wall, but also from the shock wave emitted 
from the bubble at collapse. This subject will be treated 
in the next subsection. 

Another light-scattering technique based on differen- 
tial measurement and polarization (differential light 
scattering) has been developed by Vacca et al. (1999) in 





"Sce, for instance, the work of Gaitan, 1990; Barber et al., 
1992, 1997; Gaitan et al., 1992; Lentz et al, 1995; Weninger, 
Barber, and Putterman, 1997; Matula, 1999; Gompf and Pecha, 
2000; Pecha and Gompf, 2000; Weninger et aL, 2000. 
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FIG. 22. Outgoing shock wave from a collapsing bubble: (a) 
Streak image of the emitted outgoing shock wave from the 
collapsing bubble and (b) an intensity cross section along the 
line AA’. From Pecha and Gompf (2000). 


order to measure the dynamics of the bubble radius 
With this technique a time resolution of up to 0.5 ns 
around the Rayleigh collapse has been achieved. 


F. Sound emission from the bubble 


The Rayleigh-Plesset equation predicts the response 
not only of the bubble radius, but also of the surround- 
ing liquid. This has been detected by Cordry (1995), 
Holzfuss, Riiggeberg, and Billo (1998), Matula et al. 
(1998), Wang er al. (1999), Gompf and Pecha (2000), 
Pecha and Gompf (2000), and Weninger et al. (2000). 
Matula et al. (1998) used a piezoelectric hydrophone to 
measure a pressure pulse with fast rise time (5.2 ns) and 
high amplitude (1.7 bars) at a transducer at 1-mm dis- 
tance from the bubble. Wang et al. (1999) carried out a 
systematic study of the strength and duration of the 
pressure pulses as a function of gas concentration, driv- 
ing pressure, and liquid temperature. They demon- 
strated that a probe 2.5 mm from the bubble observes 
pressure pulses with rise times varying from 5 to 30 ns as 
the driving pressure and dissolved gas concentration 
vary. The amplitude of the pressure pulses varies be- 
tween 1 and 3 bars. 

Another study of this type was carried out by Pecha 
and Gompf (2000; Gompf and Pecha, 2000). They mea- 
sured pressure amplitudes and rise times consistent with 
the other measurements, and were able to measure the 
pressure pulse much closer (within 50 um) to the 
bubble. In addition, using a streak camera and shadow- 
graph technique, they visualized the shock wave leaving 
the bubble (see Fig. 22). Pecha and Gompf (2000) found 
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that the shock velocity in the immediate vicinity of the 
bubble is as fast as 4000 m/s, much faster than the speed 
of sound c= 1430 m/s in water under normal conditions, 
but in good agreement with the results of Holzfuss, Riig- 
geberg, and Billo (1998). This high shock speed origi- 
nates from the strong compression of the fluid around 
the bubble at collapse. From the nonlinear propagation 
the pressure in the vicinity of the bubble can be esti- 
mated to be in the range 40-60 kbar. 

For large enough P, the presence of shocks in the 
liquid results from the Rayleigh-Plesset dynamics for the 
bubble wall, independent of the state of motion of the 
gas inside the bubble. Comparisons by Wang et al 
(1999) between the strength of the measured pulse and 
that predicted by the Rayleigh-Plesset equation show 
that the strength of the wave in the liquid can be ac- 
counted for without including the effects of possible 
shocks in the gas. 

Another interesting effect of the emitted sound radia 
tion is that it influences measurements of the bubble 
radius by Mie scattering. Gompf and Pecha (2000; Pecha 
and Gompf, 2000) showed that in the last nanoseconds 
around the minimum radius most of the Mie scattering is 
by the highly compressed water around the bubble (see 
Fig. 22), not by the bubble surface itself. Neglecting this 
effect leads to an overestimate of the bubble-wall veloc- 
ity. Taking this effect into account, Gompf and Pecha 
(2000) found the bubble wall accelerates to about 
950 m/s, revising previously reported values of 
1200-1600 m/s by Weninger, Barber, and Putterman 
(1997; Putterman and Weninger, 2000). 


G. Bjerknes forces 


AIl of the calculations above assume that the center of 
the bubble is stationary in space. When neglecting vis- 
cous effects, the instantaneous force on the bubble is 


given by 
jl puds, 


where n is the outward normal vector to the bubble sur- 
face, and p — p,d is the pressure in the fluid. Multi- 
plication of Eq. (23) by 6, the unit vector in the direc- 
tion from the origin to the bubble position, gives the 
force component in that direction. Using Gauss's theo- 
rem and time averaging over a driving period, we obtain 
the (primary) Bjerknes force, first described by Bjerknes 
(1909), 


Fyj7 (BE puptie)=(— $ TRV). (24) 


To leading order, we can replace Vp by Vp(r=0,t) here. 
While both p and R are periodic, the product occurring 
in Eq. (24) does not, in general, average to zero. For the 
center of the bubble to be stationary, this force must 
vanish. For bubbles at a pressure minimum or maxi- 
mum, such as in the center of a flask in an SBSL experi- 
ment, Vp=0, and indeed Fy;=0. When the bubble is 
slightly off center, it depends on the relative phase of the 





Fhubble= (23) 
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pressure at the center and the bubble-radius dynamics if 
the net effect of F is to drive the bubble back to the 
center (stabilizing it), or to drive it further away. For 
linearly oscillating bubbles, it is easy to verify that 
bubbles whose resonance frequency fy is greater than 
the driving frequency f are attracted by pressure 
maxima (antinodes) and repelled by pressure minima 
(nodes). Bubbles with a smaller resonance than driving 
frequency show the opposite behavior. Indeed, for SBSL 
bubbles fo>f, and they are driven toward the pressure 
antinode at the center of the flask, where they are driven 
maximally. 

A subtle correction to these results ori 
small buoyancy force, 


pg (7« 
^J, V(0)dt, 


which also acts on the bubble (here g is the gravitational 
acceleration, T4=1/f the period of the driving, and V 
the bubble volume). This upward force must balance the 
downward component of the Bjerknes force so that the 
resulting equilibrium position is not in the center of the 
flask (z=0), but at (Matula et al., 1997) 


Ji V(t)dt 


ates in the 





Fou (25) 





Lau Gap A —— eo 
e f V(t)sin(wr)de 
where k, is the wave number of the standing pressure 
field along the direction of gravity. Experiments by 
Matula ef al. (1997) on zequi qualitatively agree with 
equation Eq. (26). However, the theoretical prediction 
seems too small by a factor of about 10. Matula (1999) 
gives evidence that the discrepancy could be connected 
with the back reaction of the bubble on the sound field. 

Note that both the acoustic and the buoyancy forces 
are fluctuating over one period, leading to small fluctua- 
tions of the equilibrium position as well. Aspherical, 
weaker bubble collapses and fainter light emission could 
be a consequence. Matula (2000) presented evidence 
that in microgravity, SBSL is somewhat stronger than for 
normal gravity, because the bubble collapse is more 
spherical. 

For small driving pressures, the position of an SBSL 
bubble is stabilized by the Bjerknes forces (see above). 
But sonoluminescing bubbles are strongly driven, which 
leads to variations in the phase shift between driving and 
bubble dynamics. As pointed out by Cordry (1995), 
Akhatov et al. (1997), Matula et al. (1997), and Matula 
(1999), for very large forcing pressure, Fy; can become 
repulsive, driving the bubble away from the center of 
the flask, rendering SBSL impossible. The calculations 
of Akhatov etal. (1997), Matula etal. (1997), and 
Matula (1999) demonstrate that this Bjerknes instability 
occurs above pressure amplitudes of P,~1.8 bars, al- 
ready above the upper threshold where single-bubble 
sonoluminescence usually occurs. Current experimental 
data appear to indicate that shape instabilities limit the 
upper threshold of sonoluminescence, which is discussed 
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FIG. 23. The difficulty in modeling SBSL. The bubble tem- 
perature T(t) is obtained from the radius dynamics R(t) (left), 
and the spectral radiance P,(r) is in turn deduced from the 
temperature. In contrast to R(r) and P, (f), the temperature 
cannot be measured directly. 


in detail below, It should be remarked, however, that 
those calculations neglect the back reaction of the bub- 
ble's pressure field on the bubble, as well as the effect of 
water vapor, and so might overestimate the Bjerknes 
threshold in some situations. 


Ill, THE BUBBLE INTERIOR 


One of the key problems in sonoluminescence re- 
search is that direct measurements of the state of matter 
inside the bubble are extremely difficult to perform. 
Practically all information about the conditions inside 
the bubble is obtained indirectly. One can measure and 
model the bubble dynamics and then use this as a basis 
for inferring the temperatures, pressures, etc. inside the 
bubble. Or, alternatively, one starts with observations of 
the light emission and uses the spectral information, the 
intensity, and the widths of the light pulses to deduce the 
conditions inside. These two approaches to modeling 
SBSL are sketched in Fig. 23. 

The information obtained in these two ways should 
obviously be consistent in a viable theory of sonolumi- 
nescence. If this consistency condition is fulfilled, how- 
ever, it is still not clear whether both the hydrodynamic 
model for the interior of the bubble and the model of 
the light emission are correct, as modeling errors could 
compensate each other. 

‘The most crucial variable of the bubble interior for 
which direct measurement is not possible is tempera- 
ture. As will be discussed in Sec. V, light emission is 
expected to depend sensitively on this quantity. In addi- 
tion, the contents of the bubble are a complicated func- 
tion of time, Even when starting out with a certain well- 
defined gas or gas mixture inside the bubble, processes 
of gas diffusion (Fyrillas and Szeri, 1994), gas rectifica- 
tion (Lohse et al., 1997), water-vapor condensation and 
evaporation (Moss er al., 1999), and chemical reactions 
(Yasui, 1997a; Storey and Szeri, 2000) lead to variations 
in composition, both within a cycle (time scales of mi- 
croseconds) and over many cycles (time scales of sec- 
onds). All properties of the matter inside the bubble 
(the equations of state, thermal diffusivity, viscosity, etc.) 
in turn depend on both gas composition and tempera- 
ture. Unfortunately, there are few solid data for these 
important dependencies under the extreme conditions 
of sonoluminescence, conditions not approached in any 
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other lab experiment, with the possible exception of 
shock tubes (Zel'dovich and Raizer, 1966). 

A quantitative understanding of single-bubble sonolu- 
minescence requires that each of these difficulties be ad- 
dressed step by step. To the present authors, one of the 
exciting features of modern research on single-bubble 
sonoluminescence is that it is a testing ground for how 
well mathematical models can deal with such a compli- 
cated situation, 

We shall organize our discussion of the state of matter 
in the bubble’s interior into two parts: in this section, we 
shall describe the fluid mechanics of the bubble's interior 
and the various attempts to use it to infer bubble tem- 
peratures at collapse. The goal of this section is to un- 
derstand both the maximum temperature and the com- 
position of the bubble. These pieces of information can 
then be fed directly into a model of the light emission, a 
discussion of which will be deferred to Sec. V. Although 
we have chosen for reasons of presentation to break up 
our discussion into these two parts, it should be empha- 
sized that the research is not at all independent: Models 
of the light emission critically depend on the tempera- 
tures predicted from hydrodynamic calculations, while 
more sophisticated models of gas dynamics have in turn 
been developed in order to explain properties of the 
light emission. 

In Sec. ILA, we shall summarize work in which the 
full compressible gas-dynamical equations inside the 
bubble are solved. Over the years (spurred on by more 
detailed information about the light emission) the mod- 
els have incorporated more and more physical effects. 
The most important modifications of the earliest models 
concern the inclusion of dissipative and transport pro- 
cesses, in particular those involving water vapor inside 
the bubble. 

An alternative approach assumes a (nearly) uniform 
bubble interior and thus avoids the solution of the 
Navier-Stokes equations. While less accurate, such mod- 
eling is computationally inexpensive and allows for the 
calculation of temperatures for many more parameter 
combinations. Several variants of this simpler approach 
are treated in Sec. IIB. 

We briefly mention here that molecular dynamics is a 
third possibility for modeling the bubble interior. Fol- 
lowing the motion of the ~10! molecules or atoms in a 
SBSL bubble is beyond the capability of present-day 
computers, so that simulations have had to be conducted 
with a far smaller number of quasiparticles (Matsumoto 
et al., 2000; Metten and Lauterborn, 2000), limiting the 
prospect for quantitative comparison with experiment. 
One of the main problems of this type of approach is 
that, due to the reduced number of particles, the number 
of particle collisions is drastically lower than in reality, 
and therefore it is hard to achieve thermal equilibrium. 


A. Full gas dynami 





in the bubble 


Assuming local equilibrium, the motion of the gas in- 
side the bubble can be described by the Navier-Stokes 
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equations and the equations of energy and mass conser- 
vation (Landau and Lifshitz, 1987), 


Op_t o pgvi)=0, (27) 
(poi) + ONP edit pv; Ty) =0, (28) 
AE + a{(E +p g)v4]— 4(vj7y)— 4(K 9,7) =0. (29) 


Velocity components inside the gas are denoted v;; py 
and p, are the gas density and pressure, while E=p,e 
+ p,v7/2 is the total energy density, with e the internal 
energy per unit mass. T is the gas temperature and K, 
its thermal conductivity. The viscous stress tensor is 
given by 





(30) 


where 7 is the gas viscosity and the effects of the sec- 
ond viscosity have been neglected. These equations have 
to be completed with an equation of state, connecting 
density, pressure, and temperature. Depending on the 
degree of sophistication, it might also be necessary to 
include the effects of vibrational excitation, dissociation, 
ionization, and intermolecular potentials. In addition, 
the material parameters K,,7, themselves depend on 
temperature and pressure, 

Finally, one must impose boundary conditions at the 
moving bubble wall r= R(t). These can be dealt with in 
two ways: either the velocity at the bubble wall is taken 
to be that predicted by the Rayleigh-Plesset equation 
v,(rf)=R(t), or alternatively one could solve the full 
fluid-dynamical equations also in the surrounding water. 
For completeness, boundary conditions for both mass 
and heat exchange must also be formulated. 

This problem has been attacked with an increasing 
level of detail, motivated by advances in experiments. 
We review these efforts in roughly chronological order, 
grouping them into 
* inviscid models (Wu and Roberts, 1993; Moss et al., 

1994; Kondic et aL, 1995; Chu and Leung, 1997); 

* dissipative models (Vuong and Szeri, 1996; Moss et al., 

1997; Cheng et al., 1998); 

* dissipative models including phase change, in particu- 

lar that of water vapor (Storey and Szeri, 2000). 


Tjj 7 Nel Air ð; 





[D 


All of these approaches treat the bubble as spherically 
symmetric. 


1. Inviscid models 


Motivated by the measurements of Barber et al 
(1992) indicating that the width of the SBSL light pulse 
was shorter than 50 ps, early theories focused on the 
idea that shocks were important for single-bubble 
sonoluminescence (Greenspan and Nadim, 1993; Wu 
and Roberts, 1993; Moss et aL, 1994). Shock focusing 
provides a natural mechanism for producing both ex- 
tremely high temperatures and a pulse width that is 
much smaller than the time the bubble stays fully col- 
lapsed. 
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One of the first numerical solutions of the (spherical) 
gas-dynamical equations driven by the Rayleigh-Plesset 
dynamics was done by Wu and Roberts (1993). The 
most important approximations of this work were (i) vis- 
cosity and thermal diffusion are assumed negligible, (ii) 
no heat or mass exchange takes place between the 
bubble and the surrounding water, and (iii) a van der 
Waals equation of state with a polytropic exponent I 
=7/5 is assumed throughout the collapse. For a Ro 
.5 um bubble driven at P,=1275 atm and f 
6.5 kHz, Wu and Roberts (1993, 1994) found a 
spherical shock wave launching from the wall, focusing 
to the center, and reflecting outward again. Tempera- 
tures in excess of 10" K and light pulses of 1.2-ps dura- 
tion were predicted. 

The high temperatures and short pulse widths can be 
understood from the classical analytical solution of the 
equations of gas dynamics in an imploding sphere by 
Guderley (1942; see also Landau and Lifshitz, 1987). 
Guderley neglected viscosity and thermal diffusion, and 
assumed an ideal gas equation of state. His result shows 
that a converging shock wave focuses to the center of 
the sphere with a radius 

Rj(t)- (£*—0*, G1 
with an exponent a~0.6884 for l'—5/3 and a~0.7172 
for '=7/5. Here t* represents the time at which the 
shock reaches the bubble center. In the case of T 5/3, 
the temperature at the center of the shock diverges as 
R,^, with B~0.9053. When the shock reaches the 
bubble center, the temperature is mathematically infi- 
nite. With a van der Waals equation of state the same 
singularity (31) with a slightly different exponent occurs; 
Wu and Roberts (1994) show that their simulations con- 
verge onto this solution. Similar calculations were per- 
formed by Moss et al. (1994) and Kondic et al. (1995). 
Moss et al. (1994) used a more sophisticated equation of 
state for air inside the bubble, limiting the maximum 
temperature through the energy-consuming processes of 
dissociation and ionization. They also solved the full 
equations for the motion of the water around the 
bubble. 

The unphysical divergence of temperature from Gud- 
erley’s solution must, of course, be avoided in reality. 
Evans (1996) noted that the converging spherical shock 
wave will be susceptible to instabilities in its shape 
(much like the bubble shape instabilities discussed in 
detail in Sec. IV). Evans calculated a relatively slow di- 
vergence of the relative size of the perturbations, ap- 
proximately scaling as 6R,/R,%r~*° with yy«1. He 
concluded that very high temperatures would still be 
possible inside the bubble, but that the shape instability 
does set a limit to the degree of energy focusing. 

Dissipative processes also are capable of stopping the 
temperature divergence of Guderleys solution. Both 
heat and temperature diffusion potentially disrupt the 
formation of a shock wave because they counteract the 
steepening gradients at the shock front. The question 
here is primarily one of time scales: when a converging 
pressure wave travels towards the bubble center, can it 








444 Brenner, Hilgenfeldt, and Lohse: 


steepen sufficiently quickly to develop into a (strong) 
shock? Or is dissipation so strong that a front never de- 
velops before the center is reached? These ideas were 
first touched upon by Vuong and Szeri (1996, described 
in some detail below), who demonstrated the difficulty 
of generating shocks in a model where the transport co- 
efficients depend on gas density. 

In the years since the shock picture was proposed, 
there has been direct experimental evidence arguing 
against the likelihood of shock formation. Following the 
work of Ohl et al. (1998; Ohl, 2000) on laser-induced 
cavitating bubbles, which also emit light (“single- 
cavitation bubble luminescence,” SCBL), Baghdassarian 
etal. (1999, 2001) found that highly shape-distorted 
bubbles are still able to give off considerable lumines- 
cence. In such bubbles a spherical shock wave cannot 
exist. Evans (1996) showed that an asphericity of only 
5*6 in the bubble wall is sufficient to disrupt the energy- 
focusing power of a shock. 


2. Dissipative models 


Kondic et al. (1995) realized the necessity of including 
dissipation in the gas dynamics, giving some estimates 
for heat conduction and pointing out radiative transfer 
as an energy-loss mechanism (the latter process turns 
out to have a negligible effect on the total-energy bal- 
ance of the bubble). Going beyond estimates, Vuong 
and Szeri (1996) included thermal and viscous dissipa- 
tion in solving the equations of motion of the gas.“ To do 
this properly, it is essential to understand how the mate- 
rial properties (thermal conductivity and viscosity) de- 
pend on the temperature and on the pressure. At the 
extreme conditions achieved in the bubble, those depen- 
dencies are not known, and one has to try either to de- 
rive them from first principles or to extrapolate approxi- 
mate relationships known from measurements at lower 
pressures and temperatures, Vuong and Szeri (1996) and 
also Cheng et al, (1998) do the latter and assume the 
linear relation 


K 





r G2» 


EE i 
(To) 7 
between heat conductivity and temperature; the pres- 
sure dependence is neglected (see also Kamath et al., 
1993; Yasui, 1995). 

Vuong and Szeri (1996) and Cheng er al, (1998) also 
included the heat flux in the water by coupling to the 
equation 


AT, Eu, (33) 





333,7, 
iré) 





We remark that viscous effects had been previously consid- 
ered in bubble dynamics, primarily by Prosperetti (1991) and 
Kamath et al. (1993). As described in Sec. II, these authors 
also drew the important conclusion that the heat transport be- 
tween the interior of the bubble and the liquid results in a 
bubble with an isothermal boundary for most of the forcing 
cycle. 
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Single-bubble sonoluminescence 


TABLE I. Maximum temperatures achieved in a Ry=4.5 um 
pure argon bubble driven at f=26.5 kHz. A dissipative gas 
dynamics mode! was used in five runs using different driving 
pressure amplitudes P, and gas species. From Vuong and Szeri 
(1999). 











Run Gas P, (atm) Tax (K) 
1 Ar 11 20000 
n Ar 12 52000 
nm Ar 13 118000 
Iv He 13 43000 
v Xe 13 202000 








for r>R(t); Ty(r.t) and u(r,) are the temperature 
and the spherical velocity field in the water, and x; the 
thermal diffusivity. The radius R(f) is given by the 
Rayleigh-Plesset equation. At the boundary r=R(t), 
the heat flux out of the bubble and into the fluid must 
match, and the same is generally assumed for the tem- 
peratures of liquid and gas. 

The main result of Vuong and Szeri's paper shows that 
no shocks occur in argon bubbles. Though some wavy 
structures can be seen, they do not steepen to a shock, as 
can be seen in Fig. 16. The temperature profile in the 
collapsing bubble is not dramatically peaked near the 
center, but rather shows slow variations for most of the 
radius, with a strong decline to the ambient water tem- 
perature near the wall. Nevertheless, through the nearly 
adiabatic compression of the gas in the bubble, very high 
temperatures above 10° K are achieved in the simula- 
tions. The precise values depend on the type of gas em- 
ployed and on the control parameters, as seen in Table I. 
However, all values are dramatically lower than in the 
calculations without viscosity and thermal conductivity 
carried out by Wu and Roberts (1993), Moss et al. 
(1994), and Kondic et al. (1995). 

Cheng et al, (1998) confirmed the results of Vuong 
and Szeri (1996) for argon and in addition repeated the 
calculation for pure nitrogen gas. They found that, for 
nitrogen with its smaller polytropic exponent T'=7/5, 
shocks can develop for strong enough forcing, but they 
are limited to a tiny region around the center of the 
bubble, and the peak power (assuming blackbody radia- 
tion) is much less than that of argon. In particular, the 
width of the power pulse is in the range <5 ps, much 
smaller than observed by Gompf et al. (1997). Earlier, 
Moss etal. (1997) had arrived at similar conclusions 
about shock occurrence when neglecting viscosity and 
normal heat conduction, but including heat conduction 
of the ionic and electronic constituents of the ionized gas 
generated in the bubble. Their finding that “it is more 
difficult to generate a shock in Ar than in N;" thus 
seems to be true no matter what the dissipation process 
involved. 

One of the main reasons that shocks are suppressed in 
argon bubbles is that the strong heating already supplied 
through adiabatic compression (with T =5/3) results [via 
Eq. (32)] in an enormous thermal conduction, which lev- 
els temperature gradients in the bubble (Vuong and 
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Szeri, 1996), preventing the steepening of the wavy dis- 
turbances into a shock. In nitrogen (T — 7/5) this effect is 
less pronounced and shocks can develop for strong col- 
lapses (Moss ef al., 1997; Cheng et al., 1998). 

A more general statement on shock suppression is 
made by Lin and Szeri (2001), who attribute the diffi- 
culty of observing any sign of spherically converging 
shocks to the presence of adverse entropy gradients. In 
the case of sonoluminescing bubbles, the sound speed in 
the gas increases towards the center of the bubble, de- 
laying and weakening shock formation, 

‘Vuong and Szeri (1996) also analyzed the dependence 
of the peak temperature on the type of inert gas em- 
ployed. Both MBSL and SBSL intensity increase with 
increasing atomic weight of the noble gas, from He to 
Xe. Often, the lower thermal conductivity of xenon 
compared to helium has been favored as the main rea- 
son for higher temperatures and more light emission in 
Xe (Verral and Sehgal, 1988; Greenspan and Nadim, 
1993; Hickling, 1993, 1994). However, the inert gases 
also vary in diffusivity, ionization potential, and many 
other physical properties, and it is important to deter- 
mine which are crucial in the context of sonolumines- 
cence. Vuong and Szeri find that contributing to the 
higher peak central temperatures in their numerical 
simulations with xenon bubbles (see Table I) are not 
only the effect of lower thermal conductivity of Xe, but 
in addition another mechanism also due to the larger 
mass of this inert gas: waves in heavier gases carry more 
momentum and travel more slowly, leading to stronger 
compression lasting for a longer time. 

A further noteworthy complication of the gas dynam- 
ics inside the bubble originates from a weak mass segre- 
gation effect inside an oscillating bubble. Storey and 
Szeri (1999) have shown that for a mixture of two gases 
with different masses the heavier one has a slightly 
higher concentration at the bubble's edge, whereas the 
lighter one is concentrated towards the bubble’s center. 
This segregation is caused by the mass dependence of 
the diffusion coefficients, If the bubble's interior were a 
mixture of inert gases (the example analyzed by Storey 
and Szeri, 1999), this effect would have little relevance. 
But when water vapor is considered (see Sec. II. A.3), 
mass segregation can have important effects when com- 
paring helium bubbles with xenon bubbles (Storey and 
Szeri, 1999; Yasui, 2001): The water-vapor concentration 
in helium bubbles is smaller than that in xenon, since in 
helium the mass diffusivity is higher. Such a distribution 
of vapor tends to diminish the differences in peak tem- 
peratures between Xe and He (see also Sec. V.F.3). 





3. Dissipative models including water vapor 


Over time, dissipative models for the bubble’s interior 
have become more and more complex. An important 
step was the realization that the water vapor inside the 
SBSL bubble plays a crucial role in regulating the heat 
transfer to and from the bubble (Kamath et al., 1993; 
‘Yasui, 1997a; Colussi and Hoffmann, 1999; Moss et al., 
1999; Storey and Szeri, 2000; Toegel, Gompf, et al., 
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2000). This effect of liquid vapor has been well known in 
mulübubble sonoluminescence for decades (Jarman, 
1959; Flint and Suslick, 1989, 1991a, 1991b). In addition, 
the water vapor invading a collapsing bubble will un- 
dergo chemical reactions that also change the tempera- 
ture (Kamath er al., 1993; Yasui, 1997a; Gong and Hart, 
1998; Storey and Szeri, 2000). 

Initially, attempts to incorporate water vapor were re- 
stricted to simplified models (e.g., Yasui, 1997). To our 
knowledge, the first full numerical simulation of a 
sonoluminescing bubble with water vapor was per- 
formed by Moss et al. (1999), building on their earlier 
work (Moss et al., 1997), which will be discussed in detail 
in Sec. V. Moss et al. (1999) kept the amount of water 
vapor fixed and uniform during the cycle and did not 
take chemical reactions into consideration. Thermal 
conduction is included for the neutral gas inside the 
bubble, for ions, for electrons, and for the water around. 
the bubble. Due to the high pressures and densities in- 
side the collapsed bubble, ions, electrons, and neutral 
particles equilibrate on a time scale much shorter than 
the other relevant time scales, so that a single tempera- 
ture field T(r,) is sufficient (rather than having differ- 
ent temperatures for the ions, electrons, etc., as in Moss 
et al., 1997). The equation of state was obtained from a 
combination of data and theoretical work by Young and 
Corey (1995). 

The fixed water-vapor content was fit to one data set 
of Gaitan and Holt (1999) on the light intensity for 
known forcing pressure, frequency, and compression ra- 
tio (Rmax/Ro), and then extrapolated to other forcing 
pressures, frequencies, and compression ratios. For stan- 
dard conditions at T,=20°C and f~20 kHz, the fixed 
water-vapor concentration was found to be 30-40%. At 
collapse, the transfer of vapor out of the bubble is too 
slow to keep up with the bubble-wall motion, so that 
vapor becomes trapped inside the bubble (Moss et al., 
1999), This point of view was later confirmed and inves- 
tigated further by Storey and Szeri (2000; see below). 
For the extreme conditions achieved at collapse, fixing 
the amount of water inside the bubble is not as crude an 
approximation as it may sound. 

The main result of Moss et al. (1999) was that the in- 
clusion of a sufficient quantity of water vapor leads to a 
smaller (effective) adiabatic exponent y. Heating is thus 
reduced and shock waves can occur close to the center, 
in agreement with earlier work by Cheng et al. (1998) 
for nitrogen bubbles. The shock-wave heating, however, 
is very inhomogeneous. Although the center of the 
bubble can be very hot (up to 130000 K in the examples 
given), the total radiating volume within the shock front 
is so small that there is far less optical emission than for 
adiabatic heating with a larger l' which does not lead to 
shock waves (Vuong and Szeri, 1996; Cheng et al., 1998). 
A large part of the compressional energy of the collapse 
is found to go into internal degrees of freedom rather 
than into heating of the bubble. Moss et al. (1999) argue 
that this is one reason why sonoluminescing bubbles are 
brighter at low water temperature (Barber et al., 1994), 
where the water vapor pressure is less. 
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The model of Moss et al. (1999) does rather well in 
explaining the trends in the experimental data of Gaitan 
and Holt (1999). Moss et al. in particular stress that a 
pure argon bubble leads to very different heating char- 
acteristics (and thus light emission) when compared to 
an argon bubble with some water vapor in it; only the 
latter is consistent with the data. Calculations without 
water vapor would predict that a strongly driven bubble 
in the Gaitan-Holt measurements (Ro=6.0 um, R max 
=64 um) should emit 100 times as much light as a 
weakly driven bubble (Ro=2.1 um, R mas=30 um). The 
experiments actually show an intensity ratio of only 
~10, in agreement with simulations including water va- 
por. The calculated results are extremely sensitively de- 
pendent on the experimental conditions and 
parameters—or, likewise, the modeling assumptions. A. 
particularly striking example is the dependence on the 
ambient pressure Po. Reducing it from P, 1 bar to 0.99 
bar can lead to 1.6 times more light, according to the 
model. The reason for the sensitivity lies not so much in 
the hydrodynamics, but in the extremely strong depen- 
dence of the light intensity on the temperature achieved 
in the bubble, as we shall discuss in detail in Sec. V. In 
the above example the temperature of the light-emitting 
region in the bubble increases by just 5%. 

‘The restriction of Moss et al. to constant water-vapor 
content is relaxed in a remarkable paper by Storey and 
Szeri (2000), who build on their earlier work on mass 
segregation in gas bubbles containing a mixture of dif- 
ferent inert gases (Storey and Szeri, 1999). The transport 
parameters (thermal diffusion, mass diffusion, and vis- 
Cosity) are calculated from equations based on 
Chapman-Enskog theory (Hirschfelder eral, 1954); 
thermal dependencies are considered as far as they are 
known. The equations of state are of the Soave-Redlich- 
Kwong type (Gardiner, 1984; Reid et aL, 1987). Storey 
and Szeri (2000) argue that although the basic physical 
mechanisms they uncover are robust, the exact numbers 
do depend on (unknown) details of the equation of 
state. 

Evaporation and condensation are modeled using a 
kinetic theory (Carey, 1992). The rate of mass transport 
(per unit area) at which water molecules pervade the 
interface is «(pio pau). Where Pro is the partial 
pressure of water and p,,, the saturation pressure at the 
temperature of the interface. Not all water molecules 
that hit the wall actually stick to it, only a fraction v, 
thereof. This accommodation coefficient is chosen to be 
0.4, following Yasui (1997a) and Eames et al. (1997). 

As an example, Storey and Szeri (2000) study a 
bubble initially consisting of argon, with Ro=4.5 um 
and driven at P,=1.2 bars and f=26.5 kHz: Including 
water vapor, the maximum temperature is reduced from 
20900 K (cf. Vuong and Szeri, 1996) to 9700 K, due to 
the lower T. No shock waves are observed under these 
conditions. The amount of water in the bubble is not 
constant at all, as can be seen from Fig. 24: A large 
amount of water evaporates into the bubble during the 
main expansion when the pressure is low; at this stage 
vapor is in equilibrium with the vapor pressure. At 
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FIG. 24. Importance of water vapor exchange between bubble 
and fluid: (a) Bubble dynamics and (b) number of water mol- 
ecules in the bubble as a function of time for an (initially) Ro 
Sum argon bubble driven at P,-12atm and f 
26.5 kHz. Note the comparison in (b) to the constant num- 
ber of argon atoms. From Storey and Szeri (2000). 





bubble maximum, about 90% of the bubble contents is 
water, Vapor condenses out of the bubble again at col- 
lapse, but not completely, since the time scale of the 
collapse z;,,7 R/|R| becomes much smaller than the 
time scale for the transport of water vapor out of the 
bubble. The vapor transport is a two-step process, con- 
sisting of diffusion to the wall and condensation, so that 
it involves two time scales, one for vapor diffusion in the 
bubble, 
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and one for condensation at the wall, 
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where c,= JT Popga is the sound velocity of the gas in 
the initial state, and My,o and Mo are the molecular 
masses of water and the initial bubble content, respec- 
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FIG. 25. Profiles of the molar fraction of water vapor vs radius 
for several times prior to the moment of minimum radius 
(time=0, leftmost profile). In order of increasing radius, the 
Other curves correspond to times prior to collapse of 
1 —5, 15, —14, —153, — 299, —455, and — 843 ns. Calculation 
for a Ro=45 um argon bubble driven at P,=1.2 atm and f 
26,5 kHz. From Storey and Szeri (2000). 





tively. Figure 25 shows the profiles of the water-vapor 
molar fraction inside the bubble for various times: In the 
early part of the collapse, 745, Taif» Teon» SO that the 
profiles are uniform. The vapor concentration is de- 
creasing, as water has sufficient time to escape. Later, 
when 7,,47,//,7,,, (the three lowest profiles of Fig. 
25, starting from 15ns before the collapse), the water- 
vapor profile is essentially “frozen” within the bubble. 
Therefore, the bubble consists of ~14% water vapor 
(on a molar basis) through the collapse; the vapor is 
trapped in the bubble’s interior. 

Comparing the time scales 7, and Teon, Storey and 
Szeri (2000) found that the vapor transport was always 
diffusion limited (r7 Teon). Models with unrealisti- 
cally low accommodation coefficients a,~0.001, like 
that of Colussi et al. (1998), on the other hand, could 
well be condensation limited. 

Rather remarkably, when applying their numerical 
model to the more strongly forced bubbles of Moss er al. 
(1999; see above), Storey and Szeri (2000) find average 
water-vapor concentrations close to what Moss et al. ob- 
tained from fitting their numerical results to the Gaitan 
and Holt (1999) data mentioned above, a strong indica- 
tion that these numbers are in the right ballpark. 

In the second part of their paper, Storey and Szeri 
(2000) consider chemical reactions of the water vapor. 
The reaction scheme employed originates from Maas 
and Warnatz (1988) and includes 19 forward and reverse 
elementary reactions of the nine species Ar, H, Hp, O, 
Oz, OH, HO;, H,O, and H,O, (argon only enters the 
reactions as a third body). The reaction rates were taken 
from Gardiner (1984) and in the high-pressure limit, for 
which there is “considerable uncertainty,” from Bow- 
man et al. (1999). 

With the chemical reactions included, the maximum 
temperature for the standard case (a Rọ=4.5 um argon 
bubble driven at P,=12 bars and f=26.5 kHz) de- 
creases from 9700 K (with water vapor, but without re- 
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FIG. 26. Dependence of the bubble center temperature on the 
compression ratio Rmax/Rmin. O, calculations without water 
vapor; *, with water vapor, but no reactions; +, with water 
vapor, including reactions. Note the logarithmic scale on the 
axes, From Storey and Szeri (2000). 


actions) to 7000K. Most of the reactions are endother- 
mic and therefore eat up the collapse energy, which 
otherwise would be used for heating. The reaction path- 
ways can also be interpreted as additional degrees of 
freedom, further lowering the effective polytropic expo- 
nent. Note that the time scales of the chemical reactions 
are so fast (because of the high densities) that thermo- 
chemical equilibrium should prevail in the collapse up to 
the point of minimum radius. 

How does the bubble temperature depend on the 
forcing pressure? Figure 26 shows results by Storey and 
Szeri (2000) for the temperature at the center as a func- 
tion of the compression ratio R max/Rmin: No shocks oc- 
curred in these calculations, so that the center tempera- 
ture is a fairly good indication of the overall bubble 
temperature. For other parameter combinations (not 
shown in Fig. 26), Storey and Szeri (2000) do find shocks 
and very high maximum temperatures, but only right at 
the center, corresponding to ~0.1% of the total bubble 
volume. A suitably averaged temperature, representa- 
tive of the actual energy content of the bubble, will be 
very similar to those shown in Fig. 26. It is noteworthy 
that for these shocks the presence of water vapor is nec- 
essary, as expected from the work of Cheng ef al. (1998) 
and Moss er al. (1999). For a low compression ratio, 
hardly any vapor invades the bubbles of Fig. 26 and the 
three cases (without water vapor, with nonreacting water 
vapor, and with reacting water vapor) give roughly the 
same maximum temperature. The remarkable feature, 
however, is that the temperature levels off for large com- 
pression ratios when water vapor is taken into account, 
asymptoting to ~10000 K without chemical reactions 
and ~7000 K with chemical reactions included. Thus the 
initial power-law increase of the maximum temperature 
with expansion ratio does not extend into the regime of 
SBSL temperatures: the effect of the larger expansion 
ratio of stronger forced bubbles is compensated by the 
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presence of more water vapor, thus limiting the heating 
at collapse. 

This finding has important implications for the quest 
for “upscaled” sonoluminescence: To learn the limits of 
energy focusing and light emission, researchers have 
tried to induce more violent bubble collapses by, for ex- 
ample, applying nonsinusoidal driving pressures (Holz- 
fuss, Riiggeberg, and Mettin, 1998) or reduced driving 
frequencies as suggested by Hilgenfeldt and Lohse 
(1999) and Toegel, Gompf, et al. (2000). But the larger 
expansion ratios achieved with these techniques lead to 
a larger water-vapor content of the bubble, which again 
limits the heating at collapse, A theoretical study by 
Toegel, Gompf, et al. (2000) shows that these two effects 
roughly cancel each other, essentially leading to the 
same temperature and the same amount of light at lower 
driving frequencies for otherwise identical bubble pa- 
rameters, A residual upscaling effect may still be ob- 
served (as in Barber and Putterman, 1991) due to the 
possibility of stabilizing bubbles with larger Ro at lower 
driving frequencies (see Sec. IV.D). 


B. Simple models 


‘The previous section has shown that when the com- 
plex interplay of physical effects inside a sonoluminesc- 
ing bubble is meticulously included, spatial inhomogene- 
ities inside the bubble are not very pronounced. 
‘Therefore it seems reasonable to try to approximate the 
bubble’s interior by spatially constant, time-dependent 
pressure p(t) and temperature T(t). Clearly, such an 
approximation is too crude to capture some of the afore- 
mentioned effects such as mass segregation (Storey and 
Szeri, 1999). However, such modeling simplifications are 
extremely useful tools in exploring the phase space of 
sonoluminescence (see Sec. IV). Such a scanning of 
phase space is at present not possible for the complete 
models discussed in the previous section. 

Simple models assuming a uniform bubble interior 
have been developed with increasing detail; here, we 
shall discuss two types. Just as in the more elaborate 
models, one can either completely neglect heat and mass 
transfer to and from the bubbles exterior (Sec. IILB.1) 
or try to embody these effects (Sec. IILB.2). 


1. Homogeneous van der Waals gas without heat and mass 
exchange 


The simplest model is to assume an adiabatic equation 
of state for the bubble interior (Löfstedt eral, 1993; 
Barber et al., 1997), 
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FIG. 27. Temperature inside the SBSL bubble within the sim- 
plified model of Hilgenfeldt et al. (19992, 199b). (a) Bubble 
dynamics and (b) temperature dynamics of Ro=5 um bubble 
driven at 1.3-atm pressure amplitude at /—20 kHz. The tem- 
perature is calculated via Eq. (38); the full width at half maxi- 
mum of the asymmetric peak is ~1 ns. 





15 20 


equivalent to Eq. (10) when replacing y by —C,/C, 
the ratio of the specific heats. Equation (36) supple- 
ments the Rayleigh-Plesset equation and permits its so- 
lution. 

One obvious problem with Eqs. (36) and (37) is that 
an adiabatic bubble motion is assumed with no heat ex- 
changed between the bubble and the exterior. As 
pointed out in Sec. II, Péclet number estimates via Eq. 
(11) show that there is almost unrestricted heat ex- 
change for most of the oscillation cycle, and the motion 
is isothermal. Only at bubble collapse is the Péclet num- 
ber Pe larger than 1. This means that most of the time 
the ratio of the specific heat I’ has to be replaced by an 
isothermal exponent y=1. Near the cavitation collapse 
there is a change in the adiabatic value y—T. Roughly 
speaking, this transition will occur when the Péclet num- 
ber is of order unity. 

In linear approximation, Prosperetti (1977c) derived a 
transition function »(Pe(f)) from the isothermal behav- 
ior y-1 to the adiabatic behavior y-T. Hilgenfeldt 
et al. (1999a, 1999b) have employed this approach to cal- 
culate approximately the temperatures and pressures in 
SBSL bubbles. The Rayleigh-Plesset equation is then 
supplemented by a differential version of Eq. (37) with 
variable y(Pe(0). 








3RR 
T- Ly Pelt))-1] gx T- (P= Two) s / R2. 

(38) 
The last term contains the ambient water temperature 
Tyo and the thermal diffusivity of the gas y, , which at 
high densities is assumed to follow Chapman-Enskog 
theory (Hirschfelder et al, 1954). Equation (38), to- 
gether with the Rayleigh-Plesset equation, gives a 
simple model for bubble radius and bubble temperature. 
The peak temperatures obtained—see Fig. 27—are com- 
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parable to bubble temperatures resulting from the more 
sophisticated model by Storey and Szeri (2000). The lat- 
ter publication reports somewhat lower temperatures, 
though, because of the inclusion of heat loss and mass 
exchange in the model, to which we turn in the following 
section, 


2. Homogeneous van der Waals gas with heat and mass 
exchange 


In more sophisticated variants of the model, the heat 
and mass exchange between the bubble and its exterior 
are explicitly taken into account. One of the first models 
of this type was conceived by Yasui (1997a). The ingre- 
dients of his approach are (i) a Rayleigh-Plesset-type 
equation for the bubble radius with a van der Waals gas 
inside; (ii) a spatially homogeneous bubble interior with 
time-dependent pressure and temperature (the material 
constants inside the bubble such as the thermal conduc- 
tivity are modeled as a function of the temperature); (iii) 
water-vapor exchange with the bubble's exterior due to 
condensation and evaporation; (iv) heat exchange with 
the bubbles exterior modeled by an energy flux depend- 
ing on the compression, the temperature gradient, and 
condensation/evaporation; (v) a thin layer of water 
around the bubble that can be heated; (vi) 25 chemical 
reactions of the water vapor, following Kamath et al. 
(1993) and using data from Baulch et al. (1972-1976). It 
is crucial that the net effect of these reactions be con- 
sumption of thermal energy, i.e., they must be endother- 
mic. 

The central result of Yasui (1997a) is displayed in Fig. 
28: At bubble maximum, the bubble consists nearly ex- 
clusively of water vapor. Even at collapse, the bubble 
still retains some of the water (~1% of the total bubble 
contents). Because of the invading water vapor and the 
endothermic chemical reactions, the maximum tempera- 
ture inside the bubble at typical control parameters is 
only around 10000 K, in agreement with the complete 
model of Storey and Szeri (2000). 

Models like that of Yasui (1997a) are useful for exam- 
ining the energy balance at collapse. For standard pa- 
rameters of single-bubble sonoluminescence (see Fig. 
28), one finds in the last 120 ps before collapse a reduc- 
tion of the thermal energy by 1.4 nJ through chemical 
reactions and by 0.6 nJ through thermal conduction. The 
loss through photon emission is negligible, only about 
0.2 pl. 

Yasui (19974) assumed that the transport of mass 
through the boundary layer was condensation limited, 
rather than diffusion limited (mass diffusion was not ex- 
plicitly modeled, and thus assumed instantaneous). 
However, Storey and Szeri (2000) later showed that 
transport under SBSL conditions is diffusion limited 
(see Sec. II.A.3). Therefore Toegel, Gompf, et al. (2000; 
Toegel et al., 2002) took the opposite approach and de- 
veloped a simple diffusion-limited model for water- 
vapor exchange between bubble and liquid, using a 
boundary layer approximation. The diffusive change in 
the number of water molecules over time is then 
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FIG. 28. Bubble temperature and number of molecules with 
the model of Yasui (19973). The parameters are Ro=5 um, 
P.,7 135 bars, f=20.6 kHz, P 1 atm, and a water tempera- 
ture of Ty=20°C. From Yasui (19972). 
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Ni,o=47R?D. à,n|,- ac 47R?D. Uv" (39) 
where no corresponds to the equilibrium density of va- 
por molecules at the wall and n is their actual concen- 
tration. The diffusion length is obtained through dimen- 
sional analysis as /,;,7 min[(RD/R)'^, R/ v], where D is 
the gas diffusion constant. The cutoff prevents the 
boundary layer from becoming unphysically large. 

The heat flux is treated in complete analogy to the 
flux of water vapor by Toegel, Gompf, er al. (2000; Toe- 
gel et al., 2002), i-e., 








O= 4a Rg (40) 
n 

where Q is the heat content of the bubble, Tyo the equi- 
librium (ambient) temperature, and the thermal diffu- 
sion length is f,4=min[(RYmix/R)"2,R/ s] with the ther- 
mal diffusivity x,;, of the gas mixture, Yasui (1997a) 
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uses a boundary layer as well, but in addition assumes a 
temperature jump between the outer edge of this 
boundary layer and the liquid. This temperature slip is 
more usually associated with low-density systems and 
does not appear in other work like that of Toegel, 
Gompf, et al. (2000) or Storey and Szeri (2000). 

Chemical reactions can also be included in a boundary 
layer model (Toegel et al., 2002). The most important 
endothermic process is 


H,0+5.1 eVOH+H, (41) 


whose inclusion already shows the essential effects of 
more sophisticated reaction schemes. Within the ordi- 
nary differential equations (ODE) formalism of Toegel 
et al. (2002), Eq. (41) yields reaction rates in good ac- 
cord with those of Storey and Szeri (2000). 

Together with the first law of thermodynamics and a 
van der Waals equation of state, the above formalism 
yields another ODE for the temperature inside the 
bubble, 

ðE 

CoP =O- pV +h Ni o > D (42) 
where the sum is over all species X= Ar, HO, OH, and 
H. The derivatives 2E/ƏN y as well as the heat capacity 
C, take into account rotational and vibrational degrees 
of freedom in the various molecules and are thus depen- 
dent on temperature. Therefore Eq. (42) is an implicite- 
quation for T. Finally, h,, is the enthalpy of water mol- 
ecules near the (cold) bubble wall. 

Equation (42) provides closure of the model together 
with a Rayleigh-Plesset equation variant. Similar equa- 
tions were discussed by Yasui (1997a), though the ex- 
plicit form of the terms differs as indicated above, and 
the diffusion-limited character of the transport was not 
taken into account. Such a set of four first-order equa- 
tions for R(t), R(t), Nuo(t), and T(t) can be solved 
for various physical parameters, such as forcing pressure 
P,, ambient radius (without water vapor) Ro, water 
temperature Tyo, and driving frequency f, without 
much numerical effort. The scheme can also be easily 
extended to more reactions. The results of Toegel, 
Gompf, et al. (2000) agree well with the full simulations 
of Storey and Szeri (2000). 

Most models that include chemical reactions inside 
the bubble—whether based on partial differential equa- 
tions like that of Storey and Szeri (2000) or on ODE like 
those of Kamath eral (1993) and Storey and Szeri 
(2001)—seem to underestimate the temperature inside 
the bubble. With the exception of the Yasui (19972) 
model (in which water vapor is not prominent because 
its transport is very fast), all of these models seem to 
imply temperatures substantially below 10000 K in the 
SBSL regime. Assuming a thermal lightemission 
mechanism (see Sec. V), these temperatures do not pro- 
duce enough photons to comply with experimental re- 
sults, Toegel et al. (2002) address this paradox and sug- 
gest that according to the Le Chatelier- Brown principle, 
the high densities inside the bubble favor the back reac- 
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tion H+ OH—H;0, in particular because of the finite 
excluded volume of the particles. The energy-consuming 
water dissociation is thus reduced, and higher tempera- 
tures in the bubble are possible (see also Sec. V). A 
result like this suggests that one cannot extrapolate the 
reaction rates at temperatures and pressures achieved in 
other laboratory experiments to the unusual regime of 
sonoluminescence particle densities. Thus it is advisable 
to revert to a first-principles statistical physics approach 
in order to deduce reaction rates under SBSL condi- 
tions, deriving the laws of mass action directly from the 
partition function (Toegel et al., 2002). 

Another similar, simple model with both heat and 
water-vapor exchange was developed by Storey and 
Szeri (2001). Here the authors even drop the assumption 
of boundary layers for heat and mass transport and rely 
on the ratios of the relevant time scales for bubble dy- 
namics, particle diffusion, and heat diffusion. This ap- 
proach still contains the same essential physics as the 
models of Yasui and Toegel etal. Storey and Szeri 
(2001) quantitatively tested the simple model against 
their full simulations (Storey and Szeri, 2000). In both 
approaches, they find the same trends (and similar val- 
ues) for the peak temperature in the bubble, the mass 
percentage of argon in the bubble, and the number of 
reaction products, which again lends credibility to the 
simple ODE approach, 


C. How accurate are the bubble temperatures? 


It was mentioned in the introductory remarks that an 
understanding of the bubble interior and SBSL light 
emission hinges on a good understanding of the bubble 
temperature. The many models developed so far, with 
different degrees of sophistication, result in predictions 
for the peak temperature (suitably averaged over the 
bubble) ranging from 6000 K to slightly above 20000 K, 
given the same (typical) driving parameters for an SBSL 
bubble. While this factor of uncertainty of about 3 is 
significant (and leads to widely different results for the 
ensuing light emission), it is reassuring that these models 
share a number of common traits: (i) They agree that 
the interior of the bubble heats up and becomes at least. 
as hot as that measured for MBSL bubbles; (ii) water 
vapor is a major temperature-limiting factor, forming a 
substantial part of the bubble contents at collapse; and 
(iii) the temperature, when averaged over the bubble, is 
not likely to rise much above 10* K, ruling out the much 
ioo high predictions of earlier attempts at modeling 
SBSL. We shall see in Sec. V that a thermal origin of 
SBSL is deemed very likely today, and that the tempera- 
ture range presented here does support light emission in 
the experimentally observed range of intensities. In- 
deed, any nonthermal theory of SBSL has to explain 
first why thermal emission, which explains the experi- 
mental results, is suppressed. 

The uncertainties in bubble temperature and the re- 
sulting predictions for light emission do not allow for a 
direct comparison between experiment and theory for a 
specific combination of parameters (e.g., a bubble with 
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FIG. 29. Minimal bubble radius Ryjg/Ry as a function of Ry 
and P,,. The frequency is f=26.5 kHz. 


R75 um, P,=13 bars, and f=20 kHz). On top of the 
modeling uncertainties, the sensitive dependence of light 
emission on parameters severely limits the reproducibil- 
ity of measurements. A comparison to experiment over 
a whole range of parameter values is much more prom- 
ising, focusing on the experimentally robust dependence 
of bubble dynamics and light emission on various exter- 
nal parameters. The following section takes this step 
from single-parameter combinations to an overview of 
the parameter space of SBSL—starting with the ques- 
tion of under what conditions sonoluminescence can be 
observed at all. 


IV. THE PARAMETER RANGE OF SINGLE-BUBBLE 
SONOLUMINESCENCE 


In the previous section we studied the dynamics of a 
bubble under the action of a driving pressure of ampli- 
tude P,, frequency f=w/2m, and with an ambient 
bubble radius Ro, all fixed to appropriate values for 
sonoluminescence to occur. The key question now is 
what happens if these and other experimental param- 
eters are changed—will SBSL still be observed, and if 
so, will it be brighter or dimmer? What are the physical 
processes that determine the limits of the parameter re- 
gime of sonoluminescence? Experiments have shown 
that, apart from the parameters implicit in the Rayleigh- 
Plesset equation, other quantities of crucial importance 
to SBSL are the concentration of gas c., dissolved in the 
liquid, the temperature of the liquid, the type of liquid, 
and the type of gas. The goal of this section is to under 
stand these parameter dependencies quantitatively. 
Various physical constraints limit the parameter range 
in which sonoluminescence can be observed: to emit 
light, the bubble must be forced strongly enough for a 
cavitation event to occur during each cycle of the drive; 
the bubble must not break into pieces, which roughly 
translates into the requirement that viscous processes 
and surface tension be strong enough to limit the growth 
of bubble shape instabilities. For the consistent, stable 
light emission of SBSL, the number of gas molecules 
inside the bubble, averaged over one cycle of oscillation, 
must neither increase nor decrease. This requirement is 
what sets the ambient radius Ro of the bubble; it in- 
volves a subtle interplay between diffusive processes ex- 
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changing gas between the bubble interior and the out- 
side liquid, and chemical reactions. And finally, it is 
necessary that the Bjerknes forces holding the bubble 
trapped in the flask (see Sec. ILG) be strong enough to 
ensure that the center of the bubble does not move ap- 
preciably. This section presents the current understand- 
ing of each of these effects and assesses the extent to 
which the theoretical predictions agree with experi- 
ments. 


A. The Blake threshold 


Regardless of the exact mechanism of sonolumines- 
cence, it is abundantly clear that the light results from 
energy focusing during a rapid bubble collapse. There- 
fore the bubble must be forced strongly enough to in- 
duce a cavitation event of sufficient violence—in es- 
sence, the Rayleigh collapse solution of Sec. ILD must 
be fully established. Whether this happens depends on 
both the ambient bubble size (mass of gas inside the 
bubble) and the forcing pressure. Figure 29 shows the 
minimal radius during a cycle of the drive as a function 
of forcing pressure P, and ambient radius Ro. There is 
an abrupt transition Ro(P,) where the onset of Rayleigh 
collapse occurs and the gas inside the bubble gets 
strongly compressed, leading to heating. Therefore 
sonoluminescence can only occur above this threshold. 

The functional form of this threshold curve can be 
deduced from dimensional considerations: other than 
the ambient radius Ro, the parameters in the Rayleigh- 
Plesset equation are P, Po, o; p, 7, and o. Since we are 
dealing with the transition from gentle oscillation to in- 
ertial collapse, we do not expect dissipative effects to be 
important and exclude the viscosity 7 from our consid- 
erations, From the remaining parameters, two indepen- 
dent length scales can be derived: the resonant bubble 
size Rr - \3Pop Tw 7 [cf. Eq. (22)] and the capillary 
length scale a/ P. As we have seen in previous sections, 
SBSL bubbles are driven far from resonance, so Rg” is 
too large to be important here. 

The relationship between the critical ambient size and 
the pressure must therefore have the form 





(43) 





with a dimensionless function f. A closer look at the 
Rayleigh-Plesset equation allows for the specialization 


(44) 





where C is a scalar constant, A calculation like this was 
first performed by Blake (1949) for bubbles under static 
pressure. He found C=2/3 for isothermal bubble move- 
ment (ct. Brennen, 1995). The threshold in Ro-P, space 
separating gently oscillating from violently collapsing 
bubbles is therefore known as the Blake threshold. 
More detailed studies for driven bubbles arrive at C 
=4V3/9~0.77 (Leighton, 1994; Hilgenfeldt, Brenner, 
et al., 1998). 
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Only bubbles larger than R$ in Eq. (44) can show 
SBSL, so that the Blake threshold criterion cuts down 
the available parameter space considerably. When shock 
waves inside the bubble were considered crucial for the 
light emission, Hilgenfeldt et al. (1996) suggested replac- 
ing this criterion with the threshold for strong shock for- 
mation, expected when the maximum bubble speed dur- 
ing collapse exceeds the speed of sound in the gas 
(although shocks can and will form at lower speeds as 
well). Because of the abruptness of the Blake threshold, 
supersonic velocities are reached already at parameter 
combinations very close to the threshold. Thus this 
shock threshold, like any other dynamical criterion for 
strongly collapsing bubbles, will effectively yield the 
same result as the Blake threshold calculations. 

Further constraints set an upper limit to the bubble 
size, in order to ensure the bubble’s stability. Let us first 
consider diffusive stability. 


B. Diffusive stability 


Since bubble dynamics and energy focusing during 
collapse are sensitively dependent on Ro, it is crucial for 
a stable SBSL bubble to maintain the same ambient ra- 
dius, i.e., not to exchange any net mass with its sur- 
roundings. The gas exchange between the bubble and 
the liquid is affected by diffusion of gas through the liq- 
uid and by advection of this gas with the fluid velocity. 

The typical model that is used for this process starts 
with the transport equation for the mass concentration 
c(r,t) (mass/volume) of gas around a spherical bubble: 





ac+ud,c=D marae), (45) 
where D is the gas diffusion coefficient in water. The 
velocity u(r,t)=R?R/r? of the fluid a distance r from 
the center of the bubble is given by Eq. (6), with R(t) 
entering through solution of the Rayleigh-Plesset equa- 
tion. The gas in the bubble is assumed to remain in equi- 
librium with that in the liquid at the boundary of the 
bubble wall; hence the gas concentration at the bubble 
wall is given by Henry's law, 


c(R,t) e copy (Rt) Py (46) 


Finally, the gas concentration far from the bubble is 
given by the ambient concentration c., 


e(,t)=en. an 


‘The mass loss/gain of the bubble is then proportional to 
the concentration gradient at the bubble wall, 


mh=47R?Dð,e\ri) + (48) 


This model for the gas exchange is accurate as long as 
the Rayleigh-Plesset solution is valid. As was empha- 
sized above, this may not be generally true very close to 
the point of bubble collapse. At collapse, many other 
effects could also play an important role and affect the 
gas transfer, including (a) sound waves in the liquid, (b) 
breakdown of Henry's law at the interface because of its 
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fast motion, (c) chemical reactions inside the collapsed 
bubble, and (d) phase transformation (boiling) of the 
liquid surrounding the bubble. 

First, we shall examine the consequences of mass dif- 
fusion alone and see what they imply about the sonolu- 
minescence experiments. The advection of the solute in 
Eq. (45) reflects the conservation of fluid volume: a 
spherical shell of fluid around the oscillating bubble is 
pushed in or out when the bubble contracts or expands. 
For this reason, it is useful to think about the solutions 
to the equation in Lagrangian coordinates tracking the 
volume changes according to (Plesset and Zwick, 1952) 

h(r,t)= $[P- R0]. (49) 
This transformation trick has often been used (e 
Eller, 1969 or Brennen, 1995). If one was to neglect dif- 
fusion, the concentration field outside the bubble would 
be just c=c(h). However, even in the limit where the 
mass diffusion coefficient is very small, it has important 
consequences (for example, an undriven bubble eventu- 
ally dissolves by diffusion). The Henry's law boundary 
condition at the bubble wall implies that the gas concen- 
tration at the bubble wall is time dependent and in gen- 
eral different from the concentration c, in the bulk. 
‘There is therefore a boundary layer near the bubble 
wall, where the gas concentration relaxes from the value 
dictated by Henry's law to c... The numerical profile is 
shown in the paper by Hilgenfeldt et al. (1996). 

The concentration gradient in the boundary layer dic- 
tates the mass flux into the bubble. There are two time 
scales over which mass transfer occurs: (i) during a 
single cycle of driving, gas is pushed into and out of the 
bubble when the bubble expands and contracts; (ii) over 
many cycles of driving, small net gains or losses of mass 
potentially accumulate to produce significant changes in 
the number of molecules in the bubble. 

With the coordinate (49), we obtain from Eq. (45) 


4,c=Da(r*a,c). (50) 


Focusing on the boundary layer region close to the 
bubble (r~R) and redefining the time 7= f R*dt yields 


d.c— Dye, 61) 


which is a pure diffusion equation. Following Fyrillas 
and Szeri (1994), the solution to this equation can be 
written c(Ji, 7) - c. cos (It, 7) + Cam (hy). Here, Cose is 
an oscillatory solution, transporting gas back and forth 
at the frequency of the bubble oscillations, but effecting 
negligible net gas transport. The smooth part of the pro- 
file c,,, however, determines the gas exchange over 
time scales much longer than the oscillation period. For 
the number of molecules in the bubble to maintain equi- 
librium over long periods of time, it is necessary that 
(c(t,7))=cx, where (-) denotes averaging over 7. Spe- 
cializing to the surface of the bubble, this implies (be- 
cause of Henry's law) 
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small forcing pressure amplitudes, P,=10atm to P, 
14 atm (top to bottom). From Hilgenfeldt et al. (1996). 





The index 4 indicates that the average is weighted with 
RS. The mass loss from the bubble over long periods of 
time can be calculated via 


(c-c- 


-AzR P^ 





63) 





where dis the boundary layer thickness, given by 


af f" dh 
Zn 


nj [35 R*(0)] ^t 


a 


, (54) 


where the time integral spans one oscillation period Tg. 
Because of the dominance of the maximum radius in 
the integration (see also below), 6 can be approximated 
as ~Ri/Rmax (cf. Hilgenfeldt et al., 1996). The time 
scale of mass diffusion is therefore approximately 
Ripg!RnaxDco~0.1 8, with the gas density p, under 
standard conditions. 

Fyrillas and Szeri (1994) applied the method of mul- 
tiple scales (see, e.g., Hinch, 1991) to systematically de- 
rive the above formulas and demonstrate their accuracy 
to leading order in D(oR;) !. As long as the diffusive 
time scale is much longer than the bubble oscillation. 
cycle, the approximation is very reliable and can there- 
fore be used instead of the more cumbersome direct so- 
lution of Eq. (45). 

Equation (53) can be used to study the stability of the 
equilibrium ambient radius Rj: taking Ro=Rj+e and 
expanding (c)(Ro) around Rj implies that 


D dics 
prr nM 








The equilibrium point is therefore stable as long as 
d(c)4/dRo>0 or, equivalently, d(p)4/dRo>0. 

Löfstedt et al. (1995) pointed out a useful approxi- 
mate version of Eq. (52): Since the time average is 
weighted by R^, and the bubble spends a large fraction 
of the cycle near the maximum radius, the equilibrium 
condition is approximately 
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FIG. 31. Bifurcation diagrams in the Ro-P, parameter space. 
The regimes with positive slope are stable. Gas concentrations 
are (right) ca/co=0.002, (middle) c./co7042, and (left) 
c, [cy 05. To the left of the curves the bubbles shrink and 
finally dissolve, to the right of them they grow by rectified 
diffusion. From Hilgenfeldt et al. (1996). 





(55) 





What do these results imply for sonoluminescence? Fig- 
ure 30 plots (p,)4 as a function of Ro for an argon 
bubble for various forcing pressure amplitudes P,. We 
first examine small pressure, P,~1.0 atm, and gas con- 
centrations of, say, c»/co~0.3. There is an unstable 
equilibrium at R§~6 um. Smaller bubbles shrink and 
larger bubbles grow. The mechanism of growth is called 
rectified diffusion: when the bubble is large, the gas con- 
centration in the bubble is small, and vice versa, so there 
is an influx or outflux of material. The net effect is not 
zero for a nonlinear oscillation because at maximum ra- 
dius (i) the area for gas exchange is larger, and (ii) the 
diffusive boundary layer is stretched out (Brennen, 
1995). Both effects favor net growth of the bubble. Note 
that an undriven bubble dissolves because its internal 
pressure exceeds that in the liquid. Sufficiently strong 
driving will overcome this tendency to shrink, and start 
rectified diffusion. 

‘Turning our attention to higher forcing pressures (Fig. 
30), we see that the average pressure is quite diminished, 
because Rmar is very large now [see Eq. (55)]. This 
means very low gas concentrations (saturation levels 
calco) are needed to achieve diffusive equilibrium 
(around 0.005 for P,=1.3 atm). Moreover, under these 
conditions the equilibria become stable, as demonstrated 
by the positive slopes. In this example, the bubble 
reaches an equilibrium size of Rj-5 um and will not 
gain or lose any further mass. 

Figure 31 is a diagram showing the equilibrium bubble 
states in the Ro-P, parameter space for three fixed con- 
centrations (c../¢y=0.5,0.02,0.002). Stable equilibria 
have positive slope @Ro/4P,|-,; negative slope repre- 
sents unstable equilibria. To the right of each line 
bubbles grow, and to the left they shrink. 
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Figure 31 shows that there are no diffusively stable 
sonoluminescence bubbles for large gas concentrations, 
where all equilibria are unstable.” For small concentra- 
tions c../co the situation is quite different. There are 
stable equilibria at large P, and small Rọ. Only in this 
region, and for very low gas concentration cco 
~0,001-0.02 (depending on P,,) is the bubble diffu- 
sively stable, and only then is stable SBSL possible. This 
is again easily appreciated using the approximation (55) 
of Löfstedt ef al. (1995): when the bubble is small and 
strongly forced, the collapse ratio is R max/Ro~=10. Hence 
Eq. (55) implies that c, /cg~ 107° for the bubble to be in 
equilibrium, as observed in Fig. 31. 

Löfstedt er al. (1995) realized that this requirement of 
diffusive equilibrium for strongly forced bubbles can 
only be fulfilled at very small gas concentrations, in 
agreement with experiments using argon or other inert 
gas bubbles. However, the finding blatantly disagrees 
with the results obtained for air bubbles, where stability 
is achieved at roughly 100 times larger gas concentra- 
tions. Recall that Gaitan (1990) needed to degas to only 
about 40% of saturation. If he had had to go to 100times 
lower concentrations, he might never have discovered 
SBSL. Lófstedt et al. (1995) postulated an “anomalous 
mass flow mechanism” to resolve this discrepancy. In- 
deed, in order to keep a constant Ro in a liquid with 
such a high gas concentration, an air bubble would have 
to eject mass far above the usual diffusive transport rate 
in order to balance rectified diffusion. 


C. Sonoluminescing bubbles rectify inert gases 


1. The mechanism 


‘The stabilizing mass ejection mechanism of sonolumi- 
nescing bubbles is now believed to be the consequence 
of chemical reactions that occur when the bubble is com- 
pressed. The maximum temperature of the bubble 
(larger than 10000 K; see Sec. III) is large enough to 
destroy both molecular nitrogen and oxygen, so these 
molecules will be largely dissociated in the collapsed 
bubble. Moreover, as discussed in Sec. III, the high ex- 
pansion ratio of the bubble results in a substantial 
amount of water vapor. Chemical reactions between wa- 
ter vapor and dissociated nitrogen and oxygen are well 
known to atmospheric chemists dealing with acid rain: 
the reaction products are NO, NH, and ultimately 
HNO, and NH;. All of these substances (except NO) 
are very soluble in water. The idea of the argon rectifi- 
cation hypothesis (Lohse et al., 1997) is that the reaction 
products are absorbed completely into the water, deplet- 





5Strictly speaking, stable equilibria do appear due to the 
“wiggles” in the curve of equilibria that result from resonance 
effects. Brenner, Lohse, et al. (1996) speculated that the 
wiggles could describe multiple stable equilibria. However, 
when introducing further damping mechanisms into the 
Rayleigh-Plesset dynamics as done by Hao and Prosperetti 
(199%), the wiggles weaken and it is presently not clear 
whether there are multiple stable equilibria, 
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ing the bubble of nitrogen and oxygen and thus estab- 
lishing an efficient mass-loss mechanism. A sonolumi- 
nescing air bubble thus rectifies argon, the only 
substance inside the bubble that does not dissociate, 
which is contained in air with a concentration of 1%. 
This argument immediately explains the discrepancies 
between the diffusive equilibria of air and argon: only 
the inert gas in the bubble is in diffusive equilibrium 
with the bulk liquid. Therefore the relevant parameter 
for the stability of sonoluminescing bubbles in equations 
like Eq. (52) is the partial pressure (or concentration) of 
the inert gas dissolved in the liquid, and not the partial 
pressure of air. Defining the argon fraction q as. 
& 


c 


co co” 

the stability criteria for air and argon should differ by a 
factor of q~0.01, the fraction of argon in air. Indeed, 
experiments show stable sonoluminescence in air 
bubbles for c./co~0.2-0.4 (Gaitan etal, 1992) and 
stable sonoluminescence in argon for c,/c 
—0.002-0.004 (Barber etal., 1995). The hundredfold 
difference is quantitatively consistent with the argon rec- 
tification theory. The earliest measurements confirming 
consequences of argon rectification were by Holt and 
Gaitan (1996), described in detail below. 

‘As pointed out by Lohse and Hilgenfeldt (1997), an- 
other hint for the chemical activity inside the bubble is 
the isotope scrambling found by Hiller and Putterman 
(1995), who had analyzed sonoluminescence in H; and 
D; gas bubbles, in both normal and heavy water. The 
four respective spectra are grouped according to the liq- 
uid, not according to the dissolved gas. This suggests that 
water vapor invades the bubble and undergoes chemical 
reactions, ultimately leaving Ha for normal water and 
D, for heavy water, independent of the type of hydrogen. 
bubble with which one started. Hiller and Putterman's 
(1995) hydrogen bubbles were all diffusively unstable. 

The most direct verification of argon rectification was 
accomplished in a clever experiment by Matula and 
Crum (1998). They were able to precisely measure the 
transition time over which the light intensity changed 
when the driving conditions were changed. They ob- 
served (see Fig. 32) that an air bubble that has not pre- 
viously emitted light reaches maximum sonolumines- 
cence intensity after a transition time of about 10° 
driving cycles. By contrast, if a bubble initially emits 
light, when the forcing pressure is dropped below the 
light-emitting threshold and then immediately raised 
back above it, sonoluminescence turns on immediately. 
As experimental controls, Matula and Crum used pure 
argon and pure nitrogen bubbles. The transition time for 
the nitrogen bubble matched the air bubble that initially 
did not emit light, and the argon bubble matched the air 
bubble which had sonoluminesced previously. This 
therefore provides direct evidence for the hysteresis that 
would be expected with argon rectification. 

Other predictions from the rectification hypothesis 
have not been experimentally verified to our knowledge, 
e.g., the suggestion for making stable sonoluminescence 
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without degassing the liquid, preparing a percentage of 
argon in nitrogen so that cA'/cg—0.003, within the diffu- 
sive stability window. The ultimate proof of the argon 
rectification hypothesis would be the detection of the 
chemical species leaving the bubble. Lohse ef al. (1997) 
suggested detecting a pH change due to the acidic reac- 
tion products. However, pH meters are probably not 
sensitive enough to detect an effect even if the experi- 
ment is run for a long time. Therefore Lepoint et al. 
(1999) exchanged the water in the experiment for 
Weissler's reagent. The oscillating bubble produced per- 
oxide and chlorine radicals which oxidized iodide to io- 
dine, giving the distinct blue color of the iodine-starch 
reaction, Thus chemical reactions in and around a single 
bubble were conclusively demonstrated. The bubble 
triggered the reaction even at driving pressures below 
the SBSL threshold, as the temperature there can al- 
ready be sufficiently large. A thread of blue color was 
observed to emerge from the bubble, usually in either an 
upward or downward direction, 


2. Bubble equilibria with chemical reactions 


Chemical reactions modify the classical diffusive equi- 
libria described above. The arguments outlined above 
can only make a quantitative prediction for the extreme 
cases of no chemistry or complete argon rectification. To 
describe the transition between these two stages (for ex- 
ample, as the driving pressure acting upon an air bubble 
is increased), the chemical reaction rates have to be 
modeled explicitly. As described in Sec. III, this model- 
ing can be done to varying degrees of detail, trading 
accuracy for computational speed. Full simulations are 
still too expensive computationally to allow a mapping 
of the parameter space of SBSL. Therefore simple reac- 
tion models like those of Yasui (1997a), Toegel et al. 
(Toegel, Gompf, et al., 2000; Toegel et al., 2002), or Sto- 
rey and Szeri (2001) are often employed. 

In such a model, the changes in species numbers are 
given by diffusive transport [see Eq. (39) for water va- 
por], and by chemical reactions such as Eq. (41). Reac- 
tion rates under SBSL conditions must be either in- 
ferred from general principles (cf. Toegel et al., 2002) or 
extrapolated from lower-temperature reaction data (cf. 
Kamath et al., 1993; Bernstein et al., 1996; Yasui, 1997b). 

As an example, the change in the numbers Na, and 
Ny, of argon and nitrogen particles in a mixture of these 


two gases can be written as 














ARD Ak 
ac AR Dae o (56) 
Mar 
ATR Dy, ð," (om 
Nyy S TEM -ANS exp| = Tc] 
(57) 


where the first term in Eq. (57) represents diffusion 
and the second dissociation reactions, depending on 
the bubble temperature T(t). Dar, Dy, Har and ux, 
are the respective diffusion coefficients and molar 
masses. For simplicity, Lohse eral (1997) assumed 
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FIG. 32. (+) The maximum radius and (X) sonolumines- 
cence intensity of an air bubble plotted over approximately 
240 consecutive cycles. The "radius" is actually a signal level, 
proportional to the square of maximum radius (assuming a 
geometrical optics limit for Mie scattering from a sphere). (a) 
The bubble is initially below the light-emission threshold (the. 
photomultiplier tube signal level corresponds to noise). At ap- 
proximately the 120th acoustic cycle, the drive pressure is rap- 
idly increased to a value above the SBSL threshold. (b) For 
this case, the bubble is initially in a stable sonoluminescing 
state, The drive pressure is then rapidly reduced to a value 
below the light-emission threshold, and then quickly increased 
again after a time period of approximately 90 ms. In both (a) 
and (b), the drive amplitudes at the lower and upper values are 
the same. From Matula and Crum (1998). 


that the reaction rates follow a modified Arrhenius 
law, using empirical parameters from Bernstein et al. 
(1996) appropriate for nitrogen dissociation: A~6 
x10%(Ty/TY LN v) (to/ R)? cm*/(mols) gives the 
time scale of the reaction; T*~113000 K is the activa- 
tion temperature, T, is ambient temperature, and po the 
equilibrium gas density This reaction law is very 
crude—it neglects backward reactions as well as the ki- 
netics of the expulsion of reaction products; however, it 
is sufficient for a demonstrative calculation. 

The concentration fields c^(r,t),c™ (r,t) in Eqs. 
(56) and (57) separately obey Henry's law at the bubble 
wall, using the saturation concentrations c," c? specific 
for the gases (cã"=3c}?). For diffusive stability, both 
gas species must fulfill the equilibrium condition (52) 
separately, with concentrations far from the bubble 

eM ont by the percentage £j of argon in the gas 
given by the pi ge Ë ge gi 
dissolved in the liquid. 


Requiring the temporal averages of Nar and Ny, to 
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vanish yields the number of gas particles in the bubble at 
equilibrium and thus the percentage £, of argon inside 
the bubble, which is larger than £j because of dissocia- 
tion reactions removing N;. Of course, the last term in 
Eq. (57) only contributes when T(r) is large. The chemi- 
cal reaction rate therefore depends on the detailed space 
and time dependence of the temperature in the bubble. 
‘The calculations below use a simple spatially uniform 
temperature model inspired by Prosperetti (1977c) as 
discussed in Sec. IILB.1. 

The resulting equilibrium radii Rj in the Ry-P, plane 
for air (£j-0.01) at p. /P,—020 are shown in Fig. 10. 
For small forcing, the temperatures are not high enough 
to initiate chemical reactions, so that the curve A of 
(unstable) equilibria corresponds to those described 
above for diffusion alone. These bubbles either shrink or 
grow by rectified diffusion, and for them £,-£j to a 
good approximation. At high forcing (curve C), the re- 
actions burn off all the Nz, so that the bubble contains 
pure argon (£,~1); this equilibrium corresponds to the 
(stable) equilibrium at the argon partial pressure 
p? IPq— 001p. P5 0.002. 

Note that curves A and C belong to the same experi- 
mental system now, and that bubbles of low P, and large 
Ry grow by rectified diffusion, while those with high P, 
and large Ro shrink due to the mass loss through chemi- 
cal reactions (arrows in Fig. 10). There necessarily must 
be a region of intermediate forcing pressures where both 
processes cancel, and an additional equilibrium exists. 
‘These equilibria, which prove to be stable, are shown as 
curve B in Fig. 10. The onset of appreciable nitrogen 
dissociation (T~9000 K) is depicted as a thin line in 
the figure and is quite close to line B. Note that this 
temperature is much smaller than T*. The argon frac- 
tion £, along curve B varies with P,, but stays consid- 
erably smaller than 1. 

This picture predicts the following sequence of events 
as the forcing pressure is increased. At low P,, the 
bubble starts near the unstable equilibrium curve A, 
where the bubble is growing through rectified diffusion 
and eventually undergoes shape instabilities (see Sec. 
IVD below). If the forcing pressure is turned up fairly 
rapidly, the ambient radius will grow in this regime. At 
some point, the forcing P, will be large enough so that 
the bubble is trapped by the stable equilibrium B. When 
P, is now further increased, the ambient radius shrinks 
along the curve B. Upon increasing the driving even fur- 
ther, fluctuations throw the bubble onto the stable curve 
C, so the ambient radius grows again. This behavior had 
been observed previously by Barber er al. (1994); see 
Fig. 7. A direct measurement was made by Gaitan and 
Holt (1999), who measured both the maximum bubble 
radius and the ambient radius as a function of forcing 
pressure (see Fig. 33). The breakdown in Ro and Rmax 
caused by the onset of chemical reactions inside the 
bubble is clearly visible. 

Holt and Gaitan (1996) showed that at p./Po=02 
there is a relatively large forcing pressure regime P, 
712-13 atm where bubbles dissolve (cf. Fig. 10). Such 
dissolution islands do not exist with pure diffusion (as 
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FIG. 33. Experimental maximal bubble radius and phase dia- 
gram for SBSL: (a) the ratio between maximum radius Rmax 
and ambient radius Ro; (b) Rmax as a function of the driving 
pressure P,; (c) data points in the P, vs Ry phase space; O, 
stable glowing bubbles; ll, stable nonglowing bubbles; +, un- 
stable, nonglowing bubbles; dashed line, the experimentally 
found shape instability; solid lines, theoretical diffusive equi- 
libria for various argon concentrations. To compare this figure 
with theoretical phase diagrams such as Fig. 10 or other experi- 
mental results such as Fig. 12, one should exchange abscissa 
and ordinate. From Gaitan and Holt (1999). 


shown above; see Fig. 31). Gaitan and Holt (1999) dem- 
onstrated that the stable equilibria of sonoluminescing 
bubbles were in excellent agreement with the shape pre- 
dicted when assuming a much lower gas concentration. 
These observations agree with the theoretical results 
taking both diffusion and molecular dissociation into ac- 
count. 

In their comprehensive experimental study of bubble 
stability diagrams, Ketterling and Apfel (1998, 2000a) 
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FIG. 34. Experimental data in the phase diagram in the P,-Ro 
plane: +, with an air concentration of c#/eq=0.2; O, 0.15; +, 
01; M, 0.02. The data lie on the stable diffusive equilibria 
curves (solid lines) for argon concentrations of c“"/cy=0.002, 
0.0015, 0.001, and 0.0002, confirming the argon rectification 
theory. Experimental data are taken for as large a P, as pos- 
sible, to probe the shape instability threshold. From Simon 
et al. (2001). 





demonstrated that all stable and light-emitting bubbles 
indeed lie on a positively sloped equilibrium curve in the 
Ro-P, plane, corresponding to the partial concentration 
of argon in the mixture. This is true regardless of how 
much nitrogen is present. A representative phase dia- 
gram is shown in Fig. 12 for air in water at c. /co702. 
‘The solid line represents the (unstable) diffusive equilib- 
rium for the mixture, and the dashed line is the (stable) 
diffusive equilibrium for c../¢o=0.002 (the argon in the 
mixture). As the forcing pressure is increased, the 
bubble follows a stable equilibrium line above the un- 
stable diffusive equilibrium (corresponding to line B in 
Fig. 10) and then transitions to the lower-concentration 
stable equilibrium (where it starts emitting light.) A 
pure argon bubble at c../co=0,0026 follows its stable 
diffusive equilibrium (see Fig. 5 of Ketterling and Apfel, 
1998). No other stable configurations are possible. On 
the other hand, Fig. 4 of Ketterling and Apfel (1998) 
shows that a pure nitrogen bubble at c /cọ=0.1 follows 
a stable equilibrium much like curve B of the air mix- 
ture, but cannot reach a stable sonoluminescing state at 
higher P,. Taken together, these two types of behavior 
synthesize the results in Fig. 12, strongly implying that 
the air bubble is composed mainly of nitrogen along 
curve B, and mainly argon along curve C, in agreement 
with the predictions of argon rectification. 

Very recently, Simon et al. (2001) confirmed that light- 
emitting bubbles follow the stable diffusive equilibrium 
curves based on argon concentration alone (see Fig. 34). 
They employ a new experimental technique for measur- 
ing the parameters of the bubble dynamics (P, and Ro) 
based on the timing of the light flash in the acoustic 
period. They also find that the attainable conditions in- 
side a sonoluminescing bubble are more extreme at a 
low partial air pressure of 15 mm Hg as compared to 150 
mmHg, because then the bubbles are closer to the 
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Blake threshold (where the most extreme conditions are 
achieved) and can be driven with larger P,. This is due 
to the bubble shape instabilities treated in the next sub- 
section, 





D. Shape stability 


The theoretical diffusive equilibrium curves stretch to 
far larger ambient radii Ro than those observed for 
SBSL bubbles. There must be another requirement lim- 
iting the ambient size (or the total mass content) of the 
bubble. This limit is set by the onset of instabilities in 
the shape of the oscillating bubble. The analysis of shape 
stability is a classical problem in bubble dynamics, pio- 
neered by Plesset (1949), Birkhoff (1954), Plesset and 
Mitchell (1956), Strube (1971), and Prosperetti (19774). 
In this section we present the application of these ideas 
to single-bubble sonoluminescence. 


1. Dynamical equations 


To analyze the linear stability of the radial solution 
R(t), consider a small distortion of the spherical inter- 
face R(t), 


T=R(t)+0,()¥9(6,6), (58) 


where Y, is a spherical harmonic of degree n. The goal 
is to determine the dynamics a,(1) for each mode. The 
derivation of Plesset (1954) follows along the same lines 
as the derivation of the Rayleigh-Plesset equation, 
which it recovers to zeroth order in a, . A potential flow 
outside the bubble is constructed to satisfy the boundary 
condition that the velocity at the bubble wall be Ñ 
+4,Y,,. This potential is then used in Bernoulli's law to 
determine the pressure in the liquid at the bubble wall. 
If viscous effects are neglected, applying the pressure 
jump condition across the interface yields a dynamical 
equation for the distortion amplitude a,(1), 


d, By(t)n—Ag(t)aq=0, (59) 
where ,=(n—1)(n-+1)(n+2) and 
B,()=3RIR, (60) 
ROB, 
Aun-|o-ob- 85] (61) 


The stability of the spherical bubble then depends on 
whether solutions to Eq. (59) grow or shrink with time. 
Itis already apparent here that Eq. (59) has the form of 
a parametrically driven oscillator equation (Hill equa- 
tion), with the radial dynamics R(t) governing the peri- 
odic driving. 

‘A more accurate stability analysis requires taking ac- 
count of viscosity and other dissipative processes. Vis- 
cosity, treated by Prosperetti (197d), poses difficulties 
because viscous stresses produce vorticity in the neigh- 
borhood of the bubble wall, which spreads convectively 
through the fluid. Once created, the vorticity acts back 
on the dynamics of a, (r). This interaction is nonlocal in 
time, and so the problem requires solving integrodiffer- 
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ential equations for the vorticity in the liquid, coupled 
with the shape oscillations. Details can be found in the 
literature (Prosperetti, 197d; Hilgenfeldt et al., 1996; 
Hao and Prosperetti, 1999b). Here we simply summarize 
the results of a simple “boundary layer” approximation, 
which assumes that the vorticity is localized in a thin 
region around the bubble. It was again Prosperetti 
(19775) who first realized the usefulness of this approxi- 
mation. If 3 is the boundary layer thickness, the prefac- 
tors of Eq. (59) are modified to 








2v 1 
p| Benni 1)02) z| 
142g 
(62) 
3R 2v n(n*2)* 
By()= Rz At —F |. (63) 
1425 


with the kinematic viscosity v of the liquid. The viscous 
contribution to A,(¢) is not important since the ratio 
between the third and the second terms of the right- 
hand side in Eq. (62) is typically voRoa/a-10-2. How- 
ever, in the second term of the right-hand side of Eq. 
(63) it introduces a damping rate which causes exponen- 
tial damping of shape modulations. The amount of 
damping strongly depends on both the boundary layer 
thickness and on n. Brenner et al. (1995) and Hilgen- 
feldt et al. (1996) choose 6 to be the minimum of the 
oscillatory boundary layer thickness 7w and the wave- 
length of the shape oscillation R/(2n). 

‘The Hill equation (59) is driven by the strongly non- 
linear Rayleigh-Plesset dynamics R(t). Therefore, in 
contrast to the monofrequent driving of the prototypical 
Mathieu equation, instabilities in a, can be excited on 
the many different time scales of the bubble oscillation 
discussed in Sec, II. In particular, three types of instabili- 
ties can be distinguished: the parametric instability (over 
time scales of the oscillation period), the afterbounce 
instability (over time scales of the bubble afterbounces), 
and the Rayleigh-Taylor instability (over time scales of 
the Rayleigh collapse). 





2. Parametric instability 


The parametric shape instability acts over the rela- 
tively long time scale T4=27/w (period of the driving). 
If the nonspherical perturbations of bubble shape show 
net growth over one oscillation period, they will over- 
whelm the bubble after many periods. This argument 
neglects possible (nonlinear) saturation effects not con- 
tained in the linear approximation (59), which could in- 
hibit further growth of the perturbations. 

In the relevant parameter regime for the parametric 
instability, R(r) and thus also A,(r) and B,(r) are 
strictly periodic in time with period 7.;. Thus the stabil- 
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ity of the Hill equation (59) can be rigorously analyzed 
(Nayfeh and Mook, 1979). Instability occurs whenever 
the magnitude of the maximum eigenvalue of the Flo- 
quet transition matrix F, (T4) of Eq. (59) is larger than 
1. The Floquet transition matrix is defined as the propa- 
gator of the perturbation vector over one period, 


(get Ta) a,(0) 


(ar)= FAT aro 


By numerically computing the eigenvalues of the Flo- 
quet transition matrix, one can map out the phase dia- 
gram of parametric stability, i.e., identify parametrically 
stable and unstable regions. 

In the sonoluminescence parameter range of P, 
~1,2-1.5 atm and a typical frequency f=26.5 kHz, cal- 
culations with the boundary layer approximation suggest 
that parametric instability sets in for ambient radii in 
excess of R}!~4—5 um, with only a weak dependence 
on P,. Refined boundary layer models like those of 
Prosperetti and Hao (1999) or Augsdórfer et al. (2000) 
take into account higher-order terms in R/c, heat losses, 
or the varying gas density in the bubble upon collapse. 
These models find upper stability bounds for Ry about 
half a micron larger, because the additional effects result 
in less violent oscillation dynamics and smaller values of 
the (destabilizing) bubble acceleration. 

For nonsonoluminescing bubbles, stability diagrams 
of a similar type were first measured by Eller and Crum 
(1970) and later by Horsburgh (1990), These experimen- 
tal and theoretical studies examine larger bubble sizes 
with smaller forcing pressures for which sonolumines- 
cence cannot occur. Applying the above shape stability 
analysis in the regime P,~0.5-1 atm (Brenner, Hilgen- 
feldt, and Lohse, 1998; Brenner et aL, 1999; Hao and 
Prosperetti, 1999b; Augsdórfer et al., 2000) gives similar 
thresholds to those found in experiment by Eller and 
Crum (1970), Horsburgh (1990), and Gaitan and Holt 
(1998). For P,<0.9 atm, the n=3 mode takes over as 
the most unstable surface mode from the n=2 mode 
(Brenner, Hilgenfeldt, and Lohse, 1998; Brenner et aL, 
1999; Augsdörfer et al., 2000), just as was found in the 
experiments by Gaitan and Holt (1998). Overall, larger 
bubbles can be stabilized at small P, , and the threshold 
in Ro becomes strongly dependent on the driving pres- 
sure in this regime. 

Holt and Gaitan (1996; Gaitan and Holt, 1999) also 
measured shape instabilities in the sonoluminescence re- 
gime and close to it. Their experimental results were 
compared with theory in the work of Brenner, Dupont, 
et al. (1998) and Hao and Prosperetti (1999b), with good 
agreement in evidence [see Fig. 35(a)]. Similar results 
are found for the higher n modes (see Hao and Prosper- 
etti, 1999b). 

The boundary layer approximation has been criticized 
by Putterman and Roberts (1998) as underestimating 
the amount of dissipation. These authors claim that 
there is insufficient evidence for the role of shape insta- 
bilities in limiting sonoluminescence, although Brenner, 
Dupont, et al. (1998) found that the approximation gives 






(64) 
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FIG. 35. Comparison of boundary layer approximation and 
full integrodifferential equation: (a) The dark area is the cal- 
culated stability region for the n=2 mode for an air bubble in 
water at f=20.6 kHz. The open circles are the data of Holt 
and Gaitan (1996; Gaitan and Holt, 1999). The calculation em- 
ployed the boundary layer approximation. (b) Same calcula- 
tion as above, but based on the full integrodifferential equa- 
tions, rather than the boundary layer approximation: Hardly 
any difference as compared to (a) can be observed. From Hao 
and Prosperetti (1999). 





results in good agreement with experiments. It is there- 
fore important to ask how the shape stability results 
change when the boundary layer approximation is nor 
made, and the full integrodifferential equations are 
solved. Such a comparison was carried out by Wu and 
Roberts (1998) and Hao and Prosperetti (199b). A 
typical result for forcing pressure around P,~1 atm is 
shown in Figs. 35(a) and 35(b), taken from Hao and 
Prosperetti (1999b): Hardly any difference between the 
exact result and the boundary layer approximation can 
be seen. Wu and Roberts (1998) found similar agree- 
ment even for larger P, in the regime of SBSL. These 
results indicate that the boundary layer approximation 
of Prosperetti (1997c), Brenner et al. (1995), and Hilgen- 
feldt et al. (1996) is appropriate for sonoluminescence 
experiments. 
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FIG. 36. Shape distortion of a bubble. (Upper part): Time 
development of the bubble radius R(r) (lower part) and dis- 
tortion amplitude a,(t) for a Ro=4.4 um bubble driven at 
P,=1.1 atm. Note the transition from Rayleigh-Taylor (time 
scale ns) to afterbounce perturbations (time scale us) during 
the afterbounce part of the bubble dynamics. It can also be 
seen that the dynamics of the distortion a(t) has half the 
frequency of the forcing radial dynamics R(¢), as is typical for 
an instability of the Mathieu type. From Hilgenfeldt et al. 
(1996). 





3. Afterbounce instability 


During the afterbounces, the bubble oscillates close to 
its resonance frequency (see Sec. II) on a time scale 7 
=1fo~ NpRU3Ps—-03 us. It turns out that the charac- 
teristic period of shape oscillations about the spherical 
bubble is very close to this resonant time scale, namely, 
VeRV(7B,)~1 us! VB, (for the n=2 mode, VB,~3). 
This coincidence of time scales is the root cause of the 
parametric instability (which exhibits maximal growth 
when the time scale of the forcing is of the order of the 











meli : 


FIG. 37. Mie scattering data of the afterbounce dynamics of a 
bubble driven below the luminescence threshold. Parametric 
instabilities can be inferred from the occurrence of large scat- 
tering spikes (near 9 and 10.5 us in this figure), a consequence. 
of the strong shape distortions of the bubble. Direct imaging of. 
the bubble (Matula, 1999) shows nonspherical bubble shapes. 
From Matula (1999). 
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time scale of the natural oscillation frequency). Under 
the right circumstances this instability can be so violent 
that the bubble is destroyed during the afterbounces of a 
single cycle. Examples for which the bubble “survived” 
considerable distortions during the afterbounces are 
shown in Fig. 36 (theory) and Fig. 37 (experiment), the 
latter taken from Matula (1999). The distortion can 
grow so much that the bubble breaks apart during the 
afterbounce period. The growth of instabilities during 
the afterbounce phase has been directly observed by 
Gaitan and Holt (1999). 

The afterbounce instabilities must be triggered by 
noise. A good way to analyze this dependence is to 
model the thermal noise through coupling a Langevin- 
type force to the dynamical Eq. (59) for the shape dis- 
tortion (Augsdórfer et al., 2000), with a magnitude ad- 
justed to satisfy the fluctuation-dissipation theorem. 


4. Rayleigh-Taylor instability 


The Rayleigh-Taylor shape instability occurs when a 
lighter fluid is accelerated into a heavier fluid (the clas- 
sical example is the interface between two layers of liq- 
uid, the lower one being lighter and with buoyancy as 
the accelerating force). For sonoluminescing bubbles, 
this shape instability acts over the extremely short time 
scales of the final stages of Rayleigh cavitation collapse. 
Here the bubble interface decelerates in preparation for 
the reexpansion, leading to an extremely large relative 
acceleration of the gas with respect to the water in ex- 
cess of 10? g. This deceleration occurs for only a short 
time (nanoseconds); it is roughly the time a sound wave 
of speed c, needs to cross a fully collapsed SBSL bubble 
of radius ~h. For the Rayleigh-Taylor instability to be 
effective, it must destroy the bubble during this time 
period. The competition between large magnitude and 
short duration of the accelerating force determines the 
stability threshold, A shock-front-driven variant of the 
Rayleigh-Taylor instability is the Richtmeyer-Meshkov 
instability, whose occurrence does not seem likely, since 
no evidence for shock-wave passage through the bubble 
wall was found (Wang et al., 1999). 

It should be emphasized that, in contrast to the para- 
metric and afterbounce instabilities, the Rayleigh-Taylor 
instability acts on such a short time scale that the bubble 
dynamical approximations cannot be expected to be 
quantitatively correct. A proper calculation requires a 
full simulation of the gas-liquid interface dynamics down 
to the latest stages of the Rayleigh collapse 

‘The thresholds for the Rayleigh-Taylor shape instabil- 
ity do not depend on whether the boundary layer ap- 
proximation is used or not, but instead on the chosen 
dynamics for the bubble radius R(t). Hilgenfeldt et al. 
(1996) used the modified Rayleigh-Plesset dynamics 
(20), ie., without the R/c corrections of the kinetic term 
[ct Eq. (18)}, and without thermal damping. Those 
terms were included by Prosperetti and Hao (1999), who 
found the Rayleigh-Taylor instability greatly suppressed 
due to the smaller bubble-wall acceleration A resulting 
from higher-order equations (see Prosperetti and Hao, 
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FIG. 38. Parameter space restrictions for sonoluminescing ar- 
gon bubbles: The M=1 curve (long-dashed) characterizes the 
onset of strong bubble collapse and heating. The bubble grows 
thanks to rectified diffusion to the right of the diffusive stabil- 
ity curves (heavy lines, shown for c» /¢y=0.5, 0.02, and 0,002, 
left to right). The thin solid line marks the onset of the para- 
metric instability and the short-dashed line combines the 
threshold of Rayleigh-Taylor instability and afterbounce insta- 
bilities. These lines are calculated within the simplified theory 
of Hilgenfeldt et al. (1996), which slightly underestimates the 
shape stability, as discussed in the text. 





1999). The uncertainties in modeling the bubble 
dynamics—and in particular in the second derivative 
Rare substantial enough to make a quantitative de- 
scription of the Rayleigh-Taylor instability a difficult 
task. Another factor that changes the location of the 
Rayleigh-Taylor instability line in the phase diagrams is 
the diminishing density contrast between the liquid and 
the extremely compressed gas at collapse, Taking this 
effect into account in Eq. (62), Augsdórfer et al. (2000) 
and Yuan et al. (2001) found further suppression of the 
Rayleigh-Taylor instability. 

In conclusion, though experimental results and theory 
give many hints that the Rayleigh-Taylor instability sets 
the upper threshold of the SBSL regime towards large 
forcing pressures P,, the matter is not yet fully settled. 
As discussed in Sec. II, the Bjerknes force instability 
may also play a role. The parametric and afterbounce 
instabilities, on the other hand, set well-established lim- 
its for the parameter space of SBSL towards large Ro. 
Figure 36 shows that instabilities such as Rayleigh- 
Taylor and afterbounce can occur simultaneously in the 
same bubble, with the perturbation a(r) growing over 
both nanosecond and microsecond time scales. 





5. Parameter dependence of the shape instabilities 


All calculations up to this point have been described 
for the material parameters of water at 20°C and driv- 
ing frequencies around f=20 kHz. The shape stability 
thresholds strongly depend on changes in the radial dy- 
namics brought about, for example, by different liquid 
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FIG. 39. Theoretical result from Hilgenfeldt et al. (1996) on 
the drift of the phase of light emission ó,(7) for three different 
argon concentrations, ¢./¢9=0.003 95, ¢x/cy=0.0658, and 
Ca /¢o=0.26, corresponding to gas pressures overhead of px 
mm Hg, 50 mm Hg, and 200 mm Hg, respectively. Accord- 
ing to the theory of Hilgenfeldt et al. (1996), the drift in the. 
phase of light emission is a result of bubble growth through 
rectified diffusion, which is followed by a pinch-off of a mi- 
crobubble when the bubble is running into the shape instabil- 
ity. The figure resembles the corresponding experimental re- 
sult, Fig. 4 of Barber et al. (1995). Note that for air bubbles 
stable sonoluminescence (corresponding to a constant phase of 
the light pulse) is achieved for much higher gas concentration, 
ct#/co=0.2, corresponding to 150 mm Hg [see also Fig. 4 of 
Barber et al. (1995)]. 




















viscosities or driving frequencies. The dependencies may 
act differently for different stabilities: if f is decreased, 
the parametric instability is suppressed as the stabilizing 
influence of viscosity can act for a longer time to sup- 
press perturbations. The more violent collapses of low- 
frequency-driven bubbles, however, favor an earlier on- 
set of the Rayleigh-Taylor instability. Nevertheless, for 
not too high forcing pressures, larger bubbles can be 
stabilized to show sonoluminescence, possibly emitting 
brighter light pulses (“upscaled” SBSL). However, as 
described in Sec. IILA.3, water vapor becomes increas- 
ingly important at low frequencies and counteracts the 
upscaling (Yasui, 1997b; Moss et al., 1999; Storey and 
Szeri, 2000; Toegel, Gompf, et al., 2000). 

The water temperature also has an effect on the phase 
space of SBSL. As pointed out by Hilgenfeldt, Lohse, 
and Moss (1998) and Vuong ef al. (1998), the increased 
viscosity at lower water temperature means that bubbles 
can be stabilized by cooling the fluid, allowing for much 
stronger acoustical driving and thus more light, as in- 
deed experimentally observed by Barber et al. (1994). 


E. Interplay of diffusive equilibria and shape instabilities 


The conditions for diffusive equilibrium and shape 
stability must be fulfilled simultaneously for stable SBSL 
(Hilgenfeldt er al., 1996). Outside this parameter regime, 
bubbles do not necessarily perish, but can undergo dy- 
namical processes like rectified diffusion that can allow 
for unsteady sonoluminescence at a weaker level (“un- 
stable sonoluminescence”). 

When at low forcing pressure P,~1.1 atm the gas 
concentration is sufficiently large (e.g, c«/co 5096), 
rectified diffusion can overcome the tendency for disso- 
lution, and growing bubbles are possible if R lies above 
the unstable equilibrium line. Rectification then contin- 
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ues until shape instabilities limit the growth (see Fig. 
38). When Ry reaches the boundary for shape instability, 
a microbubble pinches off, decreasing Rọ. If the remain- 
ing bubble is still large enough, i.e., above the unstable 
equilibrium line in Fig. 38, the process will repeat. For 
these low forcing pressures, the allowed size of the 
bubble after the pinch-off is very restricted. If the 
pinched-off microbubble is too large, the remaining 
bubble dissolves. 

For relatively large forcing pressure (P,~1.3 atm) 
the situation is quite different. For low enough gas con- 
centration (e.g, c/c 0296 in Fig. 31), bubbles can 
grow (or shrink) and approach a stable diffusive equilib- 
rium. Since Rayleigh collapse occurs at large P4, stable 
SBSL occurs here, with a well-defined, stable Rp follow- 
ing from P, and the gas concentration. For larger gas 
concentrations, large enough bubbles will again grow up 
to the threshold of parametric shape instability where 
microbubbles pinch off. Here, in contrast to the smaller 
P, regime, the remaining bubble is very likely to end up 
in a regime where it can grow again. 

The characteristically slow growth of Ro (over the 
time scales of rectified diffusion) and sudden breakdown 
(at microbubble pinch-off) are reflected in other experi- 
mentally observable parameters as well, such as the 
phase of light emission (with respect to the driving). Ex- 
perimental measurements of the phase are presented in 
Fig. 4 of Barber et al. (1995), showing the pinch-ofi/ 
growth dynamics. From the diffusive theory outlined in 
Sec. IV.B, Hilgenfeldt et al. (1996) simulated this behav- 
ior (sec Fig. 39), finding good agreement with the ex- 
perimental result, 

The momentum of pinched-off microbubbles also 
gives the remaining bubble a recoil. As this repeats 
again and again on the diffusive time scale of ~0.1 s, the 
bubble seems to “dance,” as originally observed by 
Gaitan (1990) and later by Barber et al. (1995). The re- 
gime of dancing bubbles is indicated in the experimental 
phase diagram Fig. 12. If the bubbles in this regime are 
large enough (close to the instability line), they will also 
emit sonoluminescence light even as they undergo recti- 
fied diffusion, leading to the same pattern of slow in- 
crease and sudden breakdown in the light signal. This is 
known as (diffusively) unstable SBSL and is generally 
fainter than stable SBSL. 

For very large driving pressures, the Rayleigh-Taylor 
shape instability (and possibly the Bjerknes instability) 
make stable bubble oscillations impossible for small Ro. 
This sets the upper limit in P, for the observation of 
SBSL bubbles. 

‘After computing phase diagrams like Fig. 38 with 
many different dissolved gas concentrations, the results 
can be summarized in a new plot whose variables are the 
experimentally controllable parameters c« and P,. This 
phase diagram in c..-P,, phase space is shown in Fig. 9. 
‘The notation in that diagram is as follows: stable SL or 
unstable SL means that there are glowing bubbles of 
certain ambient radii which are diffusively stable or un- 
dergoing rectified diffusion, respectively; other, smaller 
bubbles dissolve. Only a small, crescent-shaped region in 





462 Brenner, Hilgenteldt, and Lohse: Single-bubble sonoluminescence 


this phase space allows for stable SBSL. Luckily, the nar- 
row range of extremely small argon concentrations nec- 
essary for stable SBSL is easily achievable by working 
with gas mixtures such as air, because the molecular con- 
stituents of air dissociate (Sec. IV.C). Experimental 
phase diagrams in the c, and P, phase space can be 
found in Simon et al. (2001). 


F. Other liquids and contaminated liquids 


Other parameter dependencies of SBSL involve 
changes in the material parameters of the liquid. It was 
already mentioned in Sec. IV.D.5 that increasing liquid 
viscosity allows for stable bubbles at larger Rọ. Like- 
wise, a change in surface tension can influence the loca- 
tion of the shape instability line in the phase diagrams 
(see Hilgenfeldt, Brenner, er aL, 1998), For nonsonolu- 
minescing bubbles, Asaki and Marston (1997) have ex- 
perimentally examined the effect of surfactants on both 
shape oscillations and dissolution rates; Fyrillas and 
Szeri (1995, 1996) supply theoretical understanding of 
the surfactant effect. 

A change of the liquid in which a bubble oscillates has 
profound consequences on possible SBSL light emission, 
largely because of the chemical reactions occurring in- 
side the bubble. In Secs. IILA.3 and IILB2, we dis- 
cussed the influence of water vapor invading the bubble 
on the bubble temperature and the resulting light emis- 
sion, finding that more water vapor leads to less light. 
Moreover, dissociation products of the liquid (water) are 
crucial for the radical reactions that remove molecular 
gases from the bubble (Sec. IV.C). 

For quite some time, it was thought that stable SBSL 
could be achieved in water only, and even today water is 
still considered the most "friendly" liquid for SBSL ex- 
periments, While dissolved salts (Matula et al., 1995) or 
mixtures of water and freely miscible liquids (Gaitan 
et al., 1996) do not prohibit stable SBSL, for years it 
could not be observed in other liquids, even though mul- 
tibubble sonoluminescence in nonaqueous liquids had 
been known for a long time (see, for example, Suslick 
and Flint, 1987). 

Weninger et al. (1995) showed that weakly emitting, 
unstable single sonoluminescing bubbles could be ob- 
served in pure alcohols, while the light of a stable SBSL 
bubble in water could be turned off by just adding a few 
drops of alcohol to the solution, The latter effect can be 
understood by recognizing that alcohols are surface ac- 
tive and tend to accumulate at the bubble surface. At 
collapse, they enter the bubble and reduce the heating 
because of their smaller polytropic exponent and endo- 
thermic chemical reaction. Alcohols act much like water 
vapor in this respect, but are much more efficient in re- 
ducing the temperatures. Evidence for this kind of 
model has been gathered experimentally and theoreti- 
cally for both MBSL and SBSL (Ashokkumar et al, 
2000; Toegel, Hilgenfeldt, and Lohse, 2000; Grieser and 
Ashokkumar, 2001). 

While alcohol contaminations quench the light emis- 
sion quickly, the bubbles still oscillate in a stable fashion, 
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indicating that the general mass balance of the diffusive! 
chemical equilibrium is not severely disrupted. The dif- 
ference with the Weninger et al. (1995) experiment in 
pure alcohols thus lies in the solubility of the chemical 
reaction products in the surrounding liquid, which is still 
almost pure water in the case of alcohol contamination 
experiments. Didenko et al. (2000b) therefore sought 
ideal nonaqueous liquids for SBSL, requiring (i) a low 
vapor pressure to limit vapor invasion of the bubble and 
(ii) a high content of O or N heteroatoms to facilitate 
the chemical formation of species upon collapse that will 
readily dissolve in or react with the liquid phase. Using 
liquids such as formamide or adiponitrile, Didenko et al. 
(2000b) did produce bright SBSL and were the first to 
observe molecular spectral lines (see Introduction and 
Sec. V.G). Their bubbles were, however, unstable, in the 
sense that they were moving on circular or elliptic tra- 
jectories around the pressure antinode that held them 
trapped. Didenko et al. (2000b) call this state moving 
SBSL. Possibly the chemical reactions (of which the line 
emission gives direct evidence) still produce too much 
“waste” unable to dissolve fast enough in the surround- 
ing liquid, sending the bubble into an unstable state of 
rectified diffusion and microbubble pinch-off. Revealing 
the mechanism of moving SBSL is an interesting open 
problem for future research. 


V. SONOLUMINESCENCE LIGHT EMISSION 


The previous sections have described how a 
micrometer-sized bubble in a water-filled flask can un- 
dergo oscillations of incredible violence, collapse at su- 
personic speeds, burn molecular gases, and still—under 
the right experimental conditions—maintain the stability 
of its spherical shape, showing precise repetitions of this 
highly nonlinear dynamics for millions and billions of 
driving cycles. These remarkable properties might not 
have been studied in detail without the phenomenon 
that gives sonoluminescence its name: the emission of 
visible light induced by insonation with an acoustic 
wave. 

In order to release a photon of visible wavelength, an 
atom, ion, or molecule must be excited a few eV above 
its ground state. A sound wave of 1-atm amplitude, by 
contrast, carries an energy density of typically 10 !! eV 
per particle. The required tremendous energy concen- 
tration of almost 12 orders of magnitude (Barber and 
Putterman, 1991) is precipitated by the rapid collapse of 
the sonoluminescent bubble, where the layers of water 
surrounding the bubble act as a radial piston compress- 
ing its interior. But is this the whole story? Is the heating 
of the gas inside the bubble resulting from the rapid 
collapse sufficient to explain the light? Is the light emis- 
sion completely of thermal origin, or is it necessary to 
invoke other physical processes? While researchers have 
come up with a plethora of creative ideas concerning 
light-emission processes, the results of a number of cru- 
cial experiments favor the thermal light-emission ap- 
proaches over others, and rule out some other theories 
outright. 
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A. Theories of MBSL: discharge vs hot spot theories 


When multibubble sonoluminescence was discovered 
in the 1930s by Marinesco and Trillat (1933) and Frenzel 
and Schultes (1934), different theories for its occurrence 
were soon put forth. Almost all light-emission mecha- 
nisms discussed since then can be classified under one of 
two headings: thermal or electrical processes. 

‘The first attempts at explaining the mechanism behind 
the light emission favored electric discharges. Levshin 
and Rzhevkin (1937) initially brought up the subject of 
charge separation in cavitation bubbles; Harvey (1939) 
thought of the bubble as a spherical capacitor with 
charges at the center and the wall. Upon collapse, the 
capacitance decreases and voltage increases until elec- 
tric breakdown takes place. Frenkel’ (1940) suggested 
charge separation by enhancing charge fluctuations on 
the bubble wall, In this latter theory, however, break- 
down should occur during the expansion phase of the 
bubble dynamics (Leighton, 1994), whereas the close 
proximity of bubble collapse and light pulse has been 
firmly established by the work of Meyer and Kuttruff 
(1959). 

Since a symmetric charge distribution cannot radiate 
light, discharge theories in general have to assume that 
the emitting bubble undergoes an asymmetric collapse, 
and would predict increasing intensity of light emission. 
as the asymmetry increases (Margulis, 2000). This con- 
tradicts recent systematic studies of single-bubble cavi- 
tation luminescence (Ohl eral., 1998; Baghdassarian 
et al., 1999, 2000, 2001; Ohl, 2000), in which the bubbles 
are created with strong, focused laser pulses, and in 
which, above a certain collapse asymmetry, light emis- 
sion ceases altogether. It is also in contradiction with the 
observation that single sonoluminescing bubbles, which 
collapse under controlled conditions with high symme- 
try, emit light of much higher intensity than MBSL 
bubbles of comparable size driven at comparable levels, 
with the light of a single bubble easily visible to the 
naked eye (Gaitan eral, 1992). Another example in 
which an increase in symmetry leads to (slightly) more 
light rather than less is given by the SBSL experiments 
under microgravity performed by Matula (2000). 

‘The other large group of MBSL theories have been 
characterized as “hot spot” models, in which the energy 
for the light emission is supplied by thermal energy re- 
sulting from an adiabatic bubble collapse. Noltingk and 
Neppiras (1950) were the first to use Rayleigh-Plesset 
bubble dynamics to deduce bubble internal tempera- 
tures as high as 10000 K at collapse of a spherically 
symmetric bubble. Within the hot spot models, which 
process of light emission will dominate depends on the 
actual maximum temperatures reached, e.g., recombina- 
tion of dissociated molecules at lower temperatures 
(Saksena and Nyborg, 1970), or characteristic molecular 
radiation due to electronic excitation, in particular of the 
OH radical (Sehgal et al., 1980). The latter was referred 
to as chemiluminescence by Suslick (1990) and must not 
be confused with secondary chemiluminescence, which 
may occur in the liquid as a result of chemical reactions 
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of the radical molecules generated in the bubble col- 
lapse, One example of this is luminol, which emits light 
as it reacts with OH. While chemicals such as luminol 
are useful in detecting cavitation (Negishi, 1961), their 
emission is not directly related to sonoluminescence. 

In the past decade, Kenneth Suslick and his group 
have amassed impressive evidence in favor of thermal 
molecular-emission luminescence in multibubble cavita- 
tion fields in the context of sonochemistry. From the 
presence of clearly identifiable molecular bands and the 
absence of other lines associated with discharges, Suslick 
(1990) deduces the thermal origin of the emission, Iden- 
tifying line emissions by their location in the spectrum, 
Suslick er al. (1986) used comparative rate thermometry 
to evaluate the temperature from the intensity of the 
different lines, The temperatures obtained are very con- 
sistent and have been confirmed by different methods of 
spectral analysis by Flint and Suslick (191a) and Mc- 
Namara et al. (1999). Typical maximum temperatures of 
MBSL bubbles are thus determined to vary between 
—3000 and 6000 K, depending on experimental param- 
eters (McNamara et al., 1999). Another convincing fea- 
ture of thermal hot spot theories is that they naturally 
predict higher temperatures for collapses of higher 
spherical symmetry, which is confirmed by studies such 
as those of Ohl et al. (1998) and Matula (2000), as men- 
tioned above. 


B. SBSL: A multitude of theories 


With the advent of single-bubble sonoluminescence, a 
simple case study for cavitation and the resulting light 
emission was found. It was hoped that all open questions 
could be answered by studying this “hydrogen atom of 
cavitation physics.” Initially, the discovery of single- 
bubble sonoluminescence actually increased the confu- 
sion about sonoluminescence light-emission mecha- 
nisms. The main reason for the multitude of speculative 
models was the uncertainty whether SBSL was similar to 
MBSL. SBSL led to more light than a single MBSL 
bubble, while the spectra did not reveal any structures 
such as lines or bands (Fig. 6). Both of these facts (as 
well as the apparently short duration of the light pulse) 
suggested more extreme temperatures and pressures. 

Thus, early theoretical research tried to invent mecha- 
nisms that occur for more extreme conditions than 
MBSL. Some of these theories were quite exotic. For 
example, Schwinger (1992) hinted at the dynamical Ca- 
simir effect as a potential photon-emission process at the 
noninertially moving bubble interface. When a corre- 
sponding model was developed by Eberlein (1996a, 
1996b), it soon became clear that, in order to match the 
observed light intensities, the bubble-wall speed would 
not just have to be comparable to the speed of sound, 
but exceed the speed of light (Unnikrishnan et al., 1996; 
Lambrecht etal., 1997; Milton and Ng, 1998; Brevik 
et al., 1999). This particular line of research was thus 
abandoned. 

‘A number of theories placed the location of light 
emission in the liquid, rather than inside the bubble. In 
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this vein, another attempt at an electrical breakdown 
model was made by Garcia et al. (1999). Earlier, Lepoint 
et al. (1997) speculated on sparklike discharges around 
water jets invading a bubble. 

Prosperetti (1997) also invoked an electrical mecha- 
nism for light emission (fractoluminescence) as a by- 
product of a fluid-mechanical picture of sonolumines- 
cence light emission that requires asymmetric collapse of 
the bubble. His idea was based on the fact that Bjerknes 
forces cause a sonoluminescing bubble to oscillate verti- 
cally during a bubble cycle, and that such oscillations are 
well known to cause asymmetric collapse (jets). The col- 
lision of the jet with the bubble wall would initiate frac- 
toluminescence, an effect documented in solid-state ma- 
terials, Prosperetti argued that, over the small time 
scales of collapse and jetting, water could “break” as 
well. The model made various predictions, but it fell out 
of favor when it became clear that they do not hold. 
First, it had predicted that less light should be emitted 
under microgravity, where a weaker jet is expected. Yet 
Matula (2000) found the opposite. Second, the details of 
Prosperetti’s (1997) model specifically rely on the prop- 
erties of water, while Didenko et al, (2000b) found 
bright SBSL in nonaqueous liquids. Finally, it had mean- 
while become clear that models based on a spherically 
symmetric collapse were able to give a more quantita- 
tive explanation for the experimental data. 

Other SBSL theories focused on a more ordinary ex- 
planation, the emission of photons due to the high tem- 
peratures in the bubble, akin to the hot spot theory of 
MBSL. In fact, all nonthermal models have to explain 
why their mechanism of emission would not be 
swamped by thermal radiation, Requiring spherical sym- 
metry of collapse for intense light emission, thermal 
models of SBSL are favored even by proponents of dis- 
charge models for MBSL (Margulis, 2000). In an inter- 
esting twist, Xu et al. (1999) have suggested that electric 
fields do not directly lead to light emission, but suppress 
the Rayleigh-Taylor instability (see Sec. IV) and thus 
stabilize the bubble. However, the existence of large 
enough fields for such stabilization is a subject of debate 
(Moss, 2000). 

Depending on the actual temperatures achieved dur- 
ing collapse, different excitations become dominant in 
the compressed gas, so that “thermal emission” can refer 
toa large variety of different processes. As temperatures 
increase from several hundred to many thousand kelvin, 
those processes can be, among others, molecular recom- 
bination (Saksena and Nyborg, 1970), collision-induced 
emission (Frommhold and Atchley, 1994), molecular 
emission (Didenko etal., 2000b), excimers (Hammer 
and Frommhold, 2001), atomic recombination (Hilgen- 
feldt et al., 1999b), radiative attachment of ions (Ham- 
mer and Frommhold, 2001), neutral and ion bremsstrah- 
lung (Moss et al., 1997; Xu et al., 1998; Hilgenfeldt et al., 
1999b), or emission from confined electrons in voids 
(Bernstein and Zakin, 1995). The uncertainty about the 
precise temperatures of SBSL bubbles (see Sec. III) 
would allow for most of these mechanisms. Generally 
speaking, however, very-low-energy excitations are un- 
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likely to produce a spectrum with enough visible pho- 
tons to account for SBSL. On the other hand, the bubble 
collapse cannot generate arbitrarily high temperatures, 
so that very-high-energy excitations are rare and do not 
produce a large photon flux either. 

Which of these mechanisms actually dominates de- 
pends critically on accurate measurements and calcula- 
tions of the temperature inside the bubble. For example, 
at temperatures above several thousand kelvin, pro- 
cesses like collision-induced emission lose importance as 
higher-energy emission processes such as bremsstrah- 
lung take over (Frommhold, 1998; Hammer and Fromm- 
hold, 2000a, 2000b, 2001). Calculations have shown that 
the emission intensity from a body at a temperature of 
~10* K is not inconsistent with SBSL observations; al- 
though radiative enhancement effects (e.g., collective 
emission as proposed by Mohanty and Khare, 1998) 
have not been strictly ruled out, they are not necessary 
to explain the experimental data. It should be empha- 
sized here that although all thermal processes can con- 
tribute to blackbody radiation, thermal emission does 
not necessarily result in a blackbody spectrum (as elabo- 
rated upon in Secs. V.D and VE). 

For many years, the most serious argument against 
thermal processes was the sonoluminescence pulse 
width. As measured, it was much shorter than the time 
for which the bubble is maximally collapsed (~1 ns, see 
Fig. 27). The need for ultrashort pulses fueled the popu- 
larity of “shock-wave” models, in which a focusing shock 
causes the light-producing region to be much smaller 
than the bubble size (see the Introduction). This shock 
picture (referred to as the mechanochemical mechanism 
in Leighton, 1994) was the most popular view of SBSL 
for quite some time. 


C. Narrowing down the field 


The critical event that narrowed the field was the ex- 
perimental resolution of the light pulse, showing its 
width to be much longer than previously anticipated. 
Gompf et al. (1997) used time-correlated single-photon 
counting, a powerful method from solid-state physics 
and biophysics (O'Connor and Phillips, 1984), where it 
is used to register fluorescence lifetimes of sometimes 
only a few picoseconds. Using a modification of the clas- 
sical setup, Gompf et al. (1997) employed two photomul- 
tiplier tubes to record the arrival times of single photons 
from the same SBSL light pulse in both detectors. From 
the autocorrelation function of the time differences be- 
tween the registered events, they could reconstruct the 
shape (the temporal variation of intensity) of the light 
pulse and thus measure its duration. Gompf and co- 
workers confirmed these results using a direct streak 
camera measurement (Pecha et al., 1998). Furthermore, 
both publications demonstrate that higher driving pres- 
sures lead to light pulses of both higher intensity and 
longer width. Their results were confirmed by Moran 
and Sweider (1998) and Hiller et al. (1998), the latter 
noting that xenon bubbles, long known to be the bright- 
est SBSL emitters, also yield the longest pulses, with 
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widths up to 350 ps and more. Only for bubbles in highly 
degassed water can very faint flashes of about 40-60 ps 
duration be found. 

We have seen before that the bubble dynamics near 
minimum radius occurs on a time scale (“turnaround 
time”) of ~1 ns (Fig. 27). While the typical pulse width 
of 100-200 ps is still much smaller than this value, the 
gap between the light-emission time scale and the radial 
dynamics time scale is not nearly as wide as once 
thought. Gompf et al. (1997) demonstrate that a simple 
thermal model can result in a light pulse with the essen- 
tial characteristics they measured. This made more cred- 
ible all thermal models with heating processes directly 
connected to bubble dynamics. Conversely, exotic pro- 
cesses (yielding ultrashort pulses) were no longer neces- 
sary, It was then that shock-wave theories underwent a 
closer scrutiny with respect to their consequences for 
light emission. In Sec. III, we discussed how calculations 
showed that shock waves for noble gases would be far 
less intense than for molecular gases (Moss et al., 1997) 
or might be altogether absent (Vuong et aL, 1999). The 
nonrobustness of the shock-wave phenomenon makes it 
an unlikely candidate for the heating necessary to gen- 
erate thermal SBSL emission. In the work of Bourne 
and Field (1991), a shock wave was detected in an asym- 
metrically collapsing, mm-sized air bubble. However, the 
light emission detected at collapse did not originate 
from the shock wave, but from the compressed gas. All 
this evidence indicates that, while shock waves may be 
present in collapsing bubbles under certain circum- 
stances, they do not substantially contribute to SBSL 
light emission. This can be because either (i) they are 
too weak, so that they do not focus energy to a scale 
much smaller than the bubble size, or (ii) they can only 
heat a tiny volume near the very center of the bubble to 
very high temperatures, and the ensuing light emission is 
much less intense than the thermal radiation from the 
bulk volume of the bubble. 


D. The blackbody model and its failure 


‘The seminal publication of Gompf et al. (1997) pro- 
vided more valuable information: using wavelength fil- 
ters, they determined the widths of the light pulse in 
different parts of the optical spectrum. The results 
showed (Fig. 13) that there was hardly any detectable 
difference between the pulse width in the ultraviolet 
part of the spectrum (300-400-nm wavelength) and in 
the red part (590-650 nm). This contradicted thermal 
models favoring blackbody emission, These predict that, 
since the bubble heats during collapse and cools during 
expansion, it maintains somewhat lower temperatures 
for a longer time than higher temperatures. It is there- 
fore capable of emission at long optical wavelengths for 
a longer time than at short wavelengths (the latter re- 
quiring higher temperatures). Translating typical tem- 
peratures achieved in models of collapsing bubbles into 
blackbody emission shows that the “red pulse” should 
be about twice as long as the “ultraviolet pulse” (see 
below). Again, Hiller etal. (1998) and Moran and 
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FIG. 40. Blackbody spectra with the indicated temperatures 
fitted to the data from Barber et al. (1997) for the noble gases: 
C), helium; O, neon; X, argon; +, krypton; and @, xenon. The 
best-fit temperatures decrease dramatically with increasing 
molecular weight of the gas. From Hammer and Frommhold 
(2001). 


Sweider (1998) confirmed these findings, although the 
latter group did find a slightly longer pulse in the red 
than in the UV if cooled water was used. 

‘Adiabatic heating does not necessarily result in black- 
body radiation. Even though a model of pure blackbody 
emission cannot account for the experiments, it is help- 
ful to outline this simplest of all thermal light-emission 
theories in the context of SBSL to introduce basic con- 
cepts and to point out where it fails. A blackbody calcu- 
lation is of unrivaled simplicity because there is no need 
to specify the thermal emission as resulting from one of 
the many excitation processes listed Sec. VB—only the 
temperature matters. 

A blackbody of a given temperature emits a spectral 
light intensity (energy per unit time, wavelength inter- 
val, solid angle, and projected surface area) at wave- 
length of 

2he* 
XSTexpicIkkT)- 1] 
the Planck intensity, with the Planck and Boltzmann 
constants } and ky, and the speed of light in vacuum c. 
Experiments measure the spectral radiance (emitted en- 
ergy per time and wavelength interval), for which we 


must integrate over the projected bubble surface and all 
solid angles: 


PE) Am RP TT]. (66) 


This formula reflects the fact that a blackbody is a sur- 
face emitter, i.e., its radiation comes from its surface 
alone, as photons originating in deeper layers are ab- 


if'tr]- (65) 
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FIG. 41. Theoretically calculated light emission from a Ry 
=5 um SBSL bubble driven at P,=13 atm and f=20 kHz 
under the assumption of blackbody radiation, using the tem- 
perature dynamics T(r) from Fig. 27 and Eq. (66): dashed 
lines, pulses in the ultraviolet (P,,); dotted lines, pulses in the 
red (P,) wavelength range; (a) absolute instantaneous powers; 
(b) pulses normalized to a peak value of 1 for better compari- 
son of their widths; solid line, pulse integrated over the entire 
range of detectable wavelengths (P4). From Hilgenfeldt et al. 
(19995). 


sorbed within the blackbody. It can be compared directly 
with the measured spectra of sonoluminescence, as the 
only quantities that matter are the bubble temperature 
T(t) and radius R(t). Figure 40 shows a comparison 
(Hammer and Frommhold, 2001) of data from Barber 
etal. (1997) with blackbody calculations, The spectra 
match the blackbody law well, as the temperature varies 
from 8000 K to 30000 K, depending on the inert gas. 

‘The curves in Fig. 40 result from Eqs. (65) or (66) with 
both intensity and temperature as free fit parameters. To 
put SBSL models to the test, these parameters have to 
be linked to the dynamics of bubble radius and tempera- 
ture resulting from the model. Indeed, since there is no 
direct method for measuring the temperature, this is the 
only way to correctly deduce it. Unfortunately, we shall 
see that small differences in modeling can lead to large 
differences in light intensity. 

As outlined in previous sections, there are various ap- 
proaches to obtaining bubble dynamics and tempera- 
ture, with varying degrees of sophistication and accu- 
racy. A simple starting point is to assume a spatially 
homogeneous temperature T(r) (see Sec. III-B). Specifi- 
cally, we choose Eq. (38) to supplement the radius dy- 
namics R(t) derived from the modified Rayleigh-Plesset 
Eq. (20). The time dependence of R,T results, via Eq. 
(66), in a dynamics of the emitted power P/'dA in the 
wavelength interval [A,\+dQ], ie., the temporal shape 
of the SBSL pulse as viewed through filters of different 
color. 

Experiments observe wavelengths in the detectable 
range A, «A «A,, where Au, 200 nm is the ultraviolet 
cutoff of the visible spectrum (due to the strong absorp- 
tion of smaller-wavelength light in water or other liq- 
uids), and 4,~800 nm marks the red end of the visible 
spectrum, where the spectral radiance is already quite 
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small and experimental noise begins to obscure the spec- 
trum. In experiment, filters of a certain bandwidth AX 
are used, e.g., AN~100 nm by Gompf et al. (1997) and 
AA~40 nm by Moran and Sweider (1998), or replaced 
by spectrographic means for higher resolution (Hiller 
et al., 1998). It is easy to integrate Eq. (66) over the 
corresponding wavelength ranges and compare the pre- 
diction with light pulse measurements such as those in 
Fig. 13. 

Figure 41 shows the result of this calculation for the 
argon bubble whose dynamics and transport behavior 
we have discussed previously (driven at P,—1.3 atm, 
ambient radius Rọ=5 ym). Using R(t), T(t) from Fig 
27, the emission intensities are calculated for the ultra- 
violet wavelength regime (200 nm«A «300 nm), for the 
red part of the spectrum (700 nm<)<800 nm), and for 
the complete detectable spectral range (200 nm<A 
«800 nm). Figure 41(a) shows the predominance of ul- 
traviolet emission. 

‘A calculation like this allows for comparison of at 
least three characteristic quantities with experimental 
data: 


(a) The intensity of the pulse. In this model, it is about 
two orders-of-magnitude larger than the experi- 
mental values for these driving parameters. 

(b) The duration of the pulses. These are comparable 
to the turnaround time of bubble dynamics [not 
surprisingly, as R(t) is directly translated into 
PË'(t)], and is thus at least a factor of 2 longer than 
experiments. 

(©) The wavelength dependence of the pulse width 
The length of the pulse varies dramatically with 
wavelength, in direct contradiction to the experi- 
ments of Gompf et al. (1997). The “red” pulse du- 
ration is over 1000 ps, and thus more than twice 
that of the UV pulse. 


Because of the considerable modeling uncertainties in 
bubble temperature (see Sec. IIT), the absolute values of 
pulse intensity and width produced by this model should 
not be considered quantitative predictions, Even though 
the discrepancies with experiment are very large, one 
could still imagine that they could be due to defects in 
the temperature modeling T(r) or a change in the size of 
the light-emitting region (as in some shock models) 
However, the large wavelength dependence is an intrin- 
sic problem with the blackbody model that is indepen- 
dent of the temperature law or the absolute value of 
temperature. Therefore experiments conclusively show 
that it is necessary to seek a modification of the black- 
body model. 


E. The SBSL bubble as thermal volume emitter 


Which of the assumptions of the blackbody model 
fails? One basic prerequisite is that the radiation be of 
thermal origin, with 7(r) a well-defined temperature, 
ie. local thermodynamic equilibrium must hold even 
over the short time scales of the bubble dynamics at 
collapse, which seems doubtful at first sight. However, 
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FIG. 42. Calculated power of light emission: dashed line, with 
plasma thermal conduction; dotted line, without plasma ther- 
mal conduction; short-dashed line, boundary line of the opti- 
cally thick region of a nitrogen bubble; heavy solid line, the 
bubble radius; thin line, the shock front in the model. From 
Moss et al. (1997). 


the immense particle densities (n~10% m?) and high 
temperatures (7—10* K) at bubble collapse create an 
environment in which collisions between particles are 
very frequent: a simple estimate suggests collision times 
well below a picosecond, so that local thermodynamic 
equilibrium is still well obeyed during SBSL light emis- 
sion. 

The other crucial assumption of the blackbody picture 
is that the bubble is black, i.e., it perfectly absorbs all 
wavelengths of electromagnetic radiation. This requires 
that the mean free path of photons (x!) be much 
smaller than the size of the object. Here x, is the ab- 
sorption coefficient for photons of wavelength à (when 
divided by the gas density, x, is often referred to as the 
opacity; see Unsóld and Baschek, 1991). A blackbody of 
extent (radius) R has x, R?-1. The product 7, 92x,R is 
the (dimensionless) optical thickness of the object at A. 

That opacity may play a role for sonoluminescence 
light emission had first been noted by Kamath et al 
(1993), who speculated that radicals in the outer layers 
of the bubble “could resonantly absorb and scatter the 
radiation coming from the center thus shutting off the 
observed light. This effectiveness of resonant scattering 
in rendering gases effectively opaque is well known." 
Wu and Roberts (1993) also speculated that the com- 
pressed bubble might be a thermal volume emitter. In- 
dependently, Moss et al. (1994, 1997) discussed the pos- 
sibility that the SBSL bubble might be transparent to its 
‘own photons. Using tabulated absorption coefficients of 
different gas species, they concluded that there was only 
a small region of opaque gas in the very center of the 
bubble (induced in their model by a converging shock 
wave), but that the rest of the bubble volume was char- 
acterized by 7, «1, and was thus optically thin (transpar- 
ent). Figure 42 shows the line of demarcation between 
the optically thick and optically thin parts of the bubble. 
‘The emission volume of the optically thin shell is much 
larger than that of the optically thick core. 
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Moss et al. (1997) did not explicitly calculate pulses of 
emitted light at different wavelengths, and other fea- 
tures of their model lead to pulses that are much shorter 
than observed. However, the idea of a transparent 
bubble is the only proposition put forth to date that ex- 
plains the wavelength independence of sonolumines- 
cence radiation, and it is now widely accepted as a key 
ingredient to a consistent view of SBSL light emission, 
Models built on bubble transparency are sometimes re- 
ferred to as “weakly ionized gas models” (Hammer and 
Frommhold, 2001), because a low degree of ionization 
of the gas is crucial (see below). We prefer the term 
thermal volume emission, as the radiation from the 
whole volume of a transparent body reaches the detec- 
tor, rather than only the surface emission of a blackbody. 

To illustrate that such a model leads to wavelength- 
independent emission spectra, we revert to the simple 
model discussed below, drawing R(t) and T(t) from 
Eqs. (20) and (38), respectively. For a volume emitter, 
the light intensity from a volume element of the bubble 
is dependent on the location inside the body (the depth 
s along the ray of vision to the element). Using the laws 
of absorption and emission in a medium (see, for ex- 
ample, Zel'dovich and Raizer, 1966, or Siegel and How- 
ell, 1972), 


As DRITAN- expl- e ETT. 
0<s<2R (60) 


is the intensity at wavelength A from a depth s, provided 
the temperature of the body is spatially uniform. Note 
that the source function of the emission is still /7!, be- 
cause local thermodynamic equilibrium is obeyed, and 
that for infinite absorption coefficient x,, the Planck 
emissivity is recovered. 

To obtain the spectral radiance, one must first inte- 
grate the intensity per surface area and per solid angle 
from Eq. (67) over the projected cross section of the 
bubble, the thickness s of the medium varying as s 
=2R cos(#) with the angle @ from the line of view. The 
result is the emitted power from the whole bubble crass 
section per unit solid angle. Assuming isotropy of SBSL 
radiation, the spectral radiance is then 

POA =P RETO qt On eR) 
mR 


exp(—2«,R)-1) 
MEET Jar (68) 


Note that the Planck radiance (66) is now multiplied by 
a factor whose magnitude varies between 0 and 1, de- 
pending on the optical thickness ,=2x,R. As ny, 
the Planck spectrum is recovered. In the transparent 
limit (741) the radiance reduces to 


thi 2 PI 
PY” = s s, PPLT() dh (69) 


a spectral power that is much smaller than that for 
blackbody emission. 
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FIG. 43. Dominant photon absorption and emission processes 
at the typical temperatures and densities in a sonoluminescing 
bubble: left, absorption; right, emission. The three pairs of pro- 
cesses are (i) electron-ion (inverse) bremsstrahlung, (ii) 
electron-atom (inverse) bremsstrahlung, and — (ii) 
photoionization/radiative recombination. Note the large gap 
between the ground-state energy (Egy) and the first excited 
state (E,,,) in the term scheme of (ii). 


1. Simple model for bubble opacity 


To proceed further, we need to know how the absorp- 
tion coefficients x, depend on parameters. In contrast to 
the blackbody calculation, this requires identification of 
the most significant physical processes that contribute to 
photon absorption in the gas inside a sonoluminescing 
bubble. Now not only the temperature is relevant, but 
also the physics of photon-matter interaction. The origi- 
nal work of Moss er al. (1997) extrapolated tabulated 
values for x, to the regime of the sonoluminescence ex- 
periments. Moss et al. (1999) used a full opacity model 
to extend these calculations. Here we summarize a 
simple model for «, (Hilgenfeldt et al., 19992, 1999b; 
Hammer and Frommhold, 2001). 

The model assumes that light emission predominantly 
stems from the ionization of a noble gas, and that nei- 
ther molecular gases nor liquid vapor play a substantial 
role. If it is further assumed that the peak temperature is 
of order —10000K, the literature on absorption and 
emission of radiation (e.g. Zel'dovich and Raizer, 1966) 
suggests that there are three important processes, illus- 
trated on the left-hand side of Fig. 43: the absorption of 
photons by (i) free electrons near ions (inverse brems- 
strahlung), (ii) free electrons near neutral atoms (in- 
verse neutral bremsstrahlung), and (iii) bound excited 
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electrons that reach continuum energies after photon 
absorption (i.e, photoionization). 

All three processes depend critically on the density of 
free electrons in the gas, which can be obtained from the 
Saha equation (Zel'dovich and Raizer, 1966; Unsöld and 
Baschek, 1991). At a temperature of T~10* K, a noble 
gas such as argon is only weakly ionized. This is because 
kgT<Ejon, the ionization energy of argon being Eion 
~158 eV. Zel'dovich and Raizer (1966) show that the 
degree of ionization a=n,/n is then 





exp = ZksT, 
where m, is the electron mass and u- „g are the statis- 
tical weights for the ionic and the neutral ground states, 
respectively. For argon SBSL bubbles æ typically does 
not exceed 196. 

Note that, for the three processes mentioned above, 
the frequency of the electron-ion reaction (i) is propor- 
tional to a?, while the two others are proportional to a 
alone. A small æ thus seems to favor (ii) and (iii) as 
dominant, but the interaction cross section of electrons 
and ions is much larger than that for electrons and neu- 
tral atoms, which makes (i) an important factor for pho- 
ton absorption. 

Given a, the absorption coefficients for the processes 
(1)-Gii) can be calculated (Zel'dovich and Raizer, 1966) 
for ionized noble gas atoms. The evaluation of x, for 
(ii) in principle requires knowledge of the complete 
atomic levels of the noble pas atom. Since the energy 
required for excitation to the first excited level is already 
a large fraction of the ionization energy (Fig. 43), the 
higher levels can be roughly approximated as part of a 
continuum of a hydrogenlike atom. In this limit, the only 
important parameter is the energy E of the first excited 
level or, equivalently, the corresponding photon wave- 
length \2=hc/E. Following Zel'dovich and Raizer 
(1966), one can show that the sum ic^" of the contribu- 
tions to x, from (i) and (iii) is 
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The contribution of electron-neutral atom inverse 
bremsstrahlung (ii) is rarely important in conventional 
plasma physics, where the degree of ionization is often 
high. For a weakly ionized SBSL gas, with the assump- 
tions made above, it is often the dominant factor in x, 
and can be written as (Zel'dovich and Raizer, 1966; 
Hilgenfeldt et al, 1999b) 

e? (QkyT) n? 


r= 
miT]-4 Ame hm a aM 


en) 


2 


Eion 
ZkaT) 





{ da \ | 
x| cnt SET) exp — N (72) 


Brenner, Hilgenfeldt, and Lohse: Single-bubble sonoluminescence. 469 


power/mw 
2° 


norm. power 
5 





FIG. 44. Theoretically calculated light emission from a Ro 
=5 um SBSL bubble driven at P,—13 atm and f=20 kHz 
under the assumption of thermal volume emission. We use 
T(t) from Fig. 27, and Eq. (68) to evaluate the emission. (a) 
Absolute power; (b) relative power; line styles as in Fig. dl. 
The strong dependence of the pulse width on wavelength dis- 
played in that figure is not present here, in agreement with 
experiments. The total power in the pulse is now much smaller. 
From Hilgenfeldt et al. (19990). 


Here, c, and d,, are coefficients describing the electron- 
atom transport scattering cross section, extracted from 
Brown (1966). The accuracy of Eq. (72) is limited by its 
assumptions. Hammer and Frommhold (200b, 2001) 
have evaluated the electron-neutral contribution using 
an ab initio quantum-mechanical calculation of the 
electron-atom scattering problem. At long wavelengths 
these calculations agree with the simple theory pre- 
sented above, but at short wavelengths the quantum cal- 
culations are significantly more intense. The two ap- 
proaches will be compared below. 


2. Light emission and comparison with experiment 


From Eqs. (71) and (72), the total absorption coeffi- 
cient x, — kQ- i?" is obtained. The dimensionless opti- 
cal thickness satisfies , «1 throughout the collapse for 
all wavelengths, and 7,-€1 for the dominant ultraviolet 
part of the emission. Figure 44 presents the emission 
from such a bubble according to Eq. (68), using the 
same parameters that resulted in the blackbody emission 
displayed in Fig. 41. 

Comparing with the blackbody calculation, it is seen, 
first, that the total intensity is reduced dramatically (to a 
few 10° photons per pulse). This is because the small 
optical thickness, by Eq. (69), leads to a drastic reduc- 
tion of radiance compared to an opaque bubble. Second, 
the light pulses are now considerably shorter, and finally, 
the dependence of pulse width on wavelength has al- 
most disappeared. All of these properties of the light 
pulse are now in much better agreement with those of 
experimentally observed SBSL pulses. Both the shorter 
pulse width and the wavelength independence are con- 
sequences of the exponential dependence of the degree 
of ionization on T in Eq. (70), which carries over to the 
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FIG. 45. Theoretical SBSL spectra for a Ro=5 um argon 
bubble driven at P,—1.3 atm and f=20 kHz: heavy solid line, 
total spectrum calculated from first principles; thin solid line, 
spectrum calculated from the semiclassical model outlined in 
Sec. VE; heavy and light dashed line, contribution of electron- 
neutral bremsstrahlung in the two models (more intense in the 
ab initio theory); dotted lines, electron-ion interaction contri- 
bution, the same for both approaches. From Hammer and 
Frommhold (2000a).. 


SBSL radiance. Since kgT<Ejon, this exponential 
“switch” is extremely sensitive to changes in T. Once 
the temperature has dropped only slightly as the bubble 
starts to reexpand, æ and x, are reduced severely. Con- 
sequently, the emission (69) is quenched rapidly (pro- 
ducing shorter pulse widths), independent of wavelength, 

Equation (68) also yields reasonable results for the 
spectrum of the emission, although the spectra have a 
“kink” due to modeling assumptions of the atomic en- 
ergy levels. Hammer and Frommhold (20002) have ex- 
tended the above calculations, using ab initio quantum 
mechanics. These calculations give a spectral shape 
more like the experiments for many parameter combi- 
nations (Fig. 45), with generally more intense radiation 
at the UV end of the spectrum, due to more intense 
electron-neutral bremsstrahlung. 

It is also useful to examine the predictions of this 
model when the gas is changed from argon to, say, xe- 
non. The spectrum is now much more intense (in agree- 
ment with observation) and has a maximum at Ao(Xe) 
~336 nm, quite close to what is measured (Barber et al., 
1997). For xenon, with its lower ionization energy 
Ej, (Xe) 12.1 eV, the degree of ionization o can be as 
large as 1096 (the temperatures obtained in argon and 


470 Brenner, Hilgenfeldt, and Lohse: Single-bubble sonoluminescence. 





10 15 20 25 30 35 
rel. intensity 











00 04 08 12 18 20 24 28 32 
rel. intensity 





FIG. 46. Comparison of the theoretical results from the ther- 
‘mal volume emitter model of Hilgenfeldt et al. (199b) to the 
experimental results of Hiller ef al. (1998). The figure was 
modified from Hilgenfeldt et al. (19992) to include the experi- 
‘mental results for helium (+). Filled symbols represent experi- 
ments using xenon (8), argon (A), and air (€) at the satura- 
tion concentrations displayed in the figure. Open symbols are 
theoretical results of the parameter-free theory of Hilgenfeldt 
et al. (19996). (a) Higher light intensities and pulse widths of 
xenon bubbles; (b) lower-intensity data. 


xenon bubbles are comparable). Consequently, the 
bubble is less transparent and its light emission more 
intense and more similar to those of a blackbody. As 
quantitatively shown by Hammer and Frommhold 
(2000b), xenon spectra display the characteristic kinked 
maximum at A in the quantum-mechanical version of 
the theory. This typical shape of the xenon spectrum was 
first understood by Moss et al. (1999). 

Identifying the major contributions to the absorption 
coefficient automatically identifies the main light- 
emission mechanisms of SBSL. This follows from Kirch- 
hoff’s law, which states that every absorption process 
must balance a corresponding emission process—this is 
why P, in Eq. (68) grows with increasing absorption. 
Inverting the absorption events (i)- (ii) identifies SBSL 
light emission as a combination of (i) electron-ion 
bremsstrahlung, (ii) electron-neutral bremsstrahlung, 
and (iii) radiative recombination (see the right-hand side 
of Fig. 43). None of the three processes seems to be 
dominant over the whole range of the SBSL parameter 
space. 
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One possible absorption/emission process not dis- 
cussed here is line emission (bound-bound-state transi- 
tions of electrons). Two sets of recent experiments have 
demonstrated that, in certain parameter regimes, line 
emission is important, We discuss these very recent re- 
sults below in Sec. V.G (see also the Introduction). 

We shall now use Eqs. (68), (71), and (72) to try 
a direct comparison with experimental results. To this 
end, a calculation for a single-parameter combination of 
P, ,Ro is not useful, as the emission depends sensitively 
on these parameters, and their accuracy from experi- 
ment is not very high. However, a promising set of data 
was presented by Hiller et al. (1998), graphing the width. 
of the light pulses vs their intensity. Both of these quan- 
tities can be measured reliably and with good accuracy. 
The experiments are conducted at a given relative con- 
centration c/co of gas in the liquid, and the data points 
shown in Fig. 46 are a result of changing P, through the 
whole range of stable SBSL under these conditions. The 
simple character of the light-emission model of Hilgen- 
feldt et al. (1999a, 1999) allows one to first calculate the 
P, Ro values for all stable SBSL bubbles at the given 
c. lcg (see, for example, Fig. 31) and then use these to 
calculate the light-emission intensities and pulse widths. 
Like the experimental data, the calculated light intensi- 
ties were normalized by the intensities for a standard 
argon bubble. As P, is increased in experiment, we fol- 
low the prescribed increase of Ro for stable bubbles and 
consequently obtain larger, more strongly driven 
bubbles whose light pulses become longer and more in- 
tense. Without adjusting any fit parameters, the theoret- 
ical results are in close agreement with the data (Fig. 46) 
for both xenon and argon gases. Note that the observed 
close proximity of the width vs intensity curves for Xe 
and Ar [Fig. 46(b)] is reproduced, even though the P, 
and Ro are quite different for the two pases. 

The data of Hiller et al. (1998) provide another check 
of the dissociation hypothesis discussed in Sec. IV.C. 
One of the data sets was obtained with strongly de- 
gassed air, which results in very small pulse widths and 
intensities [diamonds in Fig. 46(b)]. If the dissociation 
hypothesis is correct, these bubbles should not be air 
bubbles at c./co=3%, but argon bubbles at c./co 
3%. The corresponding calculation indeed repro- 
duces the experimental results, 

The simple theory outlined above cannot reproduce 
the absolute light-emission intensities for every param- 
eter combination, This should not be expected, since im- 
portant factors such as the influence of water vapor on 
the bubble temperature are absent from this model. 
However, this kind of model does yield the correct rela- 
tive intensities and their dependence on pulse width 
(Fig. 46), reasonable spectra (Fig. 45), and wavelength- 
independent pulse widths (Fig. 44). This agreement with 
experiment suggests that the essential physics is ad- 
equately represented and the concept of a bubble as a 
thermal volume emitter is probably correct (cf. also the 
extensive review of light-emission mechanisms by Ham- 
mer and Frommhold, 2001), 
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Even so, it should be emphasized that the good agree- 
ment between the model and experiment does not de- 
finitively settle the issue because a number of significant 
effects have not been included, in particular (i) the pres- 
ence of water vapor and chemical reactions inside the 
bubble and (ii) the additional possibility of light emis- 
sion from the molecular components of the gas 
(Didenko et al., 2000b; Young et al., 2001); again, water 
vapor may play an important role here. These and other 
processes could act to suppress or enhance the light 
emission, and to our knowledge no quantitative theory 
has incorporated all of these influences. In the theory 
formulated above, the number of photons is in reason- 
able agreement with the experiments, but it is entirely 
possible that this is because several neglected effects 
cancel each other out. The following section outlines 
some additions to the model that could affect the out- 
come of a light-emission calculation, and could lead to a 
more comprehensive model of SBSL radiation. 


F. Modeling uncertainties: additional effects 


The number of photons per oscillation cycle of 
sonoluminescence is the most conspicuous experimental 
observable; however, it is also the hardest to predict, in 
particular within a thermal volume emitter theory. The 
reason is that the light intensity depends exponentially 
on the temperature in the bubble [see Eq. (70)], which 
in turn depends strongly on modeling assumptions. The 
situation is illustrated in Fig. 23. What can be measured 
with good accuracy are the radius R(r) and the radiance 
P, (t) of SBSL light. Those two quantities are connected 
via the temperature T(t) of the SBSL bubble, which up 
to now could not be experimentally measured. There- 
fore modeling is required both for obtaining T(t) from 
R(t) and for obtaining P, (t) from T(t). Both modeling 
steps depend on the assumptions entering the respective 
theories, and often the dependence is sensitive. We now 
review some of the most relevant sources of modeling 
uncertainties. 


1. Bubble hydrodynamics 


As discussed in Sec. III, any uncertainty in modeling 
the hydrodynamics and thus the temperature inside the 
bubble will be directly reflected in some uncertainty 
about the light emission. Since the work of Plesset and 
Prosperetti (1977), it has been popular to state that the 
temperature law is initially isothermal and then switches 
abruptly to adiabatic once the bubble accelerates suffi- 
ciently. There are slightly more sophisticated versions of 
this, e.g., the approximation used by Hilgenfeldt et al. 
(1999b) in which the polytropic exponent changes con- 
tinuously. Precisely where the temperature law changes 
from isothermal to adiabatic strongly affects both the 
width of the pulse (since heating commences when this 
transition occurs) and the total number of photons emit- 
ted (since the transition controls the total amount of 
time for which heating occurs). 

Storey and Szeri (2000) have performed full numerical 
simulations showing that these approximations are rea- 
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sonably reliable for typical parameter combinations. 
However, to date there has not been a complete study of 
the full parameter space, documenting quantitatively the 
error in the simple models for the temperature, Without 
such an error estimate (as a function of parameters), itis 
difficult to assess how much of the current uncertainty 
about light-emission intensities comes from simple hy- 
drodynamics, 


2. Water vapor as emitter and quencher of light 


As an example for the above-mentioned possibility of 
canceling errors of approximations connected with the 
two modeling steps in Fig. 23, we mention the role of 
water vapor in argon bubbles: In Sec. III we have shown 
how the consideration of vapor can drastically reduce 
the temperatures achieved inside the bubble, both be- 
cause of the reduced polytropic exponent and because 
of the endothermic water dissociation. Because of the 
lower temperatures, less light should result. On the 
other hand, the reaction products of water (O and H) 
have an ionization potential (7713.6 eV for both species) 
very similar to that of argon, and therefore will also con- 
tribute to the light emission, increasing the photon out- 
put again. 

SBSL experiments with helium bubbles illustrate the 
latter effect: The large ionization energy E 
724.6 eV, together with the exponential dependence of 
the photon absorption coefficient on the ratio of 
Ej, kgT, predicts negligible light emission from He, 
about four orders-of-magnitude dimmer than for Xe. 
Yet Barber et al. (1997) observed He bubbles to be only 
about a factor of 10 less bright than Xe (see data in Fig. 
46). To a lesser degree, this discrepancy occurs for neon 
as well. The data of Barber et al. (1997) demonstrate 
that there is a pronounced decrease in SBSL intensity 
from the heavier to the lighter noble gas species in the 
order Xe—Kr— Ar, but Ne and He are not much differ- 
ent (in fact, He spectra are sometimes more intense than 
Ne spectra). 

Moss (1998) pointed out that these observations can 
be accounted for when light emission from O and H 
radicals generated from dissociated water is considered. 
These species emit light much more readily than neon 
and helium, whose ionization energies are considerably 
higher. In a Ne or He bubble, the emission intensity of 
the dissociation products of water vapor overcomes the 
tiny contribution of the light noble gases and does not 
vary much between these two gases. Thus, as pointed 
out by Hilgenfeldt et al. (1999b), the light from Ne and 
He bubbles could originate from the vapor in the bubble 
rather than from the noble gas. The measurement of 
OH lines in the spectrum of a sonoluminescing bubble 
by Young et al. (2001) (see Sec. V.G below) supports this 
view, directly demonstrating light emission from water 
dissociation products. 

Section IILA.3 detailed the important role of water 
vapor as a quencher for SBSL light. Through changes in 
polytropic exponent and endothermic chemical reac- 
tions, it can reduce the bubble temperature considerably, 
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as demonstrated by Yasui (1997a), Storey and Szeri 
(2000), and others. In bubbles consisting of heavier 
noble gases, this effect more than outweighs the role of 
water vapor as an emitter of light. Vazquez and Putter- 
man (2000) have found a nice demonstration of this ef- 
fect, in which two SBSL bubbles display almost indistin- 
guishable R(r) dynamics (very similar Ro and P,), but 
differ in light intensity by a factor of about 5 due to the 
different water-vapor pressures at different ambient 
temperatures. 


3. Further difficulties in modeling the temperature 


‘Temperature calculations that take into account 
chemical reactions suffer from considerable uncertainty, 
because the reaction rates of water-vapor chemistry are 
not well known under the extremely high densities and 
pressures achieved in the bubble. For example, as dis- 
cussed in Sec. III.B.2 above, Toegel et al. (2002) took the 
reduction of the water dissociation rate caused by the 
high densities inside the bubble into consideration and 
got about 50% higher temperatures than Storey and 
Szeri (2000). 

Another relevant effect may be the segregation of dif- 
ferent species inside the bubble, which had been sug- 
gested by Storey and Szeri (1999). Yasui (2001) put forth 
this mechanism to account for the relative brightness of 
helium bubbles compared to what would be expected in 
a thermal volume emitter model with homogeneous dis- 
tribution of all species. In a bubble with helium and wa- 
ter vapor, the lighter noble gas is driven towards the 
bubble center. A higher temperature could then be ex- 
pected inside the bubble: First, the diffusive equilibrium 
condition changes, allowing for higher driving pressures; 
second, the water vapor accumulates near the (cool) 
wall of the bubble and therefore does not consume so 
much of the collapse energy by its endothermic dissocia- 
tion; third, the condensation of water out of the bubble 
is also facilitated by this mass segregation (in a helium 
bubble), further helping to keep temperatures high. 





4. Modifications of photon-emission processes 


AII the effects listed above concern modeling uncer- 
tainties on the R(t)— T(r) side of Fig. 23. On the other 
side, less work has been done to evaluate how much 
quantitative errors affect the simplifying assumptions in- 
troduced by the light-emission formulas (70), (71), and 
(72). The ab initio calculations of Hammer and Fromm- 
hold (2000a, 2000b) show encouraging agreement with 
this simpler model, but other effects could still play a 
role, in particular those associated with the extremely 
high densities in the compressed bubble. 

‘The Saha equation may require corrections due to 
many-body effects (Chihara et al., 1999), leading to dif- 
ferent degrees of ionization and different absorption co- 
efficients, Also, the term schemes and ionization poten- 
tials of the noble gas atoms could shift under a high 
density of electrons (Zel'dovich and Raizer, 1966). Ham- 
mer and Frommhold (2001) pointed out that the latter 
effect tends to enhance the light emission, because 
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electron-electron interactions decrease the ionization 
potential, Thus more atoms get ionized, resulting in 
more light. 

‘A quantitative modeling of the magnitude of these 
effects for the case of a sonoluminescing bubble has not 
yet been attempted. It would be valuable in order to 
assess where, in the two-step scheme of Fig. 23, the most 
important modeling errors are introduced. 


5. Towards a more comprehensive model of SBSL light 
emission. 


How can we overcome the various modeling uncer- 
tainties? From our point of view the only way is to test 
the models through detailed experiments in a large pa- 
rameter domain, varying not only the forcing pressure 
and the gas type and concentration (and therefore also 
the ambient radius), but also the water temperature and 
driving frequencies. Indeed, the experiments along these 
lines by Vazquez and Putterman (2000) at low tempera- 
ture and those by Toegel, Gompf, et al. (2000) at low 
frequencies have demonstrated the relevance of liquid 
vapor. 

Many more experiments of this type will be necessary 
to improve the models connecting R(t) and T(t), on the 
one hand, and T(r) and P, (r) on the other hand. Given 
the unusual conditions inside the bubble (in particular, 
the high pressure at low degree of ionization), this is a 
research direction whose benefits lie not only in explain- 
ing SBSL, but in elucidating the physical makeup of 
matter in a state that has not been studied in any detail 
before. Understood this way, the SBSL bubble is an 
ideal microlab for high-pressure/high-density hydrody- 
namics, chemistry, and plasma physics. The focus of the 
questions asked should perhaps shift from the nature of 
light-emission mechanisms (which, from our point of 
view, is now understood quite well) to more application- 
oriented problems. 

If this model is able to capture the essential features 
of SBSL, no exotic “new physics” is necessary to explain 
the phenomenon. With classical hydrodynamics, disso- 
ciation chemistry, thermodynamics, and the theory of 
light absorption and emission in hot gases, well-known 
results from diverse fields now seem to be fitting into a 
consistent whole. 


G. Line emission in SBSL 


Line emission in multibubble sonoluminescence has 
long been known (Suslick, 1990; Flint and Suslick, 
1991b); typically it originates from neutral hydroxyl 
radicals, or other molecular species present in the liquid. 
Until very recently, molecular lines had never been ob- 
served in single-bubble sonoluminescence. Two explana- 
tions were typically given for the absence of lines in 
single-bubble sonoluminescence: either the lines are 
broadened by the immense pressure inside the bubble 
and smeared over the whole visible spectrum, or the 
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continuum radiation is just much more intense than the 
lines, which are consequently “swamped” by the con- 
tinuum spectrum. 

Very recently, the groups of Suslick (Didenko et al., 
2000b) and Kang (Young et al., 2001) have discovered 
parameter regimes in which there is line emission from 
single sonoluminescing bubbles. These were measured 
in situations where the bubble tended (i) to be unstable 
(“moving SBSL," in which the bubble gyrates around a 
pressure antinode of the ultrasound field) and (ii) to 
emit light rather weakly, close to the sonoluminescence 
threshold. 

Didenko et al. (200b) generated SBSL in various or- 
ganic liquids such as adiponitrile, a liquid with very low 
vapor pressure (to keep the effective polytropic expo- 
nent high) and a chemical structure that results in disso- 
ciation products easily soluble in the liquid. In these re- 
spects, adiponitrile and other polar aprotic liquids 
resemble water, The choice of liquids by Didenko et al. 
(2000b) was guided by the theoretical findings that liq- 
uid vapor can be responsible for quenching SBSL (see 
Sec. IILA.3) and that the chemistry of dissociation prod- 
ucts plays an important role in bubble stability (see Sec. 
IV.C). 

Figure 19 shows the observed spectra from moving 
SBSL bubbles in adiponitrile. Didenko et al. (2000b) at- 
tribute the observed lines to excited species of CN, one 
of the groups that make up adiponitrile. It can be seen 
that the line-emission component of the spectra is very 
prominent at lower driving pressure amplitudes (and 
lower total emission intensities) and becomes undetect- 
able as the driving pressure and the light intensity in- 
crease, until only the continuum part of the emission 
remains. This observation favors the idea of “swamped” 
spectral lines: at lower forcing, with lower temperatures 
in the bubble, the thermal bremsstrahlung processes 
(generating a continuum spectrum) are still weak, and 
the easily excited characteristic molecular emission from 
CN is strong enough to outshine the continuum. Higher 
driving generates higher temperatures and more ioniza- 
tion in the bubble, leading to bremsstrahlung emission 
much brighter than the contributions from lines, yielding 
the appearance of a pure continuum as observed for 
stable SBSL bubbles. While it is not clear to what extent 
the instability of the moving SBSL bubble in this experi- 
ment influences the relative intensity of the spectral 
lines (as compared to stable SBSL), this work is ex- 
tremely important because it conclusively demonstrates 
that liquid-vapor chemistry occurs inside a single lumi- 
nescing bubble. 

Working independently, Young et al. (2001) found 
spectral lines for pure noble gases dissolved in water. 
The lines were present at low driving pressures, right at 
the onset pressure for sonoluminescence. The dissolved 
gas concentration was about 20%, which corresponds to 
unstable sonoluminescence for a pure noble gas, with 
“dancing” bubbles ejecting microbubbles (see Sec. 
IVB). For pure argon dissolved in water, stable SBSL is 
not possible at the very low light intensities that are of 
interest here, To collect enough photons, Young et al. 
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(2001) needed to develop methods for collecting reliable 
data from single bubbles for up to 5 days; the dimmest 
bubble they studied had an intensity more than three 
orders-of-magnitude smaller than a typical sonolumi- 
nescing bubble. 

For bubbles with very low light intensity, spectral 
peaks were seen at both 310 nm and 337 nm (see Fig. 
18). The 310-nm line is a clear signature of the 
vibrational/rotational bands in excited OH radicals. This 
experiment therefore provides the first direct experi- 
mental indication that there is water vapor in the 
bubble. The origin of the line at 337 nm is unclear; 
Young et al. (2001) conjecture that it could be due to 
another molecular excitation of OH. 

When Young et al. (2001) increased the driving pres- 
sure, they observed, just as had Didenko et al. (2000b), 
that the continuum part of the spectrum became stron- 
ger relative to the line emission. They determined that, 
although the relative intensity of the lines decreases with 
increasing driving, the absolute line intensity still in- 
creases, again indicating that the reason spectral lines 
are not typically seen in single-bubble sonoluminescence 
is that they are overwhelmed by the continuum. Young 
et al. (2001) also used different noble gases and showed 
that the intensity of the 310-nm line increased with the 
molecular weight of the noble gas (Fig. 2 of that work), 
so that xenon had a more intense 310-nm emission than 
argon (the 337-nm line, however, seemed to be unaf- 
fected by the gas species). By lowering the temperature 
of the water, Young et al. (2001) observed that the inten- 
sity of the spectral lines increased. At 5°C, additional 
lines appeared. 

Another crucial issue was to determine precisely why 
the spectral lines were only observed in unstable or 
moving single-bubble sonoluminescence. Young et al. 
(2001) tried to find spectral lines in stable sonoluminesc- 
ing bubbles, but did not succeed at the high driving pres- 
sures needed for stable SBSL, where the continuum part 
of the spectrum is always overwhelming. On the other 
hand, growth by rectified diffusion, followed by the ejec- 
tion of microbubbles due to shape instabilities (Hilgen- 
feldt et al., 1996) could lead to quantitatively different 
water-vapor content and different chemistry inside an 
unstable bubble as compared to a stable one. Again, a 
quantitative model of line emission would settle this 
question, predicting the range of parameters for which 
line emission can outshine the SBSL continuum. 

Recent experimental results of Baghdassarian et al. 
(2001) on the spectra of single-cavitation bubble lumi- 
nescence may help in developing such a model: Beyond 
a certain laser energy threshold, OH lines show up in the 
emitted spectra, The interpretation is that for large 
enough bubbles the collapse is aspherical, so that some 
fluid can enter the bubble. This interpretation of the 
single-cavitation bubbles thus provides a connection to 
the moving SBSL (Didenko etal, 2000b), unstable 
SBSL (Young et al., 2001), and MBSL (Matula et aL, 
1995) spectra, which can also show lines. 
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VI. SUMMARY AND OUTLOOK 


This review has presented current ideas about single- 
bubble sonoluminescence, emphasizing the most impor- 
tant physical processes that govern bubble behavior and 
light emission. A combination of different physical ef- 
fects are important for the dynamics and thermodynam- 
ics that act over different time scales and, therefore, at 
different times during the bubble’s oscillation. To sum- 
marize these effects, the beginning of this section follows 
an SBSL bubble through one oscillation cycle, indicating 
the most important phenomena at work. The second 
part of this section deals with unanswered questions in 
SBSL and possible applications of the new insights 
sonoluminescence research has allowed. 


A. An SBSL bubble through its oscillation cycle 


In Fig. 4, we have numbered the various parts of the 
oscillation cycle of a typical SBSL bubble containing ar- 
gon in water [for definiteness, we use an example from 
Storey and Szeri (2000, 2001), with Ro=4,5 um, driven 
at f=26.5 kHz and P,=1.2 atm]. The cycle begins, by 
definition, when the driving pressure begins to dip into 
the negative half of its cycle, and the bubble is therefore 
allowed to expand, The surrounding liquid far from the 
bubble is degassed to some level and maintained at 
room temperature. Consequently, the bubble contains 
about 10!° argon atoms and about 2X10* water mol- 
ecules at the outset. 

(1) Expansion. The bubble expansion is comparatively 
slow and the growth is sustained for almost half a cycle 
(~15 ys). In this phase, the bubble is in both thermal 
and mass transfer equilibrium with the liquid. Because 
of the falling pressure inside the bubble, it gains large 
numbers of water-vapor molecules (evaporating from 
the wall) and also some gas molecules from the liquid. 

(2) Turnaround at maximum radius. The driving pres- 
sure begins to increase again, and the expansion comes 
to a halt. At maximum radius (Ryae~7Ro), the bubble 
contains little more than the initial 10" argon atoms, but 
has collected up to 10” water molecules. 

(3) Rayleigh collapse. As the external pressure in- 
creases, the inertial collapse of the liquid layers around 
the bubble begins. Even with the increased number of 
molecules, the internal pressure is still very low, and the 
collapse proceeds almost exactly like the classical col- 
lapse of an empty cavity treated by Lord Rayleigh 
(1917). The radius decreases quickly (over about 4 us) 
from Rmax to a value comparable to Rọ. During this 
collapse, water vapor recondenses at the wall and the 
argon atoms again become the dominant species inside 
the bubble. 

(4) Decoupling of water vapor. About 50 ns before the 
minimum radius is reached, the time scale of the bubble 
collapse (~R/|R|) becomes smaller than the time scale 
for the diffusion of water vapor. The water vapor still 
left inside the bubble is now trapped until the reexpan- 
sion, Calculations show that about 20% water vapor 
should be mixed with the argon. Up to this moment, the 
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polytropic exponent of the gas mixture has not increased 
significantly above 1, and the temperature has only risen 
to about 500 K. 

(5) Thermal decoupling. Only ~30 ns later, the accel- 
erating bubble wall becomes fast enough that heat can 
no longer escape the bubble. Until reexpansion, the 
bubble is now thermally isolated from the liquid as well, 
and the polytropic exponent y rises quickly to its adia- 
batic value. The latter is determined by the mixture of 
80% Ar and 20% H,O to y~1.6, From now on, the 
temperature increases rapidly. 

(6) Onset of dissociation reactions. Once the tempera- 
ture exceeds roughly 4000 K, water-vapor molecules 
start dissociating into OH and H radicals. While the dis- 
sociation products have lower adiabatic exponents than 
H;O, this endothermic reaction also consumes much en- 
ergy and curbs the temperature rise. At this stage, faint 
molecular band light emission is a possibility. 

(7) Onset of light emission. Despite the temperature- 
limiting influence of water vapor, about 10000 K are 
finally reached in the bubble about 100 ps before maxi- 
mum compression. At this temperature, a small fraction 
of the Ar as well as of the O and H atoms now present 
undergo ionization and release free electrons. The elec- 
trons interact with the ions and neutral atoms, and emis- 
sion of electromagnetic radiation (thermal bremsstrah- 
lung and radiative recombination) begins, which spans 
the spectrum of visible and ultraviolet wavelengths. 
Sonoluminescence is observed. 

(8) Maximum compression. At this point, the gas den- 
sity reaches (almost) solid-state values. The deceleration 
of the bubble wall down to zero speed has begun to 
enhance random shape perturbations (Rayleigh-Taylor 
instability) and leads to massive energy loss through 
acoustic wave emission. The temperature and light emis- 
sion peak, helped by the high densities that prevent fur- 
ther endothermic dissociation reactions. 

(9) Reexpansion. The bubble loses about 90% of its 
energy in the collapse, mostly due to acoustic emission. 
‘The reexpansion is much slower than the collapse. The 
Rayleigh-Taylor instability grows and may overwhelm a 
strongly driven bubble during this stage. Only a small 
increase in radius and decrease in temperature are suf- 
ficient to dramatically reduce the photon absorption co- 
efficient and quench the light emission uniformly for all 
wavelengths, about 100-200 ps after it has begun. Sub- 
sequently, reaction chemistry stops and thermal as well 
as diffusive equilibria are reestablished. 

(10) Afierbounces. The bubble rebounds to a much 
smaller size than the maximum radius before the main 
collapse, and for no parameter combination realized so 
far is there enough energy left to induce sonolumines- 
cence during the afterbounces. The afterbounces pro- 
vide a parametric excitation that can accumulate and 
render the bubble shape unstable. The radial motion is, 
however, damped rapidly until the driving pressure dips 
into its negative cycle once again, and the oscillation 
starts anew. Over the whole cycle, shape perturbations 
may have been enhanced (then the bubble is parametri- 
cally unstable), or a net gain or loss of gas may have 


Brenner, Hilgenfeldt, and Lohse: 


resulted (diffusive instability). In the correct parameter 
range, the bubble is stable with respect to both processes 
and continues to oscillate and emit light in exactly the 
same fashion. 

The above 10 steps hold for SBSL bubbles in water 
with only argon (or any other noble gas) dissolved, If in 
addition molecular gases such as nitrogen and oxygen 
are dissolved, not only water dissociates at step 6, but 
also these gases (at around 7000 K for N and O). The 
reaction products subsequently dissolve in water. 
"Cleaning" the bubble from molecular gases in this fash- 
ion can take thousands of cycles, as shown by Matula 
and Crum (1998). 


B. Unanswered questions 


Possible refinements of this simple picture abound. 
Detailed simulations of the bubble interior have been 
carried out (Moss et al., 1997, 1999; Cheng et al., 1998; 
Storey and Szeri, 1999), indicating that many effects 
could in principle play a role that were not mentioned in 
the above summary. Moss et al, (1997) include the pos- 
sible existence of two temperatures for ions and elec- 
trons, as the collision rate of the atoms may not be suf- 
ficient to thermalize the electrons. They emphasize the 
importance of this effect for a proper modeling of heat 
exchange between the ionized gas and the exterior of 
the bubble. In the same work, shock waves have been 
computed, although Vuong and Szeri (1996) and Lin and 
Szeri (2001) showed that their formation is much de- 
layed due to the entropy profiles inside the bubble. 
While it seems clear that shock waves are not necessary 
to explain SBSL light emission, they could lead to higher 
temperatures in an extremely small central region of the 
bubble. 

Detailed models of the water-vapor chemistry in the 
bubble (Yasui, 1997a; Storey and Szeri, 2000) suffer 
from uncertainties in the reaction rates for even funda- 
mental processes, as the peculiar conditions inside the 
bubble have not been probed in other experimental set- 
ups. Therefore, bubble temperature predictions from 
these reactions carry considerable uncertainty. More 
fundamental work is necessary here, A similar uncer- 
tainty concerns a more quantitative modeling of the 
bound-free contributions to the light emission, as the 
term schemes of the atoms involved could be substan- 
tially altered by the extreme densities. 

Current models of sonoluminescence use simplified 
notions about the final stages of the bubble collapse. 
‘There are two classes of results for heat transfer: simple 
models, which extrapolate the point where the heat 
transfer transitions from isothermal to adiabatic from 
linear theory, and full simulations, which solve for the 
complete heat transfer numerically. The drawback of the 
latter approach is that it is not possible to span the entire 
SBSL parameter space with large numerical simulations, 
while the drawback of the former is that more approxi- 
mations have to be made. 

‘The quest for the ideal liquid for SBSL continues. 
Didenko et al. (2000b) showed that the peculiar combi- 


Rev. Mod. Phys., Vol. 74, No. 2, April 2002 


Single-bubble sonoluminescence 475 


nation of low vapor pressure and benign dissociation 
chemistry favors certain organic liquids. Attempts by 
Baghdassarian et al. (2000) to achieve SBSL in liquid 
noble gases (e.g. for a helium bubble in liquid argon) 
have not achieved higher levels of light intensity than. 
those for the conventional case, One reason for this is 
ihe necessarily much lower temperature at which the 
bubble starts. If the compressive heating at comparable 
expansion ratios starts at, say, 10 K instead of 300 K, the 
peak temperature is automatically penalized by a factor 
of 30, even though the collapse heating itself may be 
more efficient. 

Perhaps other liquids or environments can also be 
used to probe spectral ranges of the light emission out- 
side the optical regime—the exact shape of the ultravio- 
let spectrum for light noble gases is still unknown, due to 
ihe strong absorption of those wavelengths in water 
(and most other liquids). Matula ef al. (2001) have re- 
cently measured the near-infrared part of the spectra of 
SBSL and MBSL, confirming the thermal nature of the 
radiation. 

One of the persistent puzzles in the research on 
MBSL is the origin of spectral lines from metal ions 
dissolved in the water (see Fig. 6). Do the metal ions get 
into the bubble at its asymmetric collapse? Or is the 
water around the bubble somehow playing a role in the 
light-emission process? At present, there are no answers 
to these questions. 

We might also briefly mention here the possible effect. 
of a strong magnetic field on sonoluminescing bubbles. 
While preliminary work by Young et al. (1996) suggested 
that this effect would be very pronounced, later findings 
by the same group revealed that the observed magnetic- 
field effect was mainly on the flask, and not on the 
bubble itself. 


C. Scientific uses and spinoffs 


One question arises naturally whenever SBSL is dis- 
cussed: Is there an application for sonoluminescence? A 
commercial use of the light itself seems unlikely, given 
that only a fraction of ~10~* of the energy in a collaps- 
ing bubble ends up as visual photons. Sonoluminescence 
might provide illumination for certain photographs, e. 
in dentistry where ultrasound-driven devices are use 
(Leighton, 1994). But in most situations, illumination 
can easily be achieved by other means 

Grieser and Ashokkumar (2001) have recently shown 
that sonoluminescence can be used to excite fluorescent 
molecules to emit light themselves, often at much higher 
intensity, a process which those authors dubbed sono- 
photoluminescence. This form of emission due to exter- 
nal photoexcitation changes in the same fashion as the 
sonoluminescence signal, with changes in solution con- 
tents and solution concentration, and thus provides a 
true, amplified representation of the sonoluminescence 
intensity. Being closely associated with chemical reac- 
tions due to cavitation, sonoluminescence light can also 
serve as an indicator for sonochemistry (see below) 
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when light intensities and chemical reaction yields are 
correlated (Beckett and Hua, 2001). 


D. Other applications of bubble dynamics and cavitation 


With the consistent and rather simple picture that has 
now emerged, SBSL could well fulfill its promise as an 
exemplary, well-controlled system of cavitation physics. 
Sonochemistry, in particular, will profit from it, In 
sonochemistry, ultrasound (usually of low frequency) is 
used to enhance, assist, or induce chemical reactions 
that would not occur spontaneously—the ultrasound 
works like a catalyst (Mason and Lorimer, 1988; Suslick, 
1990; Reisse et al., 1999). Collapsing bubbles generate 
temperatures and pressures not easily achieved other- 
wise, and do so repeatedly over very short time scales. It 
is thus possible to induce unique reactions. 

While it has been demonstrated that the presence of 
bubbles and cavitation is essential for sonochemistry, it 
is not clear at all where exactly the chemical reactions 
take place: in the interior of a collapsing bubble, at their 
surface, in a liquid layer around the bubble, or even at 
greater distances, mediated by the diffusion of primary 
reaction products from the bubble interior (Suslick, 
1990). The situation is complicated by the presence of 
many bubbles, inhomogeneities in the sound field, and 
the proximity of vessel walls, which all favor asymmetric 
bubble collapses. The ensuing strong shape deformation 
and/or fragmentation of the bubbles enhance mixing be- 
tween bubble exterior and interior, which makes the 
question of the reaction locus even more difficult to de- 
cide. 

On the other hand, even if a large number of bubbles 
partake in the reaction catalysis, the yields of many 
sonochemical reactions are still woefully low (Reisse 
et al. 1999). Some selected ultrasound-assisted processes 
of chemical technology are now approaching application. 
on an industrial scale, e.g., waste water treatment (Rus- 
sel, 1999). Nevertheless, the reaction efficiency will have 
io be enhanced before a more widespread use of 
sonochemistry is feasible. One approach to achieving 
higher yields is the search for an optimal frequency of 
the driving ultrasound (Beckett and Hua, 2001). 

Just like sonochemistry, a great number of other fields 
rely crucially on the energy-focusing powers of collaps- 
ing bubbles: Materials science makes use not only of the 
high temperatures and pressures, but also of the tremen- 
dous cooling rate of reexpanding bubbles far in excess of 
10” Kis (see Suslick, 1995). These cooling rates are un- 
matched by any other technique, and allow for the fab- 
ication of amorphous metal nanoclusters that prove to 
be highly effective catalysts, as shown by Suslick and 
Casadonte (1987). 

In ultrasonic cleaning, the shear and pressure forces 
around collapsing bubbles are used to rid material sur- 
faces of contaminations (see Leighton, 1994). If this sur- 
face erosion goes too far and affects the material prop- 
erties, we speak of cavitation damage (Leighton, 1994; 
Brennen, 1995). Weninger, Cho, et al. (1997) observed 
that an isolated bubble sitting on a solid surface can 
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even emit light. Ever since Lord Rayleigh's (1917) first 
study of cavitation was motivated by the study of cavi- 
tation damage on ship propellers, this has been a major 
driving force behind research on bubbles. A large indus- 
try is thus concerned with avoiding violent bubble col- 
lapses that can do damage by shear as well as with the 
shock waves they emit, or the thin, fast liquid jets 
ejected from asymmetrically collapsing bubbles (Blake 
et al., 1997; Ohl et al., 1998). Even when the forces are 
not too large, cavitation can simply disrupt laminar fluid 
flow, with serious consequences for applications such as 
fuel transport through valves and tubes or inkjet print- 
ing (Dijksman, 1999), 

In biology, cavitation can occur spontaneously in 
water-transporting vessels in plants, which are often un- 
der high tension (Holbrook et al., 2001). Another nice 
application of collapsing bubbles occurs for various spe- 
cies of “snapping shrimp" (alpheus heterochaelis and 
others), a shrimp whose most distinctive feature is one 
giant claw opposite a normal-sized one. This animal lives 
in large colonies that generate noise so loud that it dis- 
turbs submarine communication. It was believed that 
the noise is emitted when the giant claw rapidly closes 
and its two sides hit each other. However, Versluis et al. 
(2000) showed with the help of high-speed video and 
parallel sound detection that the origin of the noise in 
fact is a collapsing cavitation bubble: When rapidly clos- 
ing the pair of scissors, the shrimp emits a thin water jet 
at such high speeds that a cavitation bubble develops 
When the bubble collapses, sound is emitted in the form 
of a shock wave that stuns or even kills small prey. Very 
recently Lohse et al. (2001) showed that, at bubble col- 
lapse, not only sound but also a short light pulse (about 
10* photons) gets emitted. They called the phenomenon 
"shrimpoluminescence." 

Another example in which high water velocities lead 
to bubble cavitation is the converging fluid flow in a 
Venturi tube analyzed by Peterson and Anderson (1967) 
and later by Weininger et al. (1999). The bubbles cavi- 
tate as they flow through the device, and the resulting 
light emission from this device has many features in 
common with SBSL, for example, enhancement of emis- 
sion using xenon instead of argon gas. 

Some of the most promising applications lie in the 
field of medicine. More than 30 years ago, Gramiak and 
Shah (1968) had already suggested using micrometer- 
sized bubbles as contrast agents for ultrasound diagnos- 
tics. The bubbles are of resonant size for the MHz fre- 
quencies of diagnostic ultrasound, and therefore are 
extremely potent scatterers with cross sections several 
thousand times larger than their geometrical cross sec- 
tion (see Nishi, 1975). When injected intravenously, the 
bubbles allow for brighter images and higher contrast. 
The nonlinearity of bubble dynamics and the violent col- 
lapses help to increase the bubble response and imprint 
a distinctive signature onto the emitted sound, making it 
easier to distinguish bubble echoes from unwanted tis- 
sue reflections (de Jong and Hoff, 1993; Hilgenfeldt, 
Lohse, and Zomack, 1998; Frinking et al., 1999). Yet 
even here the potential for various kinds of “cavitation 
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damage" has to be carefully assessed: mechanical dam- 
age to living tissue (Barnett, 1986), thermal hazard from 
the absorption of high-frequency sound in tissue and 
blood (Wu, 1998; Hilgenfeldt er aL, 2000), or chemical 
hazard from the sonochemical production of radicals in- 
side the body (Barnett, 1986), 

Medicine is also beginning to discover the benefits of 
cavitation damage, using bubbles as vehicles for thera- 
peutic applications as well. It has recently been demon- 
strated by Tachibana et al. (1999) and others that the 
directed delivery of drugs as well as the transfection of 
genes through the wall of living cells is dramatically en- 
hanced in the presence of ultrasound and microbubbles. 
Doubtless, one or several of the above-mentioned pro- 
cesses of cavitation damage are at work here to render 
the cell wall permeable for drugs and genes, but without 
permanently damaging the cell (in experiment, many 
cells appear healthy after treatment, as drugs are effec- 
tive and genes are expressed). The exact mechanism re- 
mains to be uncovered. 


E. Multibubble fields: in search of a theory 


All of the applications mentioned in previous para- 
graphs will be feasible only if many bubbles undergo the 
described processes to focus the energy of the cavitation 
collapse where it is needed. Thus a theoretical descrip- 
tion of sonochemistry or ultrasound diagnostics cannot. 
rely entirely on what we have learned about single 
bubbles studying SBSL. Interactions between bubbles 
and their emitted sound fields will be very important. 
above a certain volume fraction of gas, which may be as 
low as 10-5 (Marsh et al., 1998). These interactions take 
on a plethora of different shapes, such as secondary 
Bjerknes forces (Brennen, 1995; Mettin et al., 1997), 
bubble shadowing (Marsh et al., 1998), collective bubble 
collapses (Brennen, 1995), or collective bubble transla- 
tion (streamers; Akhatov et al., 1996). In addition, mul- 
tibubble applications always have to deal with the inter- 
action of bubbles with boundaries, be they hard (as in 
materials science) or soft (as in biological and medical 
contexts). With the experimental observation and nu- 
merical simulation of jet cavitation in bubble collapses 
near a wall by Vogel and Lauterborn (19882, 1988b), 
Vogel et al. (1989), Tomita and Shima (1990), Blake 
et al. (1997), Philipp and Lauterborn (1997), Lauterborn 
and Ohl (1998), Blake et al. (1998, 1999), Tong et al. 
(1999), Brujan et al. (2001a, 2001b), research in this 
complex area has only just begun. 

Modeling the collapses of interacting bubbles in detail 
is very cumbersome, since the bubble motion is typically 
asymmetric, and the bubble interface can even change 
its topology. The liquid jet can penetrate the bubble that 
generated it, as seen in the famous photograph by Crum 
(1979). Clearly, a tractable theory cannot model each 
individual bubble in a cloud, The ultimate goal is a 
simple, effective-medium theory that describes the re- 
sponse of a bubble field hit by an intense sound wave 
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(continuous or pulsed) in terms of field variables encod- 
ing the deposited energy density and its distribution, 
This goal still seems far ahead. 

So for those with applications in mind, the original 
case of multibubble sonoluminescence from the 1930s 
(Frenzel and Schultes, 1934), deemed unimpressive and 
not very interesting, may in fact be the more challenging 
and rewarding task. Single-bubble sonoluminescence 
has taught us more about the astounding forces at work 
when a bubble collapses, and remains a beautiful and 
unique phenomenon, We now must show that we have 
learned enough to take this research program to the 
next, more general level. 
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